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Abstract 


On the basis of measurements of the size of natural aerosols, an outline is given of 
the importance of aerosols for the budget of impurities in the atmosphere close to the 


surface of the earth. 


New results concerning the chemical composition of natural aerosols are reported. 
These results have been obtained by determining the dispersion of electrons and by direct 
chemical analyses. Most of the nuclei with a radius less than I y seem to consist of am- 
monium sulphate while nuclei originating from sea salt in most cases are larger than I u. 

Furthermore, some consideration is given to the manner by which these substances 
are brought into precipitation. Finally the distribution of these substances in precipitation 


at the surface of the earth is given. 


Einleitung 

Die hier folgenden Ausführungen stellen 
einen Teil der Ergebnisse dar, über die der 
Verfasser auf der vom 13. bis 18. Oktober 
1952 in Stockholm unter Leitung von C. G. 
RossBy stattgefundenen Tagung berichten 
konnte. Zusammen mit einigen früheren Ar- 
beiten bilden sie den Versuch, einen Uber- 
blick zu gewinnen über die Natur unserer 
atmosphärischen Aerosole und deren Rolle 
im Spurenstoff haushalt der Atmosphäre. Neben 
eigenen Messungen werden weitgehend die 
Ergebnisse anderer Autoren herangezogen und 
verarbeitet. Als wesentlicher Gesichtspunkt 
scheint sich dabei zu ergeben, daß die gasför- 
migen Spurenstoffe sowohl hinsichtlich ihrer 
Menge als auch ihrer Rolle bei den Konden- 
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sationsprozessen des Wasserdampfes in der 
Atmosphäre von größerer Bedeutung sind 
als bisher angenommen. Was hier über die 
Beziehungen zwischen solchen gasförmigen 
Spurenstoffen und den Aerosolen beziehungs- 
weise den Kondensationsprodukten des Wasser- 
dampfes gesagt werden kann, ist in vieler 
Hinsicht noch Vermutung. Hier tut sich noch 
ein weites Feld für sorgfältige experimentelle 
Untersuchungen auf, und erst beim Vor- 
liegen wesentlich erweiterter Meßreihen kann 
an eine präzise Formulierung der Vorgänge 
gedacht werden. Die geäußerten Gedanken- 
gänge mögen deshalb unter dem Gesichts- 
punkt evtl. noch korrekturbedürftiger Arbeits- 
hypothesen betrachtet werden. 
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Abb. 1. Mittlere Aerosolverteilungskurven nach Mes- 

sungen in Frankfurt a.M., Feldberg/Ts. (800 m) und 

Zugspitze (3 000). Für Frankfurt und Zugspitze sind 

die Kurven durch Jonenmessungen vervollständigt 
(siehe JUNGE, 1952 c). 


I. Überblick über die Grössenverteilung der 
atmosphärischen Aerosole und die von ihnen 
repräsentierten Substanzmengen 


Für Betrachtungen des Spurenstoff haushaltes 
ist es zunächst wichtig, sich ein Bild über die 
Größenverteilung der Aerosole und damit 
über die von den Aerosolen möglicherweise 
repräsentierten Substanzmengen zu machen. 
In einer vorangegangenen Arbeit (Junge, 
1952c) wurde diese Frage eingehend unter- 
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sucht, die dort gewonnenen Erkenntnisse 
sollen hier — soweit nötig — kurz zusammen- 

efaßt und durch einige neuere Messungen er- 
gänzt werden. 

In Abb. ı sind die mittleren Verteilungs- 
kurven für drei typische Orte in Mitteleuropa 
wiedergegeben, die alles Wesentliche erken- 
nen lassen. In dem gewählten doppeltlogarith- 
mischen Koordinatennetz liegt danach das 
Maximum der Kernverteilung um den Radius 
0,05 x herum. Es läßt sich aus optischen und 
luftelektrischen Daten schließen, daß dieses 
Maximum der Kernverteilung den Bereich 0,01 
u<r< 0,1 u nur selten verläßt. Zu den klei- 
neren Kernen hin erfolgt ein mehr oder minder 
rascher Abfall zu der Kernzahl ©. Für unsere 
vorliegenden Betrachtungen sind diese klein- 
sten Kerne mit r < 0,05 u ohne Bedeutung, 
da ihr Substanzanteil verschwindend klein ist. 
Sehr wesentlich dagegen ist der Verlauf der 
Verteilungskurve der Teilchen mit Radien 
oberhalb 0,1 sw. Hier scheint sich aus noch 
unerklärten Gründen stets eine Verteilung ein- 


zustellen, für die das Produkt +18 konstant 


dlogr 

ist. Das bedeutet, daß für gleiche Abszissenab- 
schnitte dlog r das von den Teilchen einge- 
nommene Volumen von r unabhängig ist. 
Dieses wichtige Verteilungsgesetz gilt nach 
allen direkt oder indirekt gewonnenen Anga- 
ben ziemlich weltweit, wobei allerdings be- 
tont werden muß, daß es nur für Mittelwerte 
angenähert richtig ist, im Einzelfalle jedoch 
merkliche Abweichungen auftreten können. 
Insbesondere bei Kondensationsprozessen des 
Wasserdampfes wird diese Verteilung ganz 
erheblich gestört, da ja für die Kondensation 
gerade die größten Kerne herangezogen wer- 
den. In dem wolkenfreien Raum über der 
Erdoberfläche, wo dieser störende Einfluß 
auf seltene Nebelbildungen beschränkt bleibt, 
stellt der Teil der Größenverteilung mitr > o,1 
ft eine ziemlich invariante Komponente des 
Aerosols dar. Hier sind es die Aitkenkerne 
mitr < 0,1 y, die viel stärkeren Schwankungen 
unterworfen sind, denn sie entstehen bevor- 
zugt in den auf den Kontinenten zahlreichen 
Kernquellen und werden auch durch Koa- 
gulation am schnellsten wieder dezimiert. 


In Tabelle 1 sind die Volumina angegeben, 
die den in Abb. ı dargestellten Mittelkurven 
entsprechen. Wie man erkennt, haben die 
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Tabelle 1. Das vom Aerosol repräsentierte Volumen 

in 10-°cm? pro m? Luft nach den Aerosolvertei- 

lungen von Abbildung 1 bei einer oberen Teilchen- 
grösse von 20 u gemäss Abbildung 3 und 4 


Radiusbereich |%005S | 0,1Sr | 1,0Sr 
SOU Sipe | S204 
Aitken- | große | Riesen- 
Ort — | kerne}| Kerne | kerne 
Frankfurt/M 25 Mes- 
SMREM SERRES ANR 8,7 29 39 
Schwankungsbreite etwa 
(Faktor 1/s—3)....... 3—30 | 10—00 |13— 120 
Zugspitze 20 Messungen. 0,9 3,0 3,0 
Schwankungsbreite etwa 
(Faktor 1/;—5)....... 0,2—5| 0,6—15| 0,6— 15 


Aitkenkerne, obwohl zahlenmäßig weit über- 
wiegend, volumenmäßig keinen bedeutenden 
Anteil am Gesamtaerosol. Die seit vielen Jahr- 
zehnten durchgeführten Aitkenkernzählungen 
geben demnach nur ein schlechtes Bild von 
der Größe der Verunreinigung durch Aerosole. 
In der Tabelle sind neben der Zahl der Mes- 
sungen, aus denen das Mittel gewonnen wurde, 
noch die ungefähren Grenzen angegeben, 
zwischen denen die Werte streuten. 

Selbstverständlich können diese Zahlen vor- 
erst nur einen Überblick über die zu erwar- 
tenden Aerosolvolumina geben. Eingehende 
Messungen in verschiedenen geographischen 
Breiten sowie über Land und Meer sind er- 
forderlich, um sicherere Daten zu gewinnen. 
Insbesondere ist noch wenig über den Zu- 
sammenhang mit der Herkunft der Luftmassen 
bekannt. Da der dunstoptisch bestimmende 
Aerosolanteil im Bereich 0,1 u <r<ıyu 
liegt, läßt sich nur so viel sagen, daß die 
arktischen Luftmassen mit ihren großen Sicht- 
weiten besonders geringe Kernzahlen und 
damit Aerosolvolumina oberhalb 0,1 u er- 
warten lassen. 

Bei den in Tabelle 1 angegebenen Aerosol- 
volumina ist weiter zu bedenken, daß ein 
Teil dieses Volumens von einer wässerigen 
Lösung, ein anderer Teil von unlöslichen Sub- 
stanzen eingenommen wird. Nach unseren 
Untersuchungen über die Konstitution des 
natürlichen Aerosols (JUNGE, 1952 a) sind die 
einzelnen atmosphärischen Aerosolteilchen aus 
mehreren Stoffen zusammengesetzt, löslichen 
sowie schwer- oder unlöslichen (Mischkerne). 


Dies folgt aus der Art des Wachstums der 
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Teilchen mit der Feuchte. Im Mittel über 
dem Kontinent ist das dies Wachstum cha- 
rakterisierende Verhältnis 


Radius bei 95 % Feuchte 


Me Radius bei 40 % Feuchte 


etwa 1,3, während es für ein Lösungströpfchen 
aus reinen hygroskopischen Stoffen, 1,6 bis 
1,8 beträgt. Dem Wert 1,3 entspricht bei einer 
Feuchte von 75% ein Volumenanteil an un- 
löslichen Stoffen von 0,70. Der verbleibende 
Volumenanteil von 0,3 wird von der Lösung 
eingenommen, wobei das Verhältnis von 
gelöster Substanz zur Wassermenge mit der 
stofHichen Beschaffenheit variiert. 

Der Mischkerncharakter wurde bisher im 
wesentlichen durch Messungen über den 
Kontinent belegt. Die wenigen Messungen 
bei maritimer Luftzufuhr auf Föhr, sowie der 
Gang der Sichtweite mit der Feuchte für 
frische Atlantikluft (z. B. für Valentia nach 
WRIGHT, siche JUNGE, 1952 b) zeigen jedoch für 
den Radiusbereich 0,1 bis 1 yu fast die gleichen 
Verhältnisse. Die Persistenz der Aerosole dieses 
Größenbereiches läßt vermuten, daß wir es — 
auch über See — mit stark gealterten konti- 
nentalen Aerosolen zu tun haben, auf jeden 
Fall nicht mit zersprühten Meeressalzen. Das 
letztere wird im Folgenden noch durch andere 
Beobachtungen erhärtet. 

Bei unseren bisherigen Größenverteilungs- 
messungen ergab sich — ähnlich wie bei den 
Zählungen von Woopcock (1949) an den 
maritimen Kochsalzkernen — eine obere 
Kerngröße von 5 bis to w. Noch nicht ver- 
öffentlichte Messungen der Zerstreuungsfunk- 
tion (Vorz) zwangen jedoch zu dem Schluß, 
daß noch größere Teilchen anwesend sein 
müßten. Zur Nachprüfung dieser Frage wur- 
den hydrophob gemachte Glasplättchen in 
besonders hergerichteten, regensicheren Auf- 
fanghäuschen (Abb. 2) ausgelegt. In diesen 
Häuschen wird durch den natürlichen Wind 
stets für genügende Erneuerung der aserosol- 
haltigen Luft gesorgt, aus der die Teilchen 
auf die Platte sedimentieren. Durch die Form- 
gebung der Schutzwände werden die Luft- 
strömungen unmittelbar oberhalb der Platte 
ferngehalten, so daß der Sedimentationsvor- 
gang selbst nicht gestört wird. 

Wird ein solches Plättchen 2 bis 3 Tage ex- 
poniert, so ist das Sediment dicht genug, um 


Abb. 2. Schematischer Querschnitt durch das Auf- 
fanghäuschen für Sedimentationsmessungen. P = Glas- 
platte von 5 x 5 cm. 


bequem ausgezählt zu werden. Ist die Größen- 
verteilung des Sedimentes bekannt, so läßt sie 
sich über die Fallgeschwindigkeit der Teilchen 
leicht auf eine solche der Volumeinheit Luft 
umrechnen. Der große Vorteil dieser einfa- 
chen Methode besteht darin, daß der Anteil 
aller Teilchengrößen mit ungefähr der gleichen 
Genauigkeit bestimmt werden kann. Denn die 
wenigen großen Teilchen erscheinen auf dem 
Präparat infolge ihrer großen Fallgeschwindig- 
keit stark angereichert. Gehorcht z.B. die 
Teilchengrößenverteilung den r°-Gesetz, so 
erscheinen im Sediment die Teilchen in einer 
relativen Häufigkeit, die umgekehrt propor- 
tional ihrem Radius ist, da die Fallgeschwindig- 
keit angenähert mit dem Quadrat des Radius 
wächst. Diese Umrechnung setzt allerdings 
voraus, daß alle Teilchen Kugelgestalt und 
gleiche Dichte besitzen, was natürlich — be- 
sonders bei den größeren staubartigen Teil- 
chen — keineswegs der Fall ist. Nimmt 
man jedoch an, daß die Gestalts- und Dichte- 
abweichungen keinen merklichen Gang mit 
der Teilchengröße besitzen, so eliminiert sich 
dieser Einfluß weitgehend. 

In Abb. 3 ist eine Auswahl der Meßergeb- 
nisse in Frankfurt a.M. wiedergegeben. Die 
dünne Linie stellt die extrapolierte Mittellinie 
der früheren Frankfurter Messungen nach Abb. 
ı dar. Es zeigt sich zunächst, daß die mit 
dieser Methode erhaltenen Ergebnisse sich 
sehr gut den älteren Werten einfügen. Aller- 
dings ist jetzt die obere Grenzgröße auf 15 
bis 20 u Radius heraufgerückt. Es fällt weiter 
auf, wie merkwürdig konstant dieser Anteil 
des Aerosols ist, auch wenn man bedenkt, 
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daß jede Meßreihe ein zeitliches Mittel über 2 
bis 3 Tage darstellt, also die kürzeren Schwan- 
kungen ausgeglichen sind. Wenn auch noch 
weitere Messungen notwendig sind, so darf 
doch nach den bisher vorliegenden Ergeb- 
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Abb. 3. Ergebnisse einiger Frankfurter Sedimentations- 
messungen, 
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nissen vermutet werden, daß das annähernd 
erfüllte r°-Gesetz sowie die obere Grenzgröße 
von 10 bis 20 u auch anderswo über dem 
Kontinent erfüllt ist, allerdings bei anderen 
Absolutzahlen. Eine genaue Betrachtung der 
Kurven zeigt übrigens, daß bei Teilchengrößen 
um 5 4 herum der Verteilungsexponent nur 
etwa 2,5 beträgt, um dann zu den größten 
Teilchen auf über 3 anzusteigen. Es läßt sich 
noch nicht angeben, ob dies reell oder metho- 
disch bedingt ist. Aber die Qualität der Teil- 
chen spricht für das erstere. Denn es zeigte 
sich, daß die ganz großen Teilchen mitr > 5 u 
fast nur aus trockenen Staub bestanden, wäh- 
rend in Übereinstimmung mit unseren früheren 
Untersuchungen die Teilchen mit r < 5 u 
noch einen merklichen Anteil hygroskopischer 
Stoffe enthielten. Wir sind deshalb geneigt, in 
dem Teilchenbereich von etwa r=3 bis 20 u 
die „Staubkomponente” unserer Aerosole zu 
erblicken, die natürlich keine scharfe Grenze 
gegen die mehr hydroskopischen Kerne be- 
sitzt, aber sich möglicherweise durch die er- 
wähnte Abweichung vom r?-Gesetz abhebt. 

Diese hier angewandte Sedimentiermethode 
gestattet auf einfache Weise geographische 
Verteilungen des Staubhaltes unserer Atmos- 
phäre, insbesondere den regionalen Einfluß 
von Industrie- und Siedlungsgebieten zu er- 
fassen. Die Präparate lassen sich so herrichten, 
daß sie transportfähig sind und können durch 
Laien gewonnen werden, um dann an einer 
Zentralstelle ausgezählt zu werden. 

Es ist nun ein eigenartiger Zufall, daß unab- 
hängig von uns Woopcock aus ganz ähn- 
lichen Überlegungen heraus, nämlich um die 
Frage der oberen Grenzgröße der Aerosole 
nachzuprüfen, die gleichen Methoden ver- 
wandte (Woopcock, 1952) und zu ganz ähn- 
lichen Resultaten gelangte. Dabei befaßt sich 
Woopcock ausschließlich mit den Spritzwas- 
seraerosolen über dem Meer bezw. an der 
Küste. Daß es sich tatsächlich bei ihm um 
Meeressalzteilchen handelt, wird durch die 
Kristallform und die mit den Auszählungen in 
Übereinstimmung befindlichen Chloranalysen 
der Aerosolpräparate bewiesen. Er gibt die 
Teilchengrößen in Gramm an. Wir haben 
sie in: Radien umgerechnet, wobei wir eine 
Feuchte von 90 % und damit einen gewissen 
Wasseranteil zugrunde legten,-da die Salzteil- 
chen ja meist (bei Feuchten oberhalb 50%) 
als Tröpfchen vorliegen. Es macht dabei auf 
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die Größe des Radius für den vorliegenden 
Zweck nicht viel aus, ob die Feuchtigkeit 
zwischen so and go %, schwankt. 

In Abb. 4 sind zwei Mittelkurven von Woop- 
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Abb. 4. Mittlere Salzkernverteilungskurven nach Woop- 
COCK, umgerechnet auf Radien bei 90 % relativer Feuchte. 
(Eingeklammerte Werte unsicher.) 


relativer Feuchte 
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cock (1952) in unserer Darstellungsweise 
wiedergegeben. Zum Vergleich ist wieder 
wie in Abb. 3 die Mittelkurve unserer früheren 
Frankfurter Messungen extrapoliert einge- 
zeichnet. Die Kurven von WOODCOCK stim- 
men nun sehr angenähert mit unseren Frank- 
furter Messungen der Abb. 3 überein. Nicht 
nur das Verteilungsgesetz und die oberen 
Grenzgrößen sondern auch die Absolutzahlen 
sind die gleichen wie in Frankfurt, das letztere 
in dieser Genauigkeit wohl ein Zufall. Es zeigt 
sich auch hier wieder, daß die Aerosolver- 
teilungen tatsächlich recht allgemeingültigen 
Gesetzen zu folgen scheinen. 

Die Messungen von WOODCOcK werfen er- 
neut die Frage auf nach der unteren Grenze 
der maritimen Aerosolkomponente. Einen 
wichtigen Hinweis dazu gibt eine Größen- 
verteilung, die Woopcock bei einem tropi- 
schen Wirbelsturm also unter extremen Be- 
dingungen erhielt. Die entsprechend umgerech- 
nete Kurve ist in Abb. 4 eingetragen. Man 
erkennt sofort, daß die natürliche atmosphä- 
rische Meeressalzzerstäubung ganz bevorzugt 
die großen und größten Teilchen erzeugt, wo- 
bei die größten Teilchen in diesem Fall wohl 
nicht miterfaßt wurden. Extrapoliert man diese 
Kurve zu kleineren Radien hin, so schneidet 
sie die Frankfurter Mittelkurve bei etwar=1u. 
Wir sehen darin einen erneuten und unmittel- 
baren Beweis für die von uns vertretene An- 
sicht, daß die maritime Aerosolkomponente 
unterhalb r = 1 u rasch an Bedeutung ver- 


liert, and daß hier andere Aerosole, wahr-: 


scheinlich kontinentalen Ursprungs, auch über 
den Meeren vorzuherrschen beginnen. 

In Tab. 2 legen wir noch einen Auszug aus 
Woopcock’s Messungen mit Angaben über 
das Volumen der Aerosole und des daraus 
berechneten Chlorgehaltes vor. Wie nach 
Abb. 3 und 4 nicht anders zu erwarten, sind 
die Volumenwerte mit denen in Tabelle r in 
guter Übereinstimmung. 

Teilchen von der Größe r= 10 bis 20 u 
haben schon eine erhebliche Fallgeschwindig- 
keit, und es taucht die Frage auf, wie sich deren 
regelmäßige Anwesenheit erklärt. Ganz gleich, 
wie im einzelnen dic vertikale Verteilung des 
Austausches verläuft, die Sedimentation der 
Aerosole geht davon unabhängig über die 
ganze Höhe der Atmosphäre vor sich. Sofern 
kein dauernder Nachschub von Aerosolen von 
der Erdoberfläche her erfolgt, gibt es kein so- 


Tabelle 2. Auszug aus WOODCOCK’s Messungen 
(WOODCOCK, 1952) 


Volumen 
nacı be 72% CI-Gehalt 
Stns pgims FE te in pg/m 
S 25 u 
Flugzeug über See bei 
Miami Bodenwind 
Stärke 3 
Höhe 61m TOW) Ole 7022 RT 
915 m T9 Sue A cae 8,2 
I 370m OTTO > 11,5 
Oberflächenwind Stär- 
ke 3—4 sonst wie oben 
Höhe 61m Pit ST SE CLOS 1,6 
670 m 340107 2.1 
I370m CH 7,9 0,37 
Flugzeug über Flori- 
da ca 100 km land- 
einwärts 
Höhe 152 m T'ON 7/5 NO 738 
1 067 m 122 550,7 7,0 
1523 m 9,7 30,5 6,0 


genanntes Austauschgleichgewicht bezüglich 
der vertikalen Aerosolverteilung. Der Aus- 
tausch wirkt, sofern eine Abnahme der Aero- 
solkonzentration nach oben vorliegt, ledig- 
lich sedimentationsverzögernd. Um sich darü- 
ber ein Bild zu machen, wurde ein Beispiel in 
Abb. 5 durchgerechnet. Zugrunde gelegt 
wurden Teilchen, von 21 u Radius, der Dichte 
ı mit einer Fallgeschwindigkeit von s cm/sec 
und einer Austauschgröße 


z ARS 
ser, (A in gr/cm sec, z = Höhe in m) 


die etwa die mittlere beobachtete Vertikal- 
verteilung wiedergibt. Als Ausgangsverteilung 
wurde die Kurve ı gewählt, wie sie bei An- 
reicherung von Grobaerosolen durch Indu- 
striekamine oder Staubaufwirbelungen zu er- 
warten ist. Ihre Veränderung wurde durch 
schrittweise numerische Integration erhalten. 
Man sieht, daß nur an der oberen Grenze 
in ca. 100 m Höhe ein Aufwärtstransport der 
Teilchen stattfindet und im Endergebnis eine 
ziemliche Gleichverteilung in der Vertikalen 
resultiert. Diese Gleichverteilung ist auch tat- 
sächlich durch Messungen von WOODCock 
(1949, siehe auch Tabelle 2) belegt. Das durch- 
gerechnete Beispiel ist gleichzeitig für andere 
Kombinationen von Austausch-Größe und Teil- 
chenradius gültig, wie in der Abb. 5 angegeben. 
Es zeigt, wie schwer es ist, eine obere Grenz- 
größe überhaupt erklärlich zu machen. Denn 
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sie wird in recht verwickelter Weise durch 
vielerlei Einflüsse gesteuert, wie z. B. die 
Produktionsstärke an grobteiligen Aerosolen, 
die Stärke und Vertikalverteilung des Austau- 
sches und die Sedimentationszeiten. 


II. Untersuchungen über die chemische Zu- 
sammensetzung des atmosphärischen 
Aerosols mittels Elektronenbeugung 


Die Frage nach der Abgrenzung der mariti- 
men Komponente im atmosphärischen Aero- 
sol veranlaßte uns, die chemische Zusammen- 
setzung der Aerosolteilchen zu untersuchen. 
Schon im Herbst 1951 war von uns im Rah- 
men einer elektronenmikroskopischen Rei- 
henuntersuchung des atmosphärischen Aerosols 
die Methode der Elektronenbeugung ange- 
wandt worden, um die chemische Zusammen- 
setzung besonders hervortretender Kompo- 
nenten dieses Aerosols zu bestimmen. (JAcoBI, 
Junce Lippert, 1952). Später wurde diese 
Methode weiter ausgebaut, und es soll an 
dieser Stelle ein kurzer Abriß der bisher ge- 
wonnenen Ergebnisse gebracht werden; zu- 
nächst sind aber einige grundsätzliche Bemer- 
kungen bezüglich der Anwendbarkeit der 
Elektronenbeugungsmethode nötig. 


Methode 


Die Elektronenbeugung hat bei diesem Ver- 
wendungszweck gegenüber den chemischen 
Methoden den Vorteil, daß erstens sehr viel 
kleinere Substanzmengen eine Analyse er- 
möglichen, und daß man zweitens auf diese 
Weise direkt Angaben über die im Aerosol 
vorliegenden chemischen Verbindungen er- 
hält. Im Gegensatz dazu gibt die chemische 
Analyse die An- und Kationen an, aber nicht 
zu welchen Verbindungen diese zusammen- 
treten. Dieser Unterschied kommt daher, 
daß die Elektronenbeugungsmethode die Sub- 
stanzen im festen, kristallinen Zustand analy- 
siert, während bei der chemischen Analyse die 
Substanzen in wässeriger Lösung untersucht 
werden. 

Andererseits hat die Elektronenbeugungs- 
methode — abgesehen von der Unmöglich- 
keit einer quantitativen Analyse — große 
Nachteile. So treten an den abgeschiedenen 
Aerosolteilchen unter dem Einfluß des Va- 
kuums und der Erhitzung im Elektronen- 
strahl zum Teil starke Veränderungen auf, die 
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Abb. 5. Sedimentationsvorgang von Wassertröpfchen 
von 21 4 Radius bei vertikal veränderlichem Austausch. 


es erschweren, genaue Angaben iiber die 
Primärsubstanz zu machen. Als Veränderung 
kommen in Frage: 


1. Chemische Reaktionen (Oxydation und 
Reduktion, Entzug von Kristallwasser), — 


2. Physikalische Veränderungen (Schmelzen 
und Verdampfen). 


Es ist daher z.B. nicht möglich, Beugungsdia- 
gramme von Substanzen aufzunehmen, die bei 
den genannten Bedingungen leicht flüchtig 
sind. Kristalline Substanzen verlieren allmäh- 
lich ihr Kristallwasser und da jedem Hydrat 
zumeist ein anderes Gitter zugeordnet ist, 
erhält man im Elektronenmikroskop je nach 
der Stärke der vorhergegangenen Bestrahlung 
ein anderes DEBYE-SCHERRER-Diagramm. Fer- 
ner muß damit gerechnet werden, daß ein 
Gemisch von mehreren Salzen in einer Lösung 
beim Eindampfen nicht unbedingt in die ein- 
zelnen Komponenten zerfällt, sondern sich 
auch Mischkristalle bilden können, deren 
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ALb. 6. Elektronenmikroskopische Aufnahme von Aero- 
solteilchen (Zugspitze), die typisch sind für die hier 
behandelten Ammoniumsulfatkerne. 


Beugungsdiagramme nur sehr schwer oder 
gar nicht analysiert werden können. 

Es würde an dieser Stelle zu weit führen, 
diese Vorgänge bei den einzelnen im Aerosol 
möglicherweise vorkommenden Substanzen 
zu diskutieren, zumal diese Schwierigkeiten bei 
der praktischen Anwendung der Beugungs- 
methode zur Analyse teilweise umgangen wer- 
den. Denn es ist in den meisten Fällen zu kom- 
pliziert, aus einem Diagramm die Gitter- 
konstanten und damit die vorliegende Sub- 
stanz zu berechnen, insbesondere wenn Über- 
lagerungen von Diagrammen mehrerer Sub- 
stanzen vorliegen. Man vergleicht daher besser 
diese Diagramme unbekannter Substanzen — 
analog zur Methode der Spektralanalyse im 
sichtbaren Licht — mit denjenigen bekannter 
chemischer Verbindungen. Leider fehlt bisher 
jegliche Tabulierung der Spektren einzelner 
Kristalle, so daß eine Identifizierung von Sub- 
stanzen aufgrund von Tabellen nicht möglich 
ist. 

Als Beispiel für eine Umwandlung sei das 
Magnesiumchlorid (MgCl,) angeführt. Dieses 
oxidiert nach Entzug des Kristallwassers bei 
starker Bestrahlung zu Magnesiumoxyd. Wenn 
demnach in einem Diagramm des atmosphä- 
rischen Aerosols MgO nachweisbar ist, so 
kann man daraus schließen, daß primär ein 
Magnesiumsalz vorgelegen hat, das sich bei 
Temperaturen von einigen Hundert Grad C 
zersetzt. Tritt andererseits das Diagramm des 
MgO nicht auf, so beweist dies, daß innerhalb 
der Nachweisgenauigkeit kein MgCl, vorge- 
legen hat. Man ersieht an diesem Beispiel, daß 
für die Analyse eine genaue Kenntnis der 


Umsetzung im Elektronenmikroskop nicht 
unbedingt erforderlich ist. s 

Trotz aller genannten Schwierigkeiten dürfte 
das Elektronenbeugungsdiagramm des at- 
mosphärischen Aerosols schon beim gegen- 
wärtigen Stand unserer Kenntnisse ein wert- 
volles Hilfsmittel sein, da aus einer Änderung 
des Diagramms mit Sicherheit auf eine Ver- 
änderung des Aerosols geschlossen werden 
kann. 


Ergebnisse 


Es soll nun kurz auf die Ergebnisse einge- 
gangen werden, die die Methode der Elektro- 
nenbeugung bei der Analyse des atmosphäri- 
schen Aerosols erbrachte. Die Präparate des 
Aerosols wurden dabei mit Hilfe eines beson- 
ders zudiesem Zweck konstruierten Konimeters 
(JacoBt, JUNGE, LIPPERT, 1952) gewonnen, das 
die Teilchen aus dem Größenbereich 


1,5 X 1072 = ae Om oa 


direkt auf die Blenden des Elektronenmikro- 
skops abscheidet. 

Im Rahmen der erwähnten Reihenunter- 
suchung des Aerosols auf dem Taunus-Obser- 
vatorium bei Frankfurt/M hatte sich heraus- 
gestellt, daß bei maritimen Lufteinbriichen 
eine Komponente des Aerosols sehr stark her- 
vortritt. Bei späteren Versuchen zeigte es sich, 
daß diese Teilchen bei den gleichen Wetter- 
lagen auch im Stadtbezirk von Frankfurt/M 
immer wieder stark hervortraten. Das Beu- 
gungsdiagramm war stets das gleiche (Abb. 8a 
und b). Diese Teilchen erscheinen im Elek- 
tronenmikroskop als runde, wenig kontrast- 
reiche Flecke, die aus einem in Wasser leicht 
löslichem Salz bestehen (Abb. 6). Erst bei 
stärkerer Bestrahlung verdampfen die Teilchen 
teilweise, und es tritt eine körnige Struktur 
auf (Abb. 7 a und b). 

Verständlicherweise haben wir auf die Ana- 
lyse dieser einzigen charakteristischen, von der 
Wetterlage abhängigen Komponente des Aero- 
sols besonderen Wert gelegt. Unsere erste 
Vermutung ging dahin, daß es sich hierbei 
um Salze aus dem Meerwasser handelt, also 
vor allen Dingen um Kochsalz und Magne- 
siumchlorid. Es zeigte sich jedoch, daß das 
Beugungsdiagramm dieser Teilchen keinerlei 
Ahnlichkeiten mit den charakteristischen Beu- 
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7 a 


7b 


Abb. 7. Dicht bestäubte Präparate von natiirlichem Aerosol (Frankfurt) nach schwacher und starker Bestrahlung 
im Elektronenmikroskop. Nach starker Bestrahlung (b) ist der lösliche Bestandteil weitgehend verdampft. 


gungsspektren des Kochsalz und des Mag- 
nesiumoxyds hat (Abb. 8 k). Es ist demnach 
dieser Teil des Aerosols nicht maritimen Ur- 
sprungs. 

Ein Wink, um welche Substanzen es sich 
bei diesen Teilchen handeln könnte, ergaben 
einige Abscheidungen im Frühjahr 1952 auf 
der Zugspitze. Es zeigte sich nämlich, daß diese 
gleiche charakteristische Teilchenart auch in 
diesen Höhen von etwa 3 000 m weitaus den 
Hauptanteil des gesamten, noch vorliegenden 
Aerosols auszumachen scheint, und zwar bei 
fast allen untersuchten Wetterlagen. Dies 
deutete darauf hin, daß es sich hierbei um 
Lösungströpfchen handeln könnte, deren Salz- 
bestandteile aus der Gasphase stammen. Es 
wurden daraufhin Elektronenbeugungsdia- 
gramme von Salzen aufgenommen, die sich aus 
folgenden Bestandteilen zusammensetzen: 
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Kation-en: Na‘, (NH,)' 
Apionen SL NONDFALNOS ELCH EI SO jaar 
(SO,)= 


In Abbildung 8 sind die Beugungsdiagramme 
der wichtigsten der daraus zusammengesetzten 
Salze zusammengestellt. Beim Ammoniumsulfit 
treten zwei Diagramme auf, und zwar setzt 
sich das Sulfit rs schwacher Bestrahlung um, 
und man erhält ein Diagramm, das mit dem 
des Ammoniumsulfats übereinstimmt. 

Aus einem Vergleich mit den Diagrammen 
des natürlichen Aerosols (Abb. 8a hind b) geht 
einwandfrei hervor, daß die Diagramme von 
Ammoniumsulfat und -sulfit innerhalb der 
Meßgenauigkeit mit denen des atmosphäri- 
schen Aerosols genau übereinstimmen. Das 
primäre Spektrum des Ammoniumsulfits konnte 
an zwei Tagen auch beim natürlichen Aerosol 
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Io 
Atmosphärisches Aerosol: 
Taunus-Observat. 1951 
; Frankfurt 1952 l 
Zugspitze April 1952 
(o) I 2 3 4 5cm 


Atmosphärisches Aerosol 
Frankfurt 


Streu- 
b strahlung 


13. und 14. 8. 1952 
(0) | 2 À 3 4 5cm 
x (NH4)2 SO3 nach 
schwacher Bestrahlung 
O | 2 3 4 5cm 
(NH4)5 SO3 nach mittel- 
6 starker Bestrahlung und 
1 (NHq4)2 S04 
0 | 2 3 4 5cm 
6 
(NH4) NOz 
0 | 2 3 4 5cm 


Nag SOz und Nap SO4 


(Q | 2 3 4 Scm 


NaNOo Spektrum fritt 
g | erst nach stärkerer 
ur 


ge à Bestrahlung auf 
(o] I 2 3 4 5 cm 


NHa CI schwache 
Bestrahlung 


Na CI (ausgezogen) 
MgO (gestrichelt) 


4 
! 
| 
! 
! 
! 
I 


(6) | 2 3 4 5 6cm 


Abb. 8. Übersicht über die Elektronenbeugungsdiagramme natürlicher Aerosole und der wichtigsten Ver- 


gleichssalze. 


einwandfrei nachgewiesen werden, während erfaßt werden konnte. Abbildung 9 zeigt das 
in allen anderen Fällen stets nur noch das Elektronenbeugungsdiagramm eines atmos- 
Sulfat — also das Umwandlungsprodukt — phärischen Aerosols, das mit dem Diagramm 
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Abb. 9. Beispiel eines Elektronenbeugungsdiagramms 
von natürlichem Aerosol (Frankfurt), das mit dem pri- 
mären Diagramm des Ammoniumsulfits übereinstimmt. 


in Abbildung 8b übereinstimmt; Abbildung 
10 das Diagramm von Ammoniumsulfat. 

Linien anderer Substanzen treten, soweit 
nachweisbar, nicht auf. Es besteht jedoch die 
Vermutung, daß teilweise noch Natriumni- 
trit vorgelegen haben kann. Bei Natriumnitrit 
trat bei unseren Versuchen ein Beugungsspek- 
trum nämlich erst nach mittelstarker Bestrah- 
lung auf, und es ist möglich, daß diese Be- 
strahlungsintensität bei der Aufnahme der 
Diagramme nicht erreicht wurde. Hier müssen 
noch weitere Versuche abgewartet werden, 
ehe eine einwandfreie Klärung der Frage mög- 
lich ist, ob in der genannten charakteristischen 
Komponente des atmosphärischen Aerosols 
außer Ammoniumsulfit und -sulfat noch Na- 
triumnitrit enthalten ist. 

Aus dem bisher Erreichten geht hervor, daß 
man zu einer vollständigen Analyse des at- 
mosphärischen Aerosols, die auch die flüch- 
tigen Bestandteile einschließt, die chemische 
Mikroanalyse hinzunehmen muß. Aus diesem 
Grunde ist daher auch die chemische Analyse 
des Aerosols in Angriff genommen worden, 
über deren vorläufige Ergebnisse im folgenden 
Abschnitt berichtet wird. Die Elektronen- 
beugungsmethode hat aber wichtige Finger- 
zeige geben können, auf welche Substanzen 
hin besonders analysiert werden muß. Sie 
stellt eine einfache Methode dar, um charakte- 


Tellus V (1953), I 


Abb. 10. 


Elektronenbeugungsdiagramm eines Am- 
moniumsulfatpräparates. 


ristische Veränderungen in der Zusammenset- 
zung des atmosphärischen Aerosols festzustel- 
len. 


IH. Chemische Analyse verschiedener 
Kerngrössenbereiche 


Die im vorangegangenen Abschnitt durch- 
geführten Untersuchungen liessen den Wunsch 
aufkommen, die dort erhaltenen Ergebnisse 
durch direkte chemische Analysen zu ergänzen. 
Dabei veranlasste uns die Frage nach der 
Abgrenzung der maritimen Komponente im 
Aerosol, diese Analysen nach Teilchengrössen 
getrennt vorzunehmen. Die Aufgabe, genü- 
gende Substanzmengen von verschiedenen 
Kerngrößen zu sammeln, lösten wir in fol- 
gender Weise: Bekanntlich wird in den als 
Konimeter bekannten Geräten Luft durch 
einen Spalt hindurchgesaugt, wobei die 
Abmessungen des Spaltes oder die Förderge- 
schwindigkeit der Luft bestimmen, welche Teil- 
chengrößen auf der Fläche hinter dem Spalt 
abgeschieden werden. Dieser Abscheideme- 
chanismus ist im Einzelnen ein sehr kompli- 
zierter hydrodynamischer Vorgang und seine 
genaue Berechnung vorerst kaum durchführ- 
bar. Um sich ein ungefähres Bild von den 
Vorgängen zu machen, kann man folgende 
Abschätzung vornehmen (Abb. ır): 


12 
mm 
A — — 
CULT LTAAL LED 
Abb. 11. Strömungsvorgang am Spalt eines Koni- 


meters. 


Bedeutet 


o = Krümmungsradius ® = Fördermenge 
der Strombahn, 

s = Spaltbreite 

v = Geschwindigkeit 
der Luft in der 
Strombahn 


t = Zeit zum Durch- 
laufen der ge- 
krümmten 
Strombahn 


dann gilt bei vorgegebener Spaltlänge 
sw0; TRpfve sfr; vy Pjs 


also rw s?/®. 


Die Geschwindigkeit der Teilchen unter 
Einwirkung einer Kraft w ist nach STORES 


w 


c= & w/r 


6% u-T 


wo r der Teilchenradius und y die Zähigkeit ist. 
Die Zentrifugalkraft, mit der ein Teilchen aus 
der gekrümmten Strombahn der Luft seitlich 
herausgetrieben wird, ist 


mit d= Dichte der Teilchen. 
Damit wird 


d- 72+ y2 d-r2-@ 
© 0.5 
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das Produkt u- 7 ist die für die Abscheidung an 
der Platte maßgebende Größe, denn sie be- 
stimmt den Betrag, um den ein Teilchen sich 
quer zur Strombahn der Luft bewegt. Wir 
erhalten 


Bei konstantem ur d. h. gleichen Abscheide- 
bedingungen ist also der für die Abscheidung 
charakteristische Radius r, durch den Aus- 


druck bestimmt 
s 
Ve À a 


Die anschließend gebrachten Messungen be- 
stätigen angenähert diese Überlegung. 

In Abb. 12 ist ein schematischer Querschnitt 
durch das verwendete Gerät gegeben. Um in 
erträglicher Zeit möglichst große Substanz- 
mengen zu sammeln, wurde die maximale 
Förderleistung einer vorhandenen Pumpe mit 
etwa 400 cm/sck zugrunde gelegt. Dieser 
Luftstrom wird im zweistufigen Konimeter 
der Reihe nach durch die Spalte 1 und 2 
hindurchgesaugt, wobei Spalt 1 bei einer 
Spaltbreite von 0,85 mm die Teilchen mit r,> 
0,9 u, der Spalt 2 mit 0,2 mm Spaltbreite die 
Teilchen mit 0,2 u < r, < 0,9 u abscheidet. 
Die Spaltlänge betrug 11 mm. 

Zur Eichung des Gerätes wurde ein künst- 
liches Aerosol durch Zerstäuben einer Mag- 
nesiumchloridlösung hergestellt, also ein Aero- 
sol, bei dem die Dichte der Teilchen und ihre 
Gestalt (Kugel) bekannt waren. Die Größe, 
auf die es ankommt, ist der Abscheidefaktor & 
des Spaltes, d. h. der Bruchteil der im Luft- 
strom wirklich vorhandenen Teilchen, der die 
Auffangplatte berührt und dort haften bleibt, 
sowie seine Abhängigkeit vom Teilchenradius. 


Abb. 12. 


Schematischer Querschnitt durch das zwei- 
stufige Konimeter. 
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Abb. 13. Schematischer Verlauf des Abscheidefaktors 

€ (Kurve 1) und der Abscheidung im zweiten Spalt 

des zweistufigen Konimeters (Kurve 2) bei gleichen 
Spaltbreiten. 


Seine Bestimmung würde allerdings die Kennt- 
nis der Ausgangsgrößenverteilung des Aerosols 
voraussetzen, die sehr schwierig zu messen 
ist. Wir umgingen diese Schwierigkeit in fol- 
gender Weise: Haben beide Spalte die gleiche 
Breite, und ist der Verlauf des Abscheide- 
faktors durch Kurve 1 in Abb. 13 gegeben, so 
müßte hinter dem zweiten Spalt eine Teil- 
chengrößenverteilung abgefangen werden, wie 
sie der Kurve 2 entspricht. Wird diese Kurve 


a) 


ausgemessen, so läßt sich umgekehrt daraus 
Kurve ı ermitteln. Mit diesem Verfahren 
wurde die Anhängigkeit der Grenzgröße r, 
bei der Abscheidung in Abhängigkeit von 
der Fördermenge und der Spaltbreite ermit- 
telt. Als Grenzgröße r, wird dabei der Radius 
angesetzt, bei dem der Abscheidefaktor den 
Wert 0,5 besitzt. (Abb. 13.) Das Ergebnis ist 
in Abb. 14 dargestellt. Diese Eichung wurde 
mit einer Pumpe großer Förderleistung durch- 
geführt, um den Kurvenverlauf auch ober- 
halb einer Fördermenge von 400 cm/sec 
zu erfassen. Man erkennt, daß der Verlauf 
der Kurven für die Spaltbreite 1 mm und 
0,5 mm recht gut mit der oben durch 
Abschätzung gewonnenen Formel für r, in 
Einklang steht. Nur für die Spaltbreite 0,25 
mm ist die Grenzgröße stärker zu den kleinen 
Radien hin verschoben als zu erwarten war. 
Wir müssen allerdings bedenken, daß sich die 
Einschnürung des Luftstrahles bei engen 
Spalten prozentual stärker auswirkt als bei 
großen; d.h. man müßte für die tatsächlich 
wirksame Spaltbreite in dem Ausdruck für r, 
nicht den Wert s, sondern z. B. den Wert s—h 
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Abb. 


14. Links Ergebnis der r,-Messung in Abhängigkeit von Spaltbreite und Férdermege. Rechts re als Funk- 


tion der Spaltbreite bei der benutzten Fördermenge von 400 cm?/sck. 
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Abb. 15. 


einsetzen, wo h ein Maß für die Einschnürung 
des Luftstrahles ist. Auf diese Weise gelingt 
es, die festgestellte Abweichung auch quan- 
titativ verständlich zu machen. Im übrigen 
wurde den weiteren Untersuchungen die 
tatsächlich gemessene Grenzgröße zugrunde 
gelegt, so daß der Grad der Annäherung an 
die obige Abschätzung in die Resultate nicht 
eingeht. In Abb. 14 b ist der für die weite- 
ren Betrachtungen wichtige Zusammenhang 
zwischen Spaltbreite und Grenzgröße bei der 
angewandten Fördermenge von 400 cm?/sec 
dargestellt. 

Will man dies Ergebnis für die natürlichen 
Aerosole verwenden, so muß man bedenken, 
daß gemäß der Formel für r, die Dichte der 
Teilchen eingeht. Unsere Abscheidung der 
natürlichen Aerosole erfolgte nun bei rund 
80 % relativer Feuchte, bei denen die meisten 
Teilchen schon einen merklichen Teil Wasser 
enthalten. Wir glauben daher, daß der Unter- 
schied in der Dichte gegenüber den bei 
der Eichung verwendeten Teilchen aus kon- 
zentrierter Magnesiumchloridlösung nicht groß 
ist, zumal die Dichte nur mit der Quadrat- 
wurzel in das Resultat eingeht. 

Die Eichung ergab weiter, daß der Kurven- 
verlauf des Abscheidefaktors um die Grenz- 
größe r, herum nicht immer gleich ist, bezo- 
gen auf die Grenzgröße als Einheit. In erster 
Näherung gilt für die Teilchengrößen r > 
0,7 4 der gestrichelte Verlauf im linken Teil 
der Abb. 15, für die Teilchen mit r < 0,7 u 
der ausgezogene Verlauf; die Trennungslinie 


Abscheidefaktoren der beiden benutzten Spalte. 


der beiden Bereiche ist in Abb. 14a einge- 
zeichnet. Man erhält auf diese Weise für die 
verwendeten Spalte die in Abb. 15 darge- 
stellten Abscheidefaktoren. Die obere Grenze 
bei Spalt ı ist durch die obere Teilchengrenze 
selbst gegeben. 

Ähnlich wie bei optischen Filtergläsern sind 
also die ausgeblendeten Teilchenbereiche nicht 
scharf begrenzt, sondern überschneiden sich 
zum Teil. Die Wahl der Grenze r, bei 0,9 u 
geschah im Hinblick auf die dort zu erwar- 
tende untere Grenze der maritimen Aerosol- 
komponente. Die untere Grenze beim zweiten 
Spalt wurde môglichst tief gelegt, um die 
im zweiten Spalt abgeschiedene Aerosolmasse 
der im ersten Spalt vergleichbar zu machen. 
Das Resultat der Messungen wird zeigen, daß 
die Grenzgröße 0,9 u für die chemische 
Beschaffenheit der Aerosole in der Tat Bedeu- 
tung hat. Bei Verwendung der Kurven für 
den Abscheidefaktor in Abb. 15 ist noch zu 
bedenken, daß Ungleichheiten in der Dichte 
und in der Gestalt der Aerosolteilchen sich 
dahin auswirken, daß die Abscheidegrenzen 
der beiden Kerngrößenbereiche noch etwas 
verwaschener werden, ohne daß sich dabei 
allerdings der Schwerpunkt verlagert. 

Die Abscheidung der Teilchen erfolgte in 
etwa ı mm Abstand hinter dem Spalt auf 
kleinen Plexiglasscheiben. Der Luftstrom wur- 
de vorher auf ca. 80 %, angefeuchtet, denn bei 
dieser Feuchte besitzen praktisch alle Teilchen 
schon eine genügende Wasserhülle, so daß sie 
bei Berührung der Plexiglasoberfläche haften 


MelhismV Gloss) yar 


SPURENSTOFFHAUSHALT DER TROPOSPHÂRE I5 


bleiben. Nur die ganz großen Teilchen mit 
Radien oberhalb etwa 5 u, die meist aus troc- 
kenen Staubteilchen bestehen, werden ver- 
mutlich der Messung zum Teil entgehen, aller- 
dings nur zum Teil, denn sobald sich die erste 
Lage hygroskopischer Teilchen niedergeschla- 
gen hat, werden auch diese Teilchen festge- 
halten. Die Anfeuchtung des Luftstromes ge- 
schah in $ parallel geschalteten Röhren von 
ca. 4 cm Durchmesser und 80 cm Länge, die 
innen mit feuchtgehaltenem Fließpapier aus- 
gekleidet waren. Diese Rohre standen senk- 
recht, um eine Sedimentation von Teilchen 
vor dem Erreichen der Spalte zu vermeiden. 
Bei etwa einstündigem Betrieb der Pumpe 
werden ca. 1,5 m? Luft hindurchgesaugt. Die 
Plexiglasplättchen werden dann herausgenom- 
men, die Staubflecken mit wenigen Tropfen 
destilliertem Wasser betupft, damit die lös- 
lichen Bestandteile von ihnen aufgenommen 
werden können, und diese Lösungsmenge auf 
0,2 cm? aufgefüllt für die Analysen verwandt. 
Es wurden die Mikroanalysenverfahren von 
QUITTMANN und CAUER (1939) angewandt. Da 
die Herstellung eines Präparates für eine Sub- 
stanz ca. eine Stunde währte, so konnte zu- 
nächst nur eine beschränkte Anzahl von Sub- 
stanzen zur Analyse gelangen. Im Hinblick 
auf unsere Fragestellungen wurde folgende 
Auswahl getroffen, bei der noch die wesent- 
lichen Merkmale des Analysenverfahrens sowie 
die absolute Erfassungsgrenze angegeben sind. 


Substanz Analysenverfahren Br fossinergtraze 
ue 
CF Silbernitrat 0,2 
SOL Bariumchlorid 1,0 
NH‘, Nessler-Reagenz 0,1 
NO,’  Reagenz Griess-Ilosvay 0,05 


1 (Iug=1o-®gr) 


Für den Zeitraum vom 28. Juli bis 21. August 
wurde täglich das Aerosol auf diese Substanzen 
hin analysiert. Dabei zeigte sich, daß Nitrit nur 
"ganz vereinzelt die Erfassungsgrenze erreichte, 
also praktisch nicht vorhanden war. Das 
rgebnis für die anderen Substanzen ist in 
Abb. 16 wiedergegeben: 

Zunächst die Chloridwerte. Am unteren 
Rand der Figur ist die Richtung der Haupt- 
luftzufuhr sowie die Art der Luftmasse, ob 
frisch zugeführte Meeresluft oder gealterte 
Kontinentalluft, angegeben. Die Cl-Werte, 
die sich meist unter 1 ug/m® halten, sind im 
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Abb. 16. Zeitliche Darstellung der Aerosolanalysen in , 
Frankfurt. 


allgemeinen für die Riesenkerne (r = 0,9 u) 
größer als für die großen Kerne (r S 0,9 u). 
Besonders an den Tagen mit maritimer Luft- 
zufuhr sinkt der Cl’-Gehalt bei den großen 
Kernen auf o, während er an diesen Tagen bei 
den Riesenkernen sogar etwas größer ist. 
Dies ist in guter Übereinstimmung mit der 
Vorstellung, daß die maritimen Kochsalzkerne 
im wesentlichen auf die Radien oberhalb 1 u 
beschränkt sind. Der Gehalt an Cl’ bei den 
großen Kernen in den gealterten Luftmassen 
dürfte ein Einfluß der Industrien usw. sein. 
Gerade Kernsubstanzen, die das Produkt gas- 
förmiger Reaktionen sind — und das wird vor 
allem in den Industrieabgasen zu erwarten 
sein —, lagern sich bevorzugt an die kleineren 
Kerne an, wie gleich noch näher erläutert 
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Abb. 17. NH',- und SO,”-Gehalt der Kerne in Komponentendarstellung für grosse und für Riesenkerne. 


wird. Daß die Industriegebiete Lieferanten von 
Chloriden sind, ergibt sich aus dem weiter 
unten diskutierten Anstieg des Cl’-Gehaltes 
im Niederschlagswasser in solchen Gebieten. 

Am 13.8. ging kurz vor Beginn der Mes- 
sungen eine schwere Gewitterfront über Frank- 
furt hinweg, die einen starken Dunst offenbar 
von zerspritzten Regentröpfchen zurückließ. 
Der außerordentliche Anstieg aller Substanzen 
und der Teilchenzahl an diesem Morgen ist 
daher als außergewöhnlich zu betrachten. Auch 
am Morgen des 21.8. fielen während der 
Messung einzelne Schauer. 

Interessant sind die NH',- und SO,”-Analysen. 
Bis auf den 13.8. haben die Riesenkerne fast 
kein NH',, auch ihr Gehalt an SO,” ist geringer 
als bei den großen Kernen. Ein eindeutiger 
Gang mit der Luftmasse ist bei diesem ge- 
ringen Material noch nicht zu erkennen; es 
scheint, daß beide Stoffe bei maritimer Luft- 
zufuhr etwas abnehmen. Auffällig ist der 
parallele Gang der beiden Substanzen bei den 
großen Kernen. Um dies besser übersehen zu 
können, wurden die Werte von NH’, und SO,” 
gegeneinander für die großen und Riesen- 
kerne aufgetragen (Abb. 17). Bei den großen 
Kernen wurde jede Analyse durch einen Strich 
gekennzeichnet der die Genauigkeitsgrenzen 
der Analyse kennzeichnet. Man erkennt, daß 
die Mehrzahl der Analysen auf der Linie liegt, 
die der Verbindung (NH,),SO, entspricht, 
zum Teil liegen sie noch auf der Linie, des 
sauren Ammoniumsulfats. Der Gehalt an Cl 
gegen diese Werte ist gering, und man kann 
sagen, daß von den analysierten Substanzen 
das Ammoniumsulfat die vorherrschende ist, 
jedenfalls bei Radien unter 1 w. Es müßte 


noch das Ergebnis von Analysen von Sub- 
stanzkomponenten abgewartet werden, die bei 
unseren Messungen nicht erfaßt wurden, um 
sagen zu können, inwieweit das Ammonium- 
sulfat überhaupt die vorherrschende Substanz 
darstellt. Die im vorigen Abschnitt besproche- 
nen Ergebnisse an Elektronenbeugungsdia- 
grammen legen jedoch den Schluß nahe, daß 
dies tatsächlich der Fall ist. 

Ein ganz anderes Bild zeigen die Riesen- 
kerne (Abb. 17b). Man muß hier annehmen, 
daß das SO,” an andere Kationen gebunden 
oder als freie Säure vorliegt. Für das letztere 
spricht mancherlei Erfahrung, die bei anderen 
Untersuchungen in der Luftchemie gewonnen 
wurde. Nur am 13.8. fällt wieder der Wert 
heraus. 

Der Vergleich der Analysen in allen drei 
Stoffkomponenten Ci’, NH’, und SO,” zeigt 
deutlich, daß bei 1 u Radius eine Grenze in der 
Zusammensetzung der Aerosole liegt. Der 
merkwürdige Unterschied im Gehalt des 
NH, und damit des Ammoniumsulfats zwischen 
den großen und den Riesenernen ließ die 
Vermutung auf kommen,daß das Ammonium- 
sulfat das Produkt von Gasreaktionen ist. Um 
über diese Frage Aufschluß zu erhalten, seien 
die Vorgänge erörtert, die zu einer Anlagerung 
des (NH,).SO, an die Teilchen führen können. 
Man kann sich grundsätzlich folgende Prozesse 
denken, die zu einer Bildung aus den gasförmi- 
gen Komponenten NH, und SO, Anlaß geben: 

1. Absorption von SO, und NH, in dem 
Lösungsanteil der Kerne und dann Bildung 
des (NH,) SO, in der Lösung. In diesem Fall 
würde die gebildete Salzmenge dem Volu- 
men der Kerne proportional sein. 
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2. Adsorption der Gase an der Oberfläche 
der Kerne, z.B. bei geringen Feuchten, wo 
diese weitgehend fest sind und Reaktion an 
der Oberfläche. In diesem Falle wäre die ge- 
bildete Salzmenge der Oberfläche der Teil- 
chen proportional. 

3. Direkte Bildung des (NH,),SO, und 
Anlagerung an die Kerne durch Koagulation. 
In diesem Fall würde, wie die folgende Ab- 
schätzung zeigt, die Salzmenge proportional 
dem Radius der Kerne sein. 

Zur quantitativen Abschätzung dieser drei 
Möglichkeiten nehmen wir an, daß das Ver- 
teilungsgesetz der dritten Potenz gehorcht, 
und daß die obere Grenze der hygroskopischen 
Kerne bei etwa 6 u liegt. Als Grenzradien 
mögen also gelten 


SA (oo ne 6 y 

Spal aie; = << 09 u. 

it BEN? 
dlogr 


volumina 


M 


=< folgt dann für die Aerosol- 
r 


é 6,0 
boat 1 V, = fax: zar a ares - 0,825 
0,9 


ATC 


0,9 
Bpalt 2 V, = fan. AE PONS 
0,2 


Damit wird das Verhältnis der Volumina 


Pi 0,79. Das entsprechende Verhältnis der 


Vi 
Oberflächen ist 2 — 4,1. Für Mischkern aero- 


1 

sole müßte man ein Verhältnis erwarten, das 
zwischen dem der Absorption (0,79) und dem 
der Adsorption (4,1) liegt, sofern diese beiden 
Vorgänge für die Bildung des (NH,),SO, 
maßgebend sind. Denn bei den Mischkernen 
wird nur zum Teil eine feste Oberfläche auf- 
treten, während die Absorption nur in einer 
mehr oder minder dicken Lösungshaut um 
den Kern vor sich gehen kann. 

Im Vergleich dazu zeigen die Messungen, 
wenn man den abnormen Wert vom 13,8. 
ausschaltet, folgende Mittelwerte für den 
NH,-Gehalt: 


große Kerne 1,38 ug/m? 


Riesenkerne 0,116 » Verhältnis 11,9. 
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Der Wert 11,9 liegt also noch erheblich über 
dem der Adsorption. 

Nach Punkt 3 hat man sich den Vorgang so 
zu denken, daf sich aus der Gasphase Molekiile 
oder kleine Molekülaggregate von (NH,),SO, 
bilden, die sich infolge ihrer hohen Beweglich- 
keit in kurzer Zeit an das Aerosol anlagern. 
In der folgenden Tabelle sind die Halbwerts- 
zeiten angegeben, die der Anlagerung in 
Abhängigkeit vom Radius des Moleküls bzw. 
des Molekülkomplexes entsprechen (für die 
mittlere Aerosolverteilung von Frankfurt a.M.): 
Die fraglichen Molekülradien liegen um 
2—3%107® cm herum. 


Tl i 2 in NÉ SMIC OR can 


fij,= 0,71 2,84 11,3 25,5 45,4 72 sec 


Die Verteilung der Anlagerung auf verschiedene 
Teilchengrößen errechnet sich wie folgt: 

Für die Koagulation zwischen zwei Teil- 
chengrößen r, und r, mit den Zahlen n, und n, 
gilt (siehe Junge, 1952 b): 


Al Al 
I + I + 
r ik 
re i (ty + te) rh == 2 Jr 


wo K und A Konstanten und / die mittlere 
freie Weglänge sind. 

Nehmen wir die 1-Teilchen als die Mole- 
külkomplexe an und die 2-Teilchen als die 
Kerne mit Radien größer als 0,1 u, so kann 
vereinfacht werden zu 


dn ke Kl 
ee D Pre 


“fy Ho 


Für ein Teilchen der Art 2 folgt bei n, =1 


d K Ê 
Use Al - 2 = const : fo 
dt 4 112 


d.h. die Zahl der angelagerten Moleküle oder 
Molekülkomplexe ist proportional dem Radius 
der Kerne. Das Verhältnis der koagulierten 
Teilchen für die beiden Radienbereiche ist 
demnach gleichbedeutend mit der Summe der 
Radien dieser Teilchen, und ihr Verhältnis 
ergibt sich zu R,/R, = 19,8. Wenn die Koagu- 
lation also wirksam wäre, müßten die großen 
Kerne 20 mal so viel Ammoniumsulfat ent- 
halten wie die Riesenkerne. 
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Der Vollständigkeit halber sei noch das 
Verhältnis der Teilchenzahlen selbst mit N2/N, 
= 92 angeben. | 

Der gemessene Wert 11,9 für das Verhältnis 
des Ammoniumsulfatgehaltes liegt etwa zwi- 
schen dem der Adsorption und der Koagula- 
tion. Dies würde dafür sprechen, daß zum 
mindesten die Reaktion aus der Gasphase und 
die nachträgliche Koagulation bei dem frag- 
lichen Vorgang eine erhebliche Rolle spielt, 
daß daneben aber evtl. noch Oberflächenreak- 
tionen beteiligt sind. Um diese Frage zu prüfen 
wurde folgender Versuch gemacht, dessen 
Ergebnis allerdings noch durch eingehendere 
Messungen bestätigt werden müßte: In einem 
Raum von ca. 100 m? Inhalt wurden künstlich 
27 ug/m® NH, durch Verdunsten von 0,3 
cm? NH,-Lösung angereichert sowie gleich- 
zeitig durch Verbrennen von Schwefel eine 
Konzentration von 4000 ug/m? SO, erzeugt. 
Der NH',-und SO,’-Gehalt des Aerosols wurde 
vor und eine Stunde nach dieser Anreicherung 
mit Spurengasen gemessen. Es ergab sich: 


Vor der Nach der 
Anreicherung Anreicherung 
NH, SO NH, SO 
große Kerne 0,91 3,4 4,65 17 ug/m? 
Riesenkerne 0,00 1,53 0,00 Pros) 


Der Gehalt der großen Kerne an Ammonium- 
sulfat stieg demnach auf den fünffachen Wert 
an. Der Ammoniumgehalt der Riesenkerne 
dagegen blieb 0, und ihr SO,”-Gehalt stieg auf 
etwa das Doppelte an. Ein Vorversuch, bei 
dem nur NH, angereichert wurde, ergab 
übrigens, daß in diesem Fall keine Erhöhung 
des Ammoniumgehaltes der Kerne zu beobach- 
ten war. Man erkennt, daß das Ergebnis des 
Versuches die bei den natürlichen Aerosolen 
gefundenen Verhältnisse sehr gut erklärt. Wel- 
cher Art die Reaktion im einzelnen ist, läßt 
sich noch nicht sagen. Vermutlich erfolgt erst 
eine Oxydation des SO, und dann die Bildung 
des Ammoniumsulfats mit anschließender An- 
lagerung an die Kerne. Möglicherweise liegt 
aber auch gar keine einheitliche Reaktion vor. 
Wichtig für unsere Betrachtung ist die Tat- 
sache, daß in verhältnismäßig kurzer Zeit eine 
Substanzvermehrung erfolgte. 

Um Anhaltspunkte dafür zu bekommen, 
welchen Mengenanteil die hier analysierten 
Aerosolbestandteile (Cl, NH", und sO, an 
der Gesamtmenge des Aerosols haben, wurde 


vom 11.8. bis 21.8. gleichzeitig die Zahl der 
roßen und der Riesenkerne bestimmt und 
aus ihr das Volumen errechnet. 

Auf der anderen Seite kann man aus der 
Menge der analysierten Substanzen das von 
diesen eingenommene Volumen berechnen. 
Dabei ist zu bedenken, daß die löslichen Stoffe 
auch bei niederen Feuchten meist als Lösungen 
vorliegen, meist in stark übersättigtem Zu- 
stand. Um dies zu berücksichtigen, nehmen 
wir eine gesättigte Lösung von (NH,),SO, 
und für das Chlor eine solche von NACI an. 
Diese gesättigten Lösungen wären bei relati- 
ven Feuchten von 82 % bezw. 72 % imDampf- 
druckgleichgewicht. Das Na als Kation für 
das Cl wurde nur zufällig gewählt, ohne 
damit irgend eine Auffassung auszudrücken. 
Es entsprechen dann 


Ig NZ —> 1,86 cm? gesättigte Lösung 
un 


ee 


Für das (NH,),SO, wurde dabei aus unseren 
Messungen die Summe aus den NH';-undSO,”- 
Komponenten genommen. Mit diesen Annah- 
men über das von den anlysierten Stoffen ein- 
genommene Volumen kann natürlich nur 
eine ungefähre Schätzung vorgenommen wer- 
Das Ergebnis des Mengenvergleiches 
autet: 


re Sn » » » 


Große Kerne Riesenkerne 
aus dem Aerosol be- 
rechnetes mittleres 
Volumen aus ıı 
VVerten ame sono 
aus den analysierten 
Stoffen berechnetes 
Volumen 
Anteil der löslichen 
analysierten Stoffe 
am Gesamtvolumen 0,49 


21,4 x 107% 18,2. x 10-8 cm m? 
10,5 10,706, LO » 


0,57 


Man erkennt, daß nach dieser Abschätzung ein 
erheblicher Teil des Volumens durch die ana- 
lysierten, löslichen Stoffe eingenommen wird. 
Damit gewinnt die Auffassung wiederum eine 
Stütze, daß mit diesen Stoffkomponenten tat- 
sächlich ein wesentlicher Teil der löslichen 
Stofte erfaßt wird. Nach der Mischkernnatur 
der Kerne dürfte der Volumanteil der Lösung 
am Gesamtaerosol sogar nur etwa 0,3 sein, 
doch können die Differenzen durch die Unge- 
nauigkeit der obigen Abschätzung bedingt 
sein. Daß tatsächlich eine erhebliche Menge an 
unlöslicher Substanz sowohl bei den großen 
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Kernen als auch bei den Riesenkernen vor- 
handen war, bestätigte der Augenschein aus- 
nahmslos bei jedem Präparat: Wurde der 
Tropfen destillierten Wassers von dem Prä- 
parat wieder abpipettiert, so blieb ein Staubfleck 
nach, der nicht wesentlich aufgehellt erschien 
gegenüber dem ursprünglichen Staubfleck. 

Zum Abschluß seien hier noch die Mittel- 
werte der Frankfurter Aerosolanalysen zu- 
sammengestellt: 


Große Kerne Riesenkerne 
SO,” 4,38 2,30 g/m? 
INES 1,5% 0,31 » 
GI 0,16 0,56 » 


Bedenken wir, daß einmal das vermessene 
Aerosol am Rande von Frankfurt Teilchen- 
reicher ist als in ungestörten kontinentalen 
Gebieten, daß weiter die analysierten Stoffe 
wohl den wesentlichen Teil der löslichen 
Substanzen darstellen und daß schließich der 
bei der Messung nicht erfaßte Teil der Ait- 
kenkerne keinen bedeutenden Anteil liefern 
kann, so darf man erwarten, daß diese angege- 
benen Mittelwerte tatsächlich die löslichen 
Substanzmengen darstellen, die unter günsti- 
gen Umständen im kontinentalen Aerosol in 
Bodennähe enthalten sind. Wir betonen 
es deshalb, weil die aus luftchemischen Mes- 
sungen anderer Autoren gewonnenen Zahlen 
dieser Stoffe meist wesentlich höher liegen. 
Uns erscheint dieser Widerspruch nur da- 
durch erklärlich, daß ein wesentlicher Sub- 
stanzanteil in gasförmigen Komponenten vor- 
liegt. Denn die bisher geübten Verfahren der 
Luftchemie gestatten keine Trennung zwi- 
schen Aerosol und Gasphase. Für ein tieferes 
Eindringen in die Frage der Luftchemie ist 
deshalb eine Trennung dieser beiden Anteile 
von großer Bedeutung. 


IV. Der Gehalt des Regenwassers an Spuren- 
stoffen und Rückschlüsse auf deren tropo- 
sphärische Quellen 


Aus folgenden Quellen kann ein Nieder- 
schlagselement Spurenstoffe empfangen: 

1. Der ursprüngliche, „aktive“ Kondensa- 
tionskern. 

2. Die Koagulation von übrig gebliebenen, 
„inaktiven“‘ Kondensationskernen mit den 
Wolken- und Niederschlagselementen infolge 
Brown’scher Molekularbewegung. 
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3. Das Abfangen von Riesenkernen (r > 2 u) 
beim Fallen der Niederschlagselemente. 

4. Mit der Wasserdampfkondensation ein- 
hergehende Anreicherung von Spurenstoffen 
aus der Gasphase. 

Wie früher näher ausgeführt (JUNGE, 1952 c), 
sind die Quellen 3 und 4 allgemein überwie- 
gend. Durch die Fallbewegung können nach 
den Rechnungen von LANGMUIR (1948) nur 
Teilchen mit Radien > 2 u, also Riesenkerne, 
abgefangen werden. Zum Punkt 4. konnten 
neuere Beobachtungen an eingedampften 
Schneeflocken mitgeteilt werden (Junge, 
1952c) die eine eindeutige Beziehung der 
Menge des Spurenstoffrückstandes zur Was- 
sermenge der Flocke zeigten. Diese Beobach- 
tungen wurden auf der Zugspitze (3000 m) ge- 
wonnen. Sie können kaum anders gedeutet 
werden, als daß parallel mit der Sublimation 
des Wasserdampfes eine Anreicherung von 
Spurenstoffen in der Flocke erfolgt, zumal 
auch Bruchteile einzelner Schneckristalle sich 
in das Bild einfügten. Auf alle Fälle muß man 
sich von der Vorstellung freimachen, daß im 
Niederschlagswasser überwiegend die ,,akti- 
ven“ Kondensationskerne als Spurenstoffe 
vorherrschen. 

Es erhebt sich nun die Frage nach der zu 
erwartenden Konzentration der Spurenstoffe 
im Niederschlagswasser. Lassen wir zunächst 
Punkt 3, außer Betracht, so läßt sich folgende 
Überlegung anstellen: Die Wolkenluft, in der 
die Bildung der Niederschlagselemente er- 
folgt, habe eine absolute Feuchte von q g/m? 
und einen Spurenstoffgehalt von s ug/m?. Von 
diesen gg Wasserdampf möge der Bruchteil 
o in die Hüssige Phase übergehen und mit ihm 
gleichzeitig von s der Bruchteil 7. Dann ist 
die Konzentration des Niederschlagswassers in 
Milligramm /Liter 


wenn sin wg und q in g/m? eingesetzt werden. 
Sind s und q vorgegeben, so wird k entschei- 
dend durch y bestimmt. Für die Quellen r. 
und 2. kann wohl in erster Näherung x = 1 
angenommen werden, da der Verbrauch von 
„aktiven” und ,inaktiven” Kondensations- 
kernen mit der Intensität des Kondensations- 
prozesses parallel gehen dürfte. Nehmen wir 
unsere Frankfurter Messungen als Grundlage 
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und setzen g=10, so wird das Aerosol etwa 
folgende Konzentration im Niederschlagswasser 
verursachen können: 


SOF 0,67 mg/l 
NEIL Ce » 
CE 0,072 » 


Diese Werte liegen zum Teil erheblich unter 
den gefundenen (siche z. B. Tabelle 5 f.). 

Der Punkt 4. ist fiir unsere Betrachtungen 
von besonderer Wichtigkeit; leider ist über 
ihn nur sehr wenig bekannt. Er spielt unter 
anderem eine sehr wichtige, wenn nicht die 
entscheidende Rolle bei der Kondenskugel- 
methode von CAUER (1941) zur Bestimmung 
der fiir die Kondensationsprozesse wichtigen 
Spurenstoffe der Luft. CAUER setzt aufgrund 
seiner Erfahrungen und im Vergleich mit 
anderen Verfahren y = 1 bei flüssigem Konden- 
sat an der Kugel und x = etwa 1/2 bis 1/3 bei 
Reifbildung. Eine sorgfältige Überprüfung 
dieser Frage steht noch aus, aber in erster 
Näherung werden diese Werte wohl auch für 
die natürlichen Kondensationsprozesse in Wol- 
kenluft anwendbar sein. Dies geht insbesondere 
aus der guten Übereinstimmung hervor, die 
CAUER zwischen der Konzentration der Spu- 
renstoffe an seiner Kondenskugel mit der im 
Regenwasser erhält, wie die unten aufgeführten 
Werte noch zeigen werden. 

Auf die Frage nach der Herkunft der Spu- 
renstoffe im Kondensat an der Kugel werfen 
einige neuere, noch unveröffentlichte Mes- 
sungen einiges Licht, die der Verfasser gemein- 
sam mit Herrn Dr. CAUER durchführen konnte. 
In einem dichtverschlossenen, sauber verglasten 
Holzkasten von 0,4 m? Inhalt, an dessen Ober- 
seite im Innern die Kondenskugel angebracht 
war, wurde der Gang der Aerosolteilchen 
mit der Zeit vermessen. Nach einer Stunde 
ergab sich folgende Abnahme der Teilchenzahl: 


Aatkenkertie. av eee auf ca. 0,70 
Große Kerne ........ auf CA 100$ 
RisSenkemes .... 2 saute cae 0%0. 


Dies ist quantitativ gut erklärbar durch die 
große Beweglichkeit und dadurch rasche 
Diffusion an die Gefäßwände bei den sehr 
kleinen Aitkenkernen, sowie durch die schon 
recht merkliche Sedimentationsgeschwindig- 
keit der Riesenkerne. Dagegen stellen die 
großen Kerne ein sehr persistentes Aerosol 


dar, da bei ihnen sowohl Sedimentation als 
auch Diffusion gering sind. Nun wurde im 
Parallelversuch unter sonst gleichen Bedin- 
gungen die Kondenskugel in der üblichen 
Weise durch eine Eisfüllung gekühlt, und es 
ergab sich innerhalb der Fehlergrenzen der 
gleiche zeitliche Verlauf der Teilchenabnahme. 
Es muß daraus geschlossen werden, daß nur 
geringe Bruchteile des Aerosols im Tau der 
Kugel abgefangen sein konnten. Trotzdem 
zeigte das Kondensat der Kugel Spurenstoff- 
mengen, die durch diesen abgefangenen Aero- 
solanteil mengenmäßig nicht erklärt werden 
konnten. Es mußte also in erheblichem Maße 
eine Abscheidung gasförmiger Komponenten 
stattgefunden haben. 

Für den Punkt 3. ist y gegenstandslos; ein 
fallendes Tröpfchen kann bei genügender 
Fallhöhe beliebig viel Substanz aufnehmen, sie 
braucht nicht in einem festen Verhältnis zur 
Wassermenge des Tröpfchens zu stehen und 
gibt deshalb keinen Anhaltspunkt für die Größe 
s. Das Abfangen der Teilchen läßt sich be- 
rechnen. Ein Tröpfchen vom Radius R (cm) 
durchfalle eine Schicht von der Höhe H (cm). 
Sind in dieser Schicht pro cm? N (r)-Teilchen 
vom Radius r enthalten, so ist das Volumen der 
abgefangenen Teilchen 


V=H:N(n)-R:n2-.E:2, 


wo E der Abscheidefaktor ist und aus den 
Rechnungen von LANGMUIR entnommen wer- 
den kann (LANGMUIR, 1948). Legen wir die 
Acrosolverteilung aus den Frankfurter Mes- 
sungen zugrunde, und nehmen wir die Dichte 
der Aerosolteilchen näherungsweise zu 1 an, 
so erhalten wir für die Konzentration der 
Lösung 
=. CAES) 


r 


wobei über alle Größen der Aerosolteilchen zu 
summieren ist. Setzen wir H= 1000 m, so 
folgt für die Abhängigkeit der Nieder- 
schlagskonzentration k von R 


Tabelle 3 
R — 0,1 0,2 0,4 0,8 1,6 mm 
k RL LL 7,5 mg/l 


Die Werte von k sind von der gleichen Größen- 
ordnung wie die tatsächlich beobachteten, so 
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daß dem Effekt des Abfangens der Teilchen 
wohl ein sehr wesentlicher Anteil an der An- 
reicherung der Spurenstoffe im Niederschlag 
zukommt. Es sei nochmals betont, daß nach 
den Rechnungen von Lancmuir Teilchen mit 
Radien unter 24 grundsätzlich nicht abge- 
fangen werden können. 

Durch dieses Auswaschen tritt in der Schicht 
unterhalb der Wolke eine Verarmung an 
Riesenkernen ein. Die Konzentration des fallen- 
den Niederschlags wird damit eine Funktion 
der Niederschlagsmenge. Beim Durchfallen 
eines Regentropfens durch eine Schicht von 
ı cm Höhe werden 4AN—=—N:E:x:R? 
Teilchen ausgewaschen. Dem Regentropfen 
entspricht eine Niederschlagshöhe von dh = 
it TRE 

3 


wobei h in cm einzusetzen ist. 


Daraus folgt 
3-E-h 
4-R 


N=N,-e 


; N,=Antfangsteilchenzahl. 


Sind h cm Niederschlag gefallen und ist die 
Wolkenhöhe H, so ist die Konzentration an 
gelöster, abgefangener Substanz bei einer 
Dichte ı der Teilchen 


ly —N 
Pa LR BL 


h 


(1) 
Der Schwerpunkt der Teilchen, die abge- 
fangen werden können, liegt bei etwa 6—10 u, 
sagen wir 8 u. Er ist bedingt durch die untere 
Abfanggrenze 2 u und durch die obere Grenz- 
größe der Teilchen bei 15—20 u. Für 8 u- 
Teilchen gewinnen wir für die Größe E und 
die Niederschlagshöhe, bei der die Konzen- 
tration auf die Hälfte abgesunken ist, folgende 


Tabelle 4 
R = 0,1 0,2 0,4 1,0 2,0 mm 
E = 0,55 0,68 0,80 0,83 0,76 
Ei OS RO OCT © 2,5 5,4 mm 
nr Zu O0 27 NEO AO UT 2,4 mm 
2 


Außerdem ist noch die Niederschlagshöhe hy), 
angegeben, bei der die Hälfte der 8 u-Teilchen 
unterhalb der Wolke abgefangen ist. 
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Man erkennt, wie rasch der Prozeß des 
Auswaschens vor sich geht. Bei Regentropfen 
von 0,4 mm Radius ist die Konzentration be- 
reits nach 1 mm Niederschlagshöhe auf die 
Hälfte ihres Anfangswertes, wie er in Tabelle 
3 angegeben wurde, abgesunken. Sinkt die 


Konzentration noch weiter ab, so ist in der 
: : _3Eh 
obigen Formel (1) für k das Glied e 4R 


gegen 1 zu vernachlässigen, und man kann 
schreiben 


k-h= ar: No H (2) 


Physikalisch besagt dies, daß, wenn hı,-Nie- 
derschlag gefallen ist, schon der weitaus größte 
Teil der Riesenkerne ausgewaschen ist, und 
der weiterfallende Niederschlag nur noch eine 
Verdünnung der Lösung zur Folge hat. In 
Abb. 9 sind für Regentropfen mit R = 0,4 
und 2,0 mm die Kurven der Konzentrations- 
abnahme (1) und die obige Näherungsformel 
(2) eingezeichnet. 


In der Formel (2) stellt SAN - H die 


Gesamtmenge der Teilchen mit r > 2 u in der 
Schicht von der Höhe H dar. Die Größen R 
und E treten naturgemäß nicht mehr auf. Sehr 
sorgfältige Messungen der Konzentration in 
Abhängigkeit von der Niederschlagshöhe sind 
unter anderem von ANDERSON (1915) gemacht 
worden. Er findet tatsächlich sehr weitgehend 
die Näherungsformel (2) bestätigt. Die Mindest- 
regenmengen, bei denen er in seinen Kurven- 
darstellungen noch Konzentrationen angibt, 
liegen bei h =2,5—s,0 mm. Das sind Nieder- 
schlagshöhen, bei denen — mittlere Tropfen- 
größen vorausgesetzt — nach Abb. 18 die 
Abweichungen von der Formel (1) gering 
werden. © 
Man könnte in diesen rechnerischen Uber- 
legungen eine hinreichende Erklärung für die 
beobachteten Absolutwerte von k und ihre 
Abhängigkeit von der Niederschlagshöhe er- 
blicken. Indessen darf nicht vergessen werden, 
daß der Gang mit der Niederschlagshöhe etwa 
nach Formel (2) keineswegs auf dem Prozeß 
des Abfangens der großen Teilchen beschränkt 
zu sein braucht. Bei jedem Anreicherungs- 
prozeß der Spurenstoffe im Niederschlag wird 
man zwei Phasen unterscheiden können, näm- 
lich einmal die Anfangsphase, in der der 
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Konpenfration > 


2 # 6 8 10 12 1% 
mm Wiederschlag —+ 


Abb. 18. Berechnete Abnahme der Konzentration der 

Spurenstoffe in: Niederschlag mit der Höhe des Nieder- 

schlags. Gestrichelte Kurve nach Formel (1), ausgezo- 

gene Kurve für Näherungsformel (2), für Trôpfchen- 
radien von 0,4 und 2 mm. 


größte Teil der Spurenstoffe in die Nieder- 
schlagselemente hineingelangt und dann die 
zweite Phase, in der die weiterfallenden Nieder- 
schlagsmengen im wesentlichen nur noch eine 
Verdünnung bewirken, so daß die Näherungs- 
formel (2) in Kraft tritt. Wir glauben deshalb, 
daß neben dem Abfangen der Riesenkerne 
auch der oben unter Punkt 4. genannte Prozeß 
(gasfôrmige Spurenstoffkomponenten) eine 
merkliche Rolle spielt. Eine Trennung der 
beiden betrachteten Effekte aus den vorliegen- 
den Beobachtungen wird vorerst nicht mög- 
lich sein. 

Nach diesen rechnerischen Abschätzungen 
sei anhand von vorliegenden Beobachtungen 
an Niederschlägen versucht, einen Überblick 
über die Quellen der Spurenstoffe auf der Erd- 
oberfläche zu gewinnen. 


1. Der Gehalt an CI. 


Folgende Quellen des Cl-Gehaltes 
Niederschläge sind denkbar. 


a) Spray von Meeressalzen 
b) Gasförmige C/-Komponenten vom Meer 
c) Industrierauche und ähnliches. 


Über Punkt a) sind wir durch die Arbeiten 
von Woopcock unterrichtet, die am Anfan 
erläutert wurden. Punkt b) ist noch hypothe- 
tisch. Hinweise darauf geben vor allem die 
Messungen von CUNNINGHAM (1941) über den 
Gehalt von Seenebeln an Cl. Er findet einen 
deutlichen Zusammenhang mit der Tempera- 
tur des Meereswassers und der Dauer des 
Aufenthaltes, jedoch keinen mit der Wind- 
stärke (Spray-Wirkung). Einen weiteren An- 
halt bieten die Analysen von CAUER (1942) 


der 


> 


der iiber dem Meer und an den Kiisten ein 
anderes Verhältnis von Mg zu Cl’ findet 
als im Meerwasser und glaubt, daß Cl’ aus den 
Spray-Tröpfchen durch chemische Prozesse 
freigemacht wird. Er ist auch der erste, der 
die Bedeutung gasförmiger Cl-Komponenten 
für den Spurenstoff haushalt betont. Weitere 
Hinweise geben die folgenden Tabellen 5 u. 6. 


Tabelle 5. Cl/-Gehalt vom Niederschlag 


Küstenstationen 


Entfernung Dauer der Cl’-Gehalt 


ORFAE von der Beobach- 


Küste tung nes 
Barrahead/Hebriden .. 50 m 2 Jahre 760 mg/l 
ISIE Hes 0010. Holgi Doe «oo 50m T0 WG » 
Brest, abonner yaa - 2 Sy TO) 
SN ERO Lo 0 © 000 500 0 8 » a) ee) 
Plouvenezel ere er 50 » En 132500) 
Lincoln/Neuseeland... 15 » Me) 9,2 » 
Leiden/Holland....... 4 » Te) 10,0 » 
Barbados/Antillen.... ? » 5 » 8,7» 
Georgetown/Brit. 
GUITARE EE EE CCE SDS A) 
Durban/Südostafrika.. ?» 2 » BQ. 2 


Tabelle 6. Cl’/-Gehalt vom Niederschlag 


Inlandstationen 


Or & PR es Cl’ mg/Liter 

TESS EM SR eee eee 2 Jahre 8,7 mg/l 
Dortmund sere aoe eee 2 » 752 » 
BérlineDADIemeeeerEEEe 2 » 52 » 
Berlin Mitten ee 2 » 6,6 » 
DIUNCHÉDELS PEER 2 » 4,3 » 
EardBöhmene ers 1 » 4,3 » 
CUICAS ON PEN RE SE 8,20) 
Cincester/England....... 26 » 352 » 
Garforth/England........ ee) ZA) 
Rothemsted/England .... 24 » 23 
Donnersberg/Böhmen.... 4 » sie hy 
Oberschreiberhau 

Kondenskugel CAUER... Y%  » 0,7 ee) 

Bergstationen 
Of Höhe Cl mg/Liter 

On à sees 3450 m 0,12 mg/l 
PAN AS CMA SENS Odin a ound Buide 2960 » 0,99 =» 
SONDDICR EEE TRE 3106 » 0,51 » 
Bicad ul Midi stews ZI 2877 » O72 100) 
Tatra ae su ae ee Iooo » 0,41 » 


In Tabelle 5 sind z. T. aufgrund der sehr 
umfangreichen Bearbeitung von DRISCHEL 
(1940) eine Reihe von Küstenstationen zu- 
sammengestellt die zeigen, daß in kurzer Ent- 
fernung von der Küste eine rasche Abnahme 
des Cl’-Gehaltes erfolgt, insbesondere wenn 
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sehr stürmische Küsten (Hebriden) einbezogen 
werden. Von einer gewissen Entfernung an 
scheint dann die Konzentration recht konstante 
Werte anzunehmen, die sich auch weiter land- 
einwärts (Tabelle 6) nur noch wenig ändern. 
Dieser verhältnismäßig konstante Gehalt bis 
weit in die Kontinente hinein läßt den Ver- 
dacht aufkommen, daß es sich hierbei nicht 
mehr um ausgewaschene Kochsalzteilchen han- 
deln kann; denn diese nehmen wesentlich 
rascher mit dem Abstand von der Küste ab, 
wie z.B. ein Vergleich der Frankfurter Cl’- 
Analysen des Aerosols mit den Werten von 
Woopcock erkennen läßt. Wir möchten 
annehmen, daß weiter landeinwärts gasför- 
mige Cl-Komponenten in immer stärkerem 
Maße beteiligt sind. 

Sehr gut in die kontinentalen Werte fügen 
sich die Cl-Konzentrationen ein, die CAUER 
(1941) in dem Tauniederschlag an der Kon- 
denskugel in Oberschreiberhau gewann. Es 
sei betont, daß CAUER (QUITTMANN, CAUER, 
1939) diese luftchemische Erfassungsmethode 
gerade deshalb entwickelte, um die bei der 
Niederschlagsbildung eingefangenen Spuren- 
stoffe auch ohne natürlichen Niederschlag der 
Messung zugänglich zu machen. Er findet im 
einzelnen (1941) in Oberschreiberhau im 
Mittel aus 123 Analysen 


mittlere Konzentration an Cl’.. 6,7 mg/l 

mittlerer Wassergehalt der Luft 5,0 g/m? 
Daraus berechneter mittlerer Gehalt 

Berne an Cle RM 82,0 ug/m? 


Diese Umrechnung auf den Gehalt der Luft 
wurde mit einem y von etwa 1/3 vorgenom- 
men. Wenn auch dieser Wert noch unsicher ist, 
so ist der Unterschied gegenüber dem Mittel- 
wert an Cl’ im Frankfurter Aerosol (0,72 ug/m°) 
so groß, daß die Vermutung bekräftigt wird, 
es mögen gasförmige Komponenten bei den 
Cauer’schen Werten beteiligt sein. 
Interessant sind noch die zweijährigen in 
sich homogenen Meßreihen von Essen, Dort- 
mund usw. bis Müncheberg (NEHLS, 1939). Es 
scheint ein geringer Abfall mit wachsendem 
Abstand von der Küste aufzutreten, aber er 
kann durch das Ruhrgebiet vorgetäuscht sein. 
Denn dieser Abfall wird durch den Wert von 
Berlin-Mitte deutlich unterbrochen, welcher 
anzeigt, daß Großstädte bezw. Industrien eine 
nachweisbare, aber alles in allem doch recht 
unwirksame Quelle von Cl darstellen. Dies 
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ist in Übereinstimmung mit MENZL (1948), der 
auf dem Donnersberg in Böhmen ebenfalls 
eine deutliche Cl-Erhôhung fand, wenn der 
Wind zuvor Hüttenwerke passiert hatte. 

Die uns zugänglichen Werte von Berggipfeln 
zeigen eine rasche Abnahme mit der Höhe an 
und zwar auffälligerweise in ungefähr dem 
gleichen Verhältnis wie die Aerosolwerte (Ta- 


belle 6). 


2. SO,’-Gehalt. 


Aus einer großen Zahl von Messungen 
wurde wiederum nach DRISCHEL eine Reihe 
von Stationen ausgesucht, die ihrer Lage nach 
vermutlich wenig durch anthropogene Quellen 


Tabelle 7. SO,’-Gehalt im Niederschlag 


Ungestörte Orte 


(OQ) sf & Dauer SO,” 

MES Iowa USA PE CET EE 1 Jahr 6,5 mg/l 
Eneland nl der eee — Gey 
Westküste......... = AE » 
GIISNON el nse ance eee ig. AS 6,1 » 
Liberty/Indiana, USA ...... ZN Mg) 
Lincoln/Neuseeland......... 4» 2, Se) 
Ouimper/Bretaene m er. — a) 
Rothamsted/England........ 6 » Bear oh) 
‘MennesscewU SA ER cue BT Wo) 
Seenebel (nach CUNNINGHAM) GS 

Oberschreiberhau, Kondens- 

kugel nach (CAUBR er. a 1 » Pr, » 
THAUSLACOELE 4 aoe AE oso oc 30—200 » 
Großstädte br CE 10—25 » 
Kleinstädter er. meer 5—15 » 
Ländliche Gegenden ........ I—8 » 


verunreinigt sind (Tabelle 7). Solche Orte 
sind nicht zahlreich. Die SO”,-Werte schwan- 
ken zwar bis zu einem Faktor 5, aber wenn 
man die recht verschiedene geographische 
Situation berücksichtigt, so ist man doch auch 
hier wieder überrascht, welch einheitliches 
Niveau die Werte einnehmen. Sie fügen sich 
denen von CAUER (1941) und von CUNNING- 
HAM (Seenebel) gut ein. In beiden letzteren 
Fällen wird jedoch ein Auswaschen von Rie- 
senkernen keine Rolle spielen können. Wenn 
man wiederum den daraus von CAUER er- 
rechneten SO,”-Gehalt der Luft von 45,4 ug/m* 
mit unseren Frankfurter Aerosolmessungen 
vergleicht, so kommt man zu dem gleichen 
Schluß, wie bei dem Cl’-Gehalt, nämlich daß 
gasférmige Komponenten (hier vermutlich 
SO,) anwesend sind. 
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Eine der Hauptquellen des atmosphärischen 
SO,’-Gehaltes sind die Industrie- und Sied- 
lungsgebiete, wie einer Zusammenstellung 
solcher Werte am Schluß der Tabelle 7 zu 
entnehmen ist. Es sei noch bemerkt, daß, 
von den eigentlichen Quellgebieten abgesehen, 
die Cl- und SO,’-Gehalte von annähernd 
gleicher Größe sind, scheinbar auch über den 
Meeren. Wir dürfen deshalb vermuten, daß 
neben den anthropogenen auch noch sehr 
wirksame natürliche Quellen für das SO,” vor- 
handen sind. Als solche könnten genannt 
werden die Vulkane und die Zersetzung des 
Eiweißes organischer Substanzen. 


Tabelle 8. NO,- und NH',-Gehalt im Nieder- 


schlag 
Kiistennahe Stationen 

Oi Dauer ne 

mg/l mg/l 
Barbados/Antillen........ 25 Jahre 1,53 0,079 
COVION ER ro ro 2 » 0,310 0,236 
Caracas/Venezuela........ I » 0,300 1,88 
Georgetown/Brit. Guinea.. 25 » 0,314 0,039 
Madras/Vorderindien...... 3 » 0,600 0,164 
Mississippi/Golf v. Mexiko.. ı » 0,330 0,285 
BIebriden@ gr. Re: 6 » 0,146 0,047 
Trondjem/Norwegen...... I » 0,432 0,118 
Lincoln/Neuseeland....... 2 » 0,590 0,161 
Westküste England....... — 0,370 0,100 
Hamamatsu/Japan........ I » — 0,230 
Mittel 0,490 0,30 


Inlandstationen, ungestört 


CairuS Australien Fe 2 Jahre 0,460 0,097 
Debra Dum/Himalaya-Fuß ı » 0,310 0,125 
Flahult/Schweden......... 6 » 0,782 0,555 
GoodwelliUSA Mrs » 0,540 0,342 
Manhattan/Kansas........ 3 » 0,680 0,475 
KRemerziBohmen sea ae — 0,820 0,233 
Oberschreiberbau een. — 0,470 0,755 

Mittel 0,580 0,369 


Kondenskugel CAuEr, Oberschreiber- 


Mittel aus Großstädten und dicht be- 
siedelten Gebieten 


3. NH',- und NO, -Gehalt. 


In Tabelle 8 sind, wieder weitgehend nach 
DriscueL, für den NH',- und NO,-Gehalt 
zunächst küstennahe Stationen, die zum Teil 
auch für die Verhältnisse über See als reprä- 
sentativ gelten können, sowie möglichst unge- 
störte Inlandstationen zusammengestellt. Beide 
Stoffkomponenten zeigen — siehe die Mittel- 
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werte — keine sehr großen Unterschiede 
zwischen Land und See. Eine Erhöhung über 
Land deutet darauf hin, daß hier wohl we- 
sentliche Quellen zu suchen sind. Für NH, 
sind es vor allem die sich zersetzenden orga- 
nischen Substanzen, die natürlich auch auf 
dem Meere, wenn auch nur vereinzelt, anzu- 
treffen sind. Der Unterschied Land-Meer 
würde übrigens ausgeprägter, wenn die beiden 
herausfallenden Werte (NO, Barbardos und 
NH’, Caracas) unberücksichtigt blieben. Im 
Gegensatz zum Sulfat steigen die NH',- und 
NO,-Werte in Industrien und Siedlungsge- 
bieten nur wenig an, so daß die anthropogenen 
Quellen dieser Stoffe im Haushalt der Atmos- 
phäre wohl keine wesentliche Rolle spielen, 
ähnlich wie beim Cl. 

Die von CAUER an der Kondenskugel ge- 
wonnene Konzentration ist zwar etwas hoch, 
fügt sich aber dennoch den anderen Messungen 
ein. Er berechnet daraus einen Gehalt der 
Luft an NH’, von 16 ug/m°, während wir in 
Frankfurt 1,88 ug/m? fanden. Nun ist es 


‘sehr naheliegend, daß in der Luft stets ge- 


ringe Mengen von NH, anwesend sind, wo- 
mit der Unterschied leicht erklärlich würde, 
wenn man annimmt, daß dies gasförmige NH, 
irgendwie im Kondensat gebunden wird. 
Miyake und SUGIURA (1950) geben entspre- 
chende Werte für Tokyo: Im Regenwasser 
1,04 mg/l, in der Luft 17 ug/m®. 


Alls Quellen für das Nitrat führt DRISCHEL an 


1. Synthese der Stickoxyde aus N, und O, 
durch Blitze oder stille Entladungen 


2. photochemische Synthese der Stickoxyde 
3. Photooxydation des Ammoniaks. 


Punkt 3. scheint ihm die wichtigste Quelle zu 
sein. Jedoch sind die Verhältnisse wenig durch- 
sichtig. In einigen Städten ist hoher NH',- 
Gehalt mit niedrigen NO,-Gehalt verbunden 
und umgekehrt, was sich eventuell noch 
durch den langsamen Ablauf der Reaktion 
erklären ließe. Aber auch bei Durchsicht der 
Tabelle ist kein merklicher Parallelgang zwi- 
schen NH‘, und NO,’ vorhanden, wie man 
nach Punkt 3. erwarten sollte. Auch ist eine 
Abnahme mit wachsender geographischer 
Breite kaum angedeutet. Für Punkt 3. dage- 
gen würde sprechen, daß nach einer sieb- 
zchnjährigen Reihe in Ottawa (nach 12) die 
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Amplituden des jährlichen Ganges von NH, 
und NO,, d. h. das Verhältnis 


Sommer — Herbst 
Winter 


folgende, ungefähr gleichen Werte hat: 


NH, 
NO, 


etwa 3 
etwa 2. 


Gleichzeitige Messungen des NO,’-Gehaltes 
im Niederschlag und in der Luft liegen von 
Myaxe und SUGIURA aus Tokyo aufgrund ein- 
jähriger Beobachtungen vor: 


Niederschlag..... 0,27 mg/l 


athe eeu Me 35 „g/m? 
Also auch hier wieder sehr hohe Werte im 
Vergleich mit Aerosolmessungen. 
4. NO,’-Gehalt. 


Die Nitritwerte im Niederschlag liegen we- 
sentlich tiefer als die vom Nitrat. MIJAKE und 
SuGIuRA finden z. B. 


Sommer Winter 
lies ee aE 0,0178 0,0350 mg/l 
DEE Seen RE 0,0085 0,0187 » 
Eamamatsüis. er or 0,0079 0,0II8 » 


Das Minimum im Winter deutet auf Ver- 
brennungsvorgänge als Quelle des Nitrits hin, 
wie sie es auch durch Versuche glaubhaft 
machen konnten. Es besteht darüber hinaus 
die Möglichkeit, daß bei photosynthetischen 
Prozessen das Verhältnis von Nitrit zu Nitrat 
solchen Schwankungen unterworfen ist, daß 
die Nitritwerte im Winter absolut höher 
liegen (ANDERSON, 191$). 

Cauer fand in Oberschreiberhau im Mittel 
im Kondensat der Kugel 0,048 mg/l, woraus er 
einen Gehalt der Luft am NO, zu 0,72 ug/m? 
berechnet. In den Frankfurter Aerosolmes- 
sungen lag der Gehalt an NO, unter der Nach- 
Weisgrenze von 0,05—0,1 ug/m?. 


Mit dieser Auswahl von Substanzen und den 
Angaben über ihre geographische Verteilung 
im Niederschlagswasser wollen wir uns in 
diesem Zusammenhang begnügen. Aus der 
sehr großen Zahl’ solcher Messungen sollte 
hier nur der wesentliche Zug herausgestellt 
werden. In der umfangreichen Literatur sind 
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noch viele interessante Einzelheiten aufge- 
führt, aber es ist schwer, lokales von allgemein 
gültigem zu trennen, wobei „lokal“ in diesem 
Zusammenhang schon eine recht weite Um- 
gebung bedeuten kann. f 

Man könnte weiter neben den wenigen, hier 
angeführten Werten über den Gehalt der 
Spurenstoffe in der Luft noch größeres Mate- 
rial, insbesondere auch von CaAuer (1951), 
beibringen. Im großen und ganzen zeigt sich 
immer das gleiche Bild auch bei Benutzung 
unterschiedlicher Methoden, nämlich daß ein 
wesentlich höherer Gehalt als in unseren Frank- 
furter Aerosolmessungen vorzuliegen scheint. 
Auf die Bedeutung dieses Tatbestandes für 
unsere Auffassung über die Rolle der gasför- 
migen Komponenten haben wir mehrfach 
hingewiesen. Um diesen Punkt nochmals zu 
unterstreichen, bringen wir zum Abschluß 
noch eine Zusammenstellung aller spektro- 
skopisch bisher in der Atmosphäre nachge- 
wiesenen Gase nach Kuiper (1949). Wir haben 
neben dem Volumenanteil noch die Konzen- 
tration in ug/m? angegeben, wobei voraus- 
gesetzt wurde, daß die Atmosphäre vertikal 
gleichmäßig durchmischt ist. Dies wird für 
photochemisch entstandene Spurengase sicher 
nicht zutreffen, so daß hier die troposphäri- 
schen Werte tiefer liegen dürften, wie z. B. 
beim Ozon. Man erkennt, daß N,O immer- 
hin mit fast 1000 ug/m? vertreten ist, desglei- 
chen CH,. Die Tabelle lehrt zweierlei: Ein- 
mal, daß Spurengase in erheblichen Konzen- 
trationen vorliegen, die, wie CH,, bisher 
luftchemisch noch gar nicht erfaßt sind und 
daß zum anderen die spektroskopischen Ver- 
fahren in ihrer Empfindlichkeit nicht an die 
luftchemischen heranreichen, so daß von dieser 
Seite für unsere Problemstellung kaum eine 
Hilfe zu erwarten ist. 


Tabelle 9. Gaszusammensetzung der Atmosphäre 
nach KUIPER 


Höhe in cm 


Gas Dr bei 760 mmHg er a 

14 und 0° C = ugim 
CO; 200 10m 220 cm 593000 
Ne 18 14 16200 
He 5,2 4,3 9300 
CH, 1,5 12 1073 
Kr I 0,8 3773 
N,0 0,5 0,4 990 
H, 0,5 0,4 44,8 
O; 0,4 0,3 887 
Xe 0,08 0,06 467 
H,0 10%—104 
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Cloud Forms of the Jet Stream 


By VINCENT J. SCHAEFER®, General Electric Research Laboratory 


(Manuscript received October 15, 1952) 


Abstract 


A combination of four cloud types including specific forms of cirrus, cirrocumulus and 
altocumulus have been successfully correlated with the reported location of the major 
axis of the jet stream. Recognition of these cloud patterns from the ground may prove 
to be a simple method for locating the high velocity component of the jet. Field obser- 
vations of these spectacular clouds have been made at various parts of the United States. 
More than soo feet of time-lapse motion pictures have been obtained of such clouds and 
show remarkable high velocity, shear, turbulence and waves which typify portions of 
this high velocity stream of air. These cloud structures display the coherency and the high 
rate of motion which would be expected of the jet stream. Preliminary studies have 
indicated that the proximity of the major axis of the jet stream sometimes modifies air 
motions at low altitude and mav be responsible for excessive winds on mountain summits 
and anomalies of forest fire behavior. Detailed studies by time-lapse photography should 
establish more of the physical properties of this important component of the atmosphere. 


During the past year observations have been 
made by the writer which indicate the prob- 
ability that the jet stream may often be locat- 
ed from the ground by visual observation 
of cloud structures. 

This possibility was suggested recently by 
SAWYER (1951) after making a statistical study 
of cloud types reported on synoptic maps 
during passage of the jet stream. 

For nearly ten years a spectacular pattern 
of clouds has been observed by the writer 
at occasional intervals in the vicinity of 
Schenectady, in eastern New York. These 
were observed and sometimes photographed. 

When the high altitude isotach charts 
prepared by WBAN Analysis Center in 
Washington were added to the Facsimile net- 
work in the spring of 1952, the opportunity 
arose for correlating field observations of 
cloud forms with the observed location of 
the jet stream. This paper reports the pre- 
liminary results of a phase of field research 


1 Presented Oct. 17, 1952, at Symposium on Cloud 
Physics at Institute of Meteoroloy of Stockholm 
University. 

2 The field research studies reported in this paper were 
conducted as part of the program of basic research in 
meteorology sponsored by The Munitalp Foundation, Inc. 
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conducted during the past summer by the 
Munitalp Atmospheric Research Project. 
Within a few weeks, it was discovered 
that a combination of specific and rather 
spectacular cloud types invariably coincided 
with the close proximity of the jet stream. 
Figure 1 illustrates the position of the jet 
stream as shown by the WBAN charts during 
several different field observation periods 
during the spring and summer of 1952. Ob- 


OBSERVATION SITE 


CENTER OF JET STREAM AT 500 mb. AS DRAWN ON 
WBAN 1000 MAPS FOR DAYS WHEN TYPICAL CLOUDS 
OCCURRED 
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Fig. 3. Cirrocumulus. 
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Fig. s. Altocumulus billows. 
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servations made on these days and more than 
twenty other periods when the reported 
position of the jet stream coincided with the 
site of field observations have established the 
following cloud types as forms which may be 
utilized to establish the location of the jet 
stream. 

1. High Clouds (Hi and H:)—Cirrus 
streamers of great complexity moving at 
high velocity and showing long tufted stream- 
ers, complex shear lines, and massive whorls. 
Illustrated by Fig. 2 and Fig. 6 a, paced: | 

2. High Clouds (H,)—Cirrocumulus in 
blanket-like masses scattered in a random 
fashion although sometimes in a line showing 
evidence of being at the crests of undulations 
in the stream. Clouds sometimes changing in 
character, shifting rapidly to cirrus streamers 
or showing fine structured waves at very 
high altitude. Some blankets show high order 
Tyndall spectra in green, red, and other colors 
when near the sun. Illustrated by Fig. 3 and 
RE eut 

3. Middle Clouds (M, and M,)—Altocu- 
mulus lenticular wave clouds sometimes in 
great profusion with large lateral dimensions 
in the direction of flow of the stream. Often 
with considerable vertical depth and piling 
up in many layers. Such clouds show little 
apparent relationship to ground topography, 
although they are basically “standing clouds” 
and, therefore, do not exhibit rapid move- 
ment except when snow is shed from them. 
When this occurs, long streamers may extend 
downwind for many miles to emphasize the 
high velocity nature of the air. When near 
the sun, high order Tyndall spectra colors 
are commonplace. Illustrated by Fig. 4 and 
Eig.8 a, bene 

4. Middle Clouds (M, and M;)—Altocu- 
mulus. A billow type cloud which may ex- 
tend from horizon to horizon with the waves 
in parallel bands at right angles to the air 
flow. At times the cloud sheet may appear as 
a relatively thin layer with the units more 
cellular in form. Illustrated in Fig. 5 and Fig. 
Ova, b,c, ds 

In addition to the four basic cloud types 
which seem to be closely related to the main 
axis of the jet stream, the following three 
other characteristics have been noted. 

High Velocity—In all instances the clouds 
show evidence of the high velocity charac- 
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teristic of the jet stream. This is most easily 
seen by using as a reference point some distant 
object. No trouble is experienced in observ- 
ing the high velocity, except with the len- 
ticular forms which do not move rapidly. 
Extremely rapid cloud changes may occur 
with the cloud cover shifting from 1/10 to 
9/10 and back to 1/10 in less than an hour. 
In several instances, the writer has watched a 
storm form along the apparent axis of the jet. 
The storm development proceeded rapidly 
and suggested frontogenesis. In most cases I 
have observed, precipitation under the axis 
of the jet stream has been limited to sporadic 
sprinkles of rain or snow. As a rule, the 
weather under the major axis and to the north 
has been typified by cool crisp air, unlimited 
visibility, mostly blue skies, and persistently 
fine weather. 

Coherency—In every instance when cloud 
forms have been related to the jet stream, 
they have occurred in a coherent pattern 
often extending from horizon to horizon. 
Figure 8 c is an example of a common forma- 
tion where the jet passes in the vicinity of a 
mountain ridge. 

Gustiness—in about half of the cases ob- 
served to date, a peculiar persistent gustiness 
in the air at ground level has been noted. 
Although this feature has not yet been studied 
in detail, such winds seem to be most common 
under conditions having a super adiabatic lapse 
rate. The gustiness at ground level which 
seems to be related to the effect of the jet 
stream may be related to the instability of the 
air. In a number of cases, small stratocumulus 
clouds with tops less than 8,000 feet above the : 
ground have displayed profound modification 
due to the jet. 

A preliminary check of the incidence of 
high winds at the summit of Mt. Washington, 
New Hampshire made by Mr Raymond Fal- 
coner of our laboratory shows a high corre- 
lation between such high winds and the close 
proximity of the main axis of the jet stream 
to the mountain. It is quite likely that the 
world record in recorded wind velocity at a 
ground station measured at Mt. Washington 
of 231 m.ph. and similar excessive winds 
of the past will show a high correlation with 
the presence of the jet stream. 

In this respect the writer has recommended 
to the U.S. Forest Service that anomalies in 
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Fig. 7. Cirrocumulus. 
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Fig. 8. Altocumulus lenticular. 


Fig. 9. Altocumulus billows. 
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forest fire behavior may have a similar rela- 
tionship. 


Detection of the jet stream 


A preliminary “rule of thumb” method 
for using the foregoing cloud types and related 
characteristics in establishing the near prox- 
imity to the observer of the jet stream requires 
the appearance of at least three of the cloud 
types described, along with evidence of high 
velocity movement and coherency. If all four 
cloud types occur, the evidence may be 
considered complete. 

Since May 1, 1952, the writer has been 
accumulating photographic evidence of cloud 
types related to the jet stream. Nearly soo 
feet of time lapse motion pictures have been 
secured thus far in eastern and central New 
York, western Minnesota, western Montana, 
northern Idaho, western Wyoming, and 
southern Colorado. The general characteristics 
of the type clouds in the various regions are 
identical, although local variations may be 
noted. 

The use of clouds to establish the location 
of the main axis of the jet stream is limited, 
of course, to those periods when lower clouds 
do not obscure the sky or when enough 


fo) 
moisture is in the jet to permit cloud forma- 
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tion. It has been noted that the hour after 
sunrise and before sunset is most favorable 
for revealing high thin clouds related to the jet 
at times when very little moisture is in the air. 

An hour by hour study of the cloud forms 
related to the jet stream illustrates many inter- 
esting features of its high velocity, mean- 
dering nature, and “micro” structure. This is 
particularly interesting when the jet is rela- 
tively narrow in width. By employing a vis- 
ual method supplemented by time lapse 
motion pictures over an extensive network, it is 
quite possible that a much better understanding 
of the basic dynamic properties of the jet 
stream may be obtained. 

The relatively dense network of meteoro- 
logical stations throughout the world could 
be easily utilized for pin-pointing the location 
of the jet stream on an hourly basis following 
the cloud identification procedure described 
in this paper. This type of data, supplemented 
by aircraft pilot reports, should lead to a 
better concept of the physical nature of the 
jet stream. 


INEIREISEINIGE 
SAWYER, J., 1951: The distribution of medium and high 


cloud near the jet stream. Meteorological Magazine, 80, 
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The Relative Humidity in Radiation 


Fog 


BANC "BEST, Meteorological Office, London 


(Manuscript received October 27, 1952) 


Abstract 


The relative humidity in radiation fog is governed by the rate at which water condenses 
onto the hygroscopic nuclei and the rate of fall of temperature. An equation describing this is 
developed and from this equation it is deduced that, with the usual rates of fall of temperature 
and salt content of the air, the relative humidity is unlikely to exceed 100.1 per cent. 


I. Introduction 


It is commonly assumed that the relative 
humidity in radiation fog never exceeds 100 
per cent by more than a very small amount, 
of the order of a fraction of one per cent. This 
assumption seems very reasonable but it does 
not appear to have been demonstrated, either 
theoretically or empirically, that the assump- 
tion is valid. Experimental demonstration 
would probably be a matter of considerable 
difficulty and accordingly it has seemed worth 
while examining the question theoretically. 

Fog occurs when sufficient water vapour 
condenses onto hygroscopic nuclei suspended 
in the atmosphere. The balance of evidence 
favours the view that these nuclei consist of 
particles of sea salt. Condensation starts when 
the relative humidity is about 70 per cent, 
the raising of the vapour pressure over the 
droplet by the curvature effect being more 
than balanced by the lowering of the vapour 
pressure because the droplet consists of a 
salt solution. As the relative humidity rises 
the droplet increases in size but while the rela- 
tive humidity remains less than a critical 
value, which exceeds 100 by a very small 


amount, there will be a drop size at which 
the resultant vapour pressure over the drop 
equals the ambient vapour pressure. Having 
reached this size, the drop ceases to grow 
unless the relative humidity continues to rise. 
After the relative humidity exceeds the crit- 
ical value further growth of the drop leads to a 
lowering of the vapour pressure over the drop 
because the effect of curvature in raising the 
vapour pressure is now very small. Accord- 
ingly the drop will now grow indefinitely 
unless the relative humidity falls. The critical 
value of the relative humidity and the corre- 
sponding critical drop size are governed large- 
ly by the mass of the nucleus but are in- 
fluenced to some small extent by the tempera- 
ture. If we assume the nucleus consists of so- 
dium chloride it can be shown that at a tem- 
perature of 283° A the critical relative humid- 
ity and drop radius for a nucleus mass of 
107/? grams are 100.0036 per cent and 21 u 
respectively. 


In the early stages of the growth of the 
drop, while the relative humidity is substan- 
tially below 100 per cent, the rate of growth 
of the drop to its equilibrium size is extremely 
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rapid (see e. g. Best, 1951) and the drop can 
be considered to be in instantaneous equili- 
brium with the environment. As the relative 
humidity approaches and finally exceeds 100 
per cent the rate of growth slows down 
enormously. 

Radiation fog forms when the air tempera- 
ture falls sufficiently. If the amount of water 
—i.e. vapour plus liquid—per unit volume is 
conservative the relative humidity will rise 
as the temperature falls unless the rate of 
condensation of water onto the fog droplets 
is adequate to prevent it. The rate of conden- 
sation will obviously be proportional to the 
number of fog droplets per unit volume and 
will increase as the excess of ambient vapour 
pressure over the vapour pressure at the drop 
surface increases. It will also depend upon the 
properties of each droplet, i.c., upon the 
temperature and the mass of the nucleus. On 
the other hand the rate at which water must 
condense to prevent the relative humidity 
from rising substantially above 100 per cent 
will depend upon the rate of fall of tempera- 
ture. It is the balance between these factors 
which determines the amount of supersatura- 
tion which will occur. 

We shall formulate a differential equation 
for the variation of relative humidity in the 
presence of a number of fog droplets and with a 
steadily falling ambient temperature. This 
equation will be based upon the assumption 
of conservation of total water per unit volume 
and the rate of growth of a drop of salt solu- 
tion in air having a specified relative humi- 
dity. It is not possible to obtain an explicit 
solution to this equation but it can be shown 
that there is an upper bound to the relative 
humidity. The value of this upper bound 
depends upon the mass of the nucleus of the 
drop, the rate of fall of temperature and the 
temperature and will be evaluated for a num- 
ber of typical cases. 


2. The rate of variation of the relative 


humidity in radiation fog 


The rate of growth of a droplet of salt 
solution has been discussed elsewhere (Best, 
1951) and from the equations given there it is 
easily shown that the rate of growth can be 


described by the equations 
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dr 
ne A(H— y) (1) 
M Dkpo 
# ROk-+ADy(dpe/dO) (2) 
y = exp (Plan) — 2, (3) 
Pe Avene (4) 


Po = saturation vapour pressure at tempera- 
ture © (dynes cm7?) 


ON temperature ot ait’ (4) 

r = drop radius (cm) 

T = surface tension (74 dynes cm! for water) 
M = molecular weight of water (18) 

R = gas constant per mol (8.314 + 107 erg 


deg 1 mol- 2 


m = mass of nucleus in one drop (grams) 

À = latent heat of evaporation (cal mol-?) 

D = coeflicient of diffusion of water vapour 
in air (cm? secr}) 

k = thermal conductivity of air (cal cm”! 
secnudesy') 


100 H = relative humidity 


If now we assume that there are N droplets, 
all of radius r, per unit volume and that the 
total amount of water per unit volume (liquid 
plus vapour) remains constant we have 


constant (5) 


ie naar 2 


where go is the saturation vapour density at 
temperature ©. Differentiating equation 5 
with respect to time, substituting for dr/dt 
from equation 1 and writing « for —dO/dt 
(the rate of fall of temperature) we get 


dH ire d0o 


09 re 10 — A NrAUS y) (6) 
This is the differential equation describing 
the variation of H with time. Since pe, r and 
y are also functions of f an explicit solution 
cannot be obtained. We can, however, deter- 
mine a condition that dH/dt shall be negative. 
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3. An upper bound to the relative humidity 


We see from equation (6) that dH/dt will be 
negative, i.c., H will decrease, provided 


rp 
Res (7) 
where 
a d0e 
= = 8 
ax NA dO (8) 


In order to determine whether or not H will 
increase the right hand side of the inequality (7) 
should be evaluated for the value of r appro- 
priate to the existing value of H. Failing a 
solution to equation (6) we cannot do this but 
we can show as follows that the right hand 
side of inequality (7) is bounded and hence 
that there is a value of H which cannot be 
exceeded and may not be reached. 

Consider first the quantity K. K is a function 
of A and A is a function of y. Now physically 
y represents the ratio of the vapour pressure 
over the drop of salt solution to the saturation 
vapour pressure over a plane water surface 
at the same temperature as the drop (see Best, 
1951). As the drop size increases the variation 
of y follows a well known pattern (see e.g. 
SIMPSON, 1941) rising to a maximum value and 
thereafter falling asymptotically to unity. Thus 
y varies with r. The range of y in the condi- 
tions associated with radiation fog is small 
and it has been shown (Best, loc. cit.) that 
with a nucleus as small as 10714 grams the 
maximum value of y is only about 1.0004. 
Now at 283° A, for example, A can be written 
1.5 X 10-8/(0.68 + y) so that we can justifiably 
replace y by unity in this expression for A. 
Doing this we find that A varies from 6 x 1077 
to 12 X 1077 as © varies from 273° A to 293° A. 
Over the same temperature range d0o/dO 
varies from about 3X10-7 to 10x10”. If 
now we give « the large value of 8.3 x ro-4°C/ 
sec (3° C/hour) and N the small value of 10 
we get a value of 5 x 1076 for K. This is much 
smaller than any value of r likely to be sig- 
nificant in fog. 

Now consider the variation of the right 
hand side of inequality (7) as r increases. Since K 
is small compared with r the variation of the 
right hand side of the inequality as r increases 
will have the same general form as the variation 


of y, that is it will increase to a maximum 
value and then decrease again. This maximum 
value will depend upon the value of K and 
will constitute an upper bound to H. We 
can determine this upper bound to H by 
differentiating the right hand side of the in- 
equality with respect to r and equating to 
zero. In this process A is treated as a constant 
for the reasons given above. | 

We may note also that numerical substitu- 
tion shows that P is of the order 2x 10°”. 
This is much smaller than r and accordingly 
after differentiating the right hand side of the 
inequality we may expand the exponential 
term and neglect second and higher powers 
of P/2r. We then find that we are left with the 
equation 


2 y Kr 
a (K4 =| to, oe 


(9) 


This can be solved graphically, the result sub- 
stituted in the right hand side of the in- 
equality (7), and thus an upper bound to H 
determined. 


4. Some results for typical cases 


Of the coefficients in equation (9) the value 
of K is influenced by the rate of fall of tem- 
perature «, the actual temperature and by N. 
Substitution of numerical values shows that 
P can be written 6.4x 10-5/@. Table 1 gives 
some values of 100 Hy, where H, is the 
upper bound to H determined in the way 
described, for various values of the param- 
eters involved. 


Table i. Values of the upper bound to the 
relative humidity 
(©) 
ae er Mm = EMI Sog. 
Temp. | fall of | No. of 
° 4 temp. | drops |Case Case 
°C/hour|per cm?) No. A Ne Hy 
288 I 105) 77. 700.22 72 012.9702,272 
3 Io 2 |101.089 | 2 A |112.795 
3 100 | 3 |100.045 | 3 À |100.447 
I 1000 | 4 |100.003 | 4 A [100.030 
3 1000 | 5 [100.0045) 5 À [100.045 
278 I 10 | 6 |100.169 |6 A |101.725 
3 10 | .7. |100.8174 7.89 T09.7206 
I 1000 | 8 100.003 | 8 A [100.031 
3 1000 9 [100.004 | 9 À |100.042 
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We see from table ı that temperature has 
only a small effect on H, and that to a first 
approximation the order of magnitude of 
the supersaturation represented by Hy (ie. 
Hy—1) varies inversely as m!/?. As would be 
expected H, increases with « and decreases 
as N increases. We are primarily concerned, 
however, with the question whether the re- 
lative humidity represented by 100 H, can 
exceed 100 to a significant extent. Several 
entries in table 1 suggest that this may occur 
but it is necessary to consider the reality of the 
conditions postulated. According to WRIGHT 
(1940) the values selected for m are reasonable 
but we must consider also the implied total 
salt concentration in the air. Wright quotes 
observations which vary from 4x 1071? to 
3 x 1o-§ for this quantity. In table 1 the prod- 
uct Nm varies from 10-18 to 107%. In cases 
IA, 2A, 3A, 6A and 7A we have Nm < 10722. 
Excluding these five improbable cases we 
find that in only two cases, viz. numbers 2 


and 7, does the supersaturation exceed 0.5 
per cent. These two cases occur with a very 
small salt concentration and a rapid fall of 
temperature. In two other cases, numbers 1 
and 6, the supersaturation exceeds 0.1 per cent 
and again the salt concentration is very small. 

Bearing in mind that H, is an upper bound 
to H and that the relative humidity may not 
reach the corresponding value it is fair to 
conclude that with the salt concentrations and 
rates of fall of temperature which are in- 
dicated by observation the relative humidity 
in a radiation fog will seldom exceed 100.1 
per cent and will probably be appreciably 
nearer to 100.0 per cent. If there is any other 
factor operating which tends to remove water 


vapour from the air, e.g. upward diffusion, 


this conclusion will be reinforced. 
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À New Method of Studying the Fine Structure of Air Movements 
in the Free Atmosphere 


By J. WARNER and E. G. BOWEN 


(Manuscript received December 9, 1952) 


Abstract 


A new method of measurement of the fine structure of air movements in the atmosphere is 
described. It depends on the use of radar equipment to follow and track the position of freely 
falling strips of light metal foil which have been distributed by aircraft or balloon in the region 
under investigation. The method is capable of direct measurement of vertical air velocity, 
divergence, and the fine structure of wind in the free atmosphere. The results of some preli- 
minary investigations of the method are described. 


I. Introduction 


There are many aspects of meteorology in 
which observations of air movement are 
necessary. In particular, a knowledge of the 
horizontal and vertical air motions in and 
around clouds is of the utmost importance in a 
study of the physical processes that result in 
cloud formation and the release of precipita- 
tion. Few direct measurements have been 
made, however, of any of the commonly 
experienced atmospheric motions other than 
wind at heights greater than a few tens of 
metres above the earth’s surface. For the most 
part these motions are deduced from other 
observations and the accuracy obtainable is 
not high. 

Many different methods of measurement 
have been used in the past. These have included 
observation, both visual and by radio or radar, 
of balloons; visuil observation of smoke 
either in the form of continuous trails or pufts; 
visual observation of the behaviour of soaring 
birds; observations both in gliders and powered 


! Division of Radiophysics, Commonwealth Scien- 
üfic and Industrial Research Organization, Australia. 


aircraft of the vertical accelerations experienced 
in turbulent air; observations both visual and 
radio that depend on the fluctuations in 
atmospheric refractive index due to turbulence; 
and radar observations of the motion of freely 
falling raindrops in clouds. In addition deduc- 
tions can be made of air movements from 
measurements of pressure and temperature 
fluctuations. 

Since early 1949 the authors have been in- 
vestigating the possible applications of a new 
radar method of studying air movements. In 
certain cases the method used appears to have 
marked advantages over those previously 
employed. This paper describes the method 
and its application to certain problems of 
cloud physics, and outlines a number of other 
ways in which the method may prove useful 


in the future. 
2. The method 


The method depends on the use of radar 
equipment to follow and track the position 
of freely falling strips of light metal foil, 
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known during the last war as “window”. The 
foil is distributed either by aircraft or balloon 
in the area in which the air movement is to 
be investigated. Single bundles of foil may be 
dropped and observed, or the motion of a 
continuous line examined. The most valuable 
technique appears to be that in which a series 
of bundles is dropped in some definite pattern 
and their relative motion observed. The 
minimum separation of the individual bundles 
is made such that they can be resolved by the 
radar equipment used. 

The length of each strip of foil is about half 
the wavelength of the radio signals radiated 
by the radar equipment. There are some tens 
of thousands of strips of foil in each bundle 
which weighs between 10 and 100 grams. The 
echo from a single bundle approximates to 
that from a sphere 1 metre in diameter and it 
can be followed by the radar out to ranges of 
some tens of kilometres. The foil falls at about 
100 cm/sec. in still air and the radar can be 
used to track accurately both its horizontal 
and vertical moticn. The foil used by the 
authors spins about its long axis which remains 
horizontal during its descent. Due to the spin 
it does not fall vertically but spreads out side- 
ways. The resulting dispersion eventually re- 
sults in a weaker radar echo, and the maximum 
period for which a single bundle of foil can 
be followed accurately with the radar at pres- 
ent in use is between ıs and 20 minutes. The 
foil is stiffened by means of paper backing 
which becomes wet in cloud and rain allowing 
the foil to become limp and fall at a greatly 
increased rate. This fact limits its use to in- 
vestigations in clear air or outside clouds. 

The radar equipment used operates on a 
wavelength of 10 cm and is accurate to about 
30 metres in range and 1$ minutes in angle. 
It is possible to track the foil accurately out 
to about so km, though it is visible to over 
100 km. The antenna can be rotated con- 
tinuously about a vertical axis, either manually 
or by means of a motor drive, and the echo is 
displayed on a cathode ray tube by intensity- 
modulating the trace which is rotated in 
synchronism with the acrial. This is the fa- 
miliar Plan Position Indicator type of display. 
The antenna can also be rotated about a 
horizontal axis, thus scanning the sky from the 
horizon to the zenith. Both elevation angle 
and azimuth angle are continuously indicated 


Tellus V (1953), I 


on dials. When the azimuth is held constant, 
and the antenna is scanned in elevation, a 
second intensity-modulated display can be 
used to depict a vertical cross-section through 
the atmosphere. 

Single bundles of foil can be tracked with 
precision radar of the gun-laying type but it 
was found difficult to locate them initially 
and the form of display was unsatisfactory. 
It was concluded that rapid scanning in eleva- 
tion and azimuth and an intensity-modulated 
display were both essential features for satis- 
factorily tracking the foil. 


3. Applications of the method 


Three main applications of the method have 
been studied so far. These are the measurement 
of vertical air velocity, divergence or con- 
vergence of air in the horizontal plane, and 
the fine structure of wind. 


(a) Vertical velocity 

It is clearly possible to track the foil in the 
vertical plane and calculate its rate of fall. The 
vertical air velocity in the region can then be 
deduced provided that the velocity of fall of 
the foil in still air is known. 

The. first experiments in the use of the 
technique were made during an investigation 
in which it was desired to know the vertical 
air velocity in one particular location. Indi- 
vidual bundles of foil were dropped, either 
from aircraft or released from a balloon by a 
temperature or pressure controlled mechanism, 
at a desired altitude and their subsequent 
motion plotted. 
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Fig. 1. Path of a single bundle of foil dropped from 


an aircraft in clear air. The rate of fall is very close 
to the laboratory determined value of 100 cm/sec. 
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Fig. 2. Path of a single bundle of foil dropped from 

an aircraft in the clear space between two extensive 

cloud layers. The rate of rise implies an upward air 
velocity of about 200 cm/sec. 


On most occasions the observed rate of 
fall of the foil in the atmosphere corresponded, 
within the accuracy of the experiment, with 
the value obtained in still air in the laboratory. 
However, on other occasions abnormal results 
were obtained and the foil either rose more or 
less uniformly or followed a complicated path. 
Figure I represents a normal fall and is an 
observation made in clear air of the motion 
of a single bundle of foil dropped from an 
aircraft. Figure 2 shows the result obtained 
when a bundle was dropped in clear air 
between two extensive layers of cloud. The 
rate of rise of the foil corresponds to an upward 
air velocity of some 200 cm/sec. An observa- 
tion in which the foil followed a curved 
path is shown in Figure 3 which represents a 
drop made in clear air from a balloon. 

In some instances a knowledge of the rela- 
tive velocities of adjacent air parcels is of more 
interest than the absolute velocity. Here an 
examination of the motion of a line of foil, can 
give the necessary information. The foil is 
dropped from an aircraft which descends at a 
rate equal to the terminal velocity of the foil 
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Fig. 3. Path of a single bundle of foil released in clear 
air from a balloon by a temperature actuated catch; 
the still air rate of fall is superimposed on the figure. 


so that the latter is initially in one horizontal 
plane. One result obtained when using this 
form of distributing the foil is shown in 
Figure 4, which is a photograph of the cathode 
ray tube display of a radar mounted in an 
aircraft. The aircraft was used to sow a line of 
foil above a cumulus cloud and was then 
flown at right angles to the line and above it. 
The figure shows the line some minutes after 
sowing when it had developed a downward 
kink in the region of the cloud. The cloud 
itself had meanwhile started to sink and even- 
tually disappeared. This form of distributing 
foil may also prove useful in an examination 
of orographic effects and in displaying picto- 
rially the presence of lee waves behind moun- 
tain ranges. 

The above observations relate to low alti- 
tudes, but the method has obvious applica- 
tions to high altitude measurements where 
other methods are difficult. A few observations 
have been made of single bundles of foil 
released from a balloon by a pressure-actuated 
catch at altitudes above 12 km. From this 
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Fig. 4. Photograph of the display tube of an airborne 
radar showing a line of bundles of foil dropped above 
the top of a dissipating cumulus cloud. 
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Fig. 5. Photograph of the radar display showing a 
series of bundles of foil dropped by an aircraft around 
a convective cloud. 


level down to 10 km, where the windspeed 
Was at its maximum value, the wind shear 
was high and could be gauged by the way in 
which the foil spread out in the horizontal 
direction. 


(b) Motion in the horizontal plane 
A powerful method of measuring relative 


air movements in a horizontal plane is to 
observe the motion of a number of individual 
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Fig. 6. Plot of the position of the centres of the 
echoes from six bundles of foil at successive three 


minute intervals. The foil was dropped in clear air from 

an aircraft at a constant height. The heavy lines indi- 

cate the paths of the individual bundles and the light 

closed figures whose area is measured and 
plotted in Figure 7. 
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Fig. 7. The variation of enclosed area. as a function 
of time for the experiment shown in Figure 6. The 
plotted points show both the measured area and the 
probable error of measurement. Clearly there was no 
significant change of area with time. 


bundles of foil dropped from an aircraft in 
the form of a circle or a polygon. To keep the 
figure horizontal the aircraft descends during 
the drop as described before. The change of 
arca of this figure with time then gives a 
measure of the “divergence or convergence at 
that level. 

The appearance of the PPI radar display 
during the course of such an experiment made 
in clear air is shown in Figure 5. The position 
of the centre of the echo from each bundle is 
shown plotted at successive 3-minute intervals 
in Figure 6. The variation of the enclosed area 
as a function of time is shown in Figure 7, 
from which it is seen that no significant 
change occurred. 


30r 


(2250M) 


£ 28r 
œ 
2 
4 
© 
© 26h 
o 
ul 
2 (1750M) 
D 24 23 July 1952 
zZ 
w 

22 | En 1 rl 1 Me SN ib ene en 

5 10 15 
TIME (mins) 


Fig. 8. Results of a drop round a cumulus cloud 1100 
metres deep with its base at 1600 metres. The figure 
shows the variation of enclosed area with time. The 
altitude of the circle of foil is shown by 


approximate 
the figures in italics. 
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Two results which are of marked interest 
were obtained from experiments performed on 
23 and 29 July, 1952. In these cases bundles of 
foil were dropped from an aircraft at intervals 
of about r minute while circling clouds at a 
short and relatively constant distance from 
them. Figures 8 and 9 show the variation as a 
function of time of the area bounded by the 
individual bundles. The air in which the 
foil was dropped on 23 July was clearly 
converging as shown in Figure 8. À second 
experiment was performed shortly after that 
shown in the figure but in clear air away from 
the cloud. No significant change in the en- 


closed area occurred, so that the region of. 


convergence was restricted to the locality of 
the cloud. 

In the case of the experiment on 29 July, 
1952, represented by Figure 9, the area of the 
figure first increased and then decreased to 
about the original value. During the first 
phase the cloud was observed to be dissipating, 
then when convergence started it built up 
again. Subsequently the cloud was observed 
to rain. 

These observations show, in a quantitative 
way, the presence of areas of convergence and 
divergence near convective clouds. It is 
proposed to supplement them later with 
simultaneous observations of the vertical 
motion of the individual bundles of foil. The 
information then available should lead to 
direct confirmation, or otherwise, of theories 
anal: the entrainment of air into cumulus 
cloud. 


(c) Fine structure of wind. 


An examination of the motion of the indi- 
vidual bundles of foil, distributed by an air- 
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Fig. 9. Result of a drop round a cumulus cloud 600 
metres deep with its base at 900 metres. 
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Fig. 10. The fine structure of wind in the region sur- 

rounding a convective cloud. The mean wind vector 

is shown and the difference between actual and mean 

wind at the position of each bundle of foil averaged 
over a 3 minute interval. 


craft in lines or other patterns, enables measure- 
ments to be made of the wind vector at any 
desired point. Thus the wind structure through- 
out an area of some hundreds of square kilo- 
metres can be readily measured. The minimum 
separation of individual bundles of foil such 
that they can be resolved on the radar is 
about 1 km, and the maximum distance 
beyond which accurate measurement is im- 
possible due to other equipment limitations 
is about 100 km, thus the scale of wind meas- 
urements possible through the use of this 
technique is from 1 to 100 km. 

The fine structure of wind in the region 
surrounding a convective cloud, as deter- 
mined by the use of this technique, is shown 
in Figure 10. The wind vector at the position 
of each bundle of foil was first determined 
and then the mean wind vector calculated. The 
vector departure from the mean wind at each 
position was then calculated and plotted. 
Evidence of convergence is again indicated. 


4. Conclusion 


The preliminary investigations that have 
been made of air movements, by tracking 
with radar bundles of metal foil distributed 
by aircraft or balloons, suggest that this new 
technique will prove a powerful tool for 
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research in atmospheric physics. The method 
is capable of detecting and measuring vertical 
and horizontal air movements over con- 
siderable areas, and appears to be very useful 
in studies of convective clouds. It is likely 
that the method will also be of use in the 
investigation of the waves formed in the lee 
of mountains or of atmospheric oscillations. 
À study of the rate of dispersal of the foil, by 
an examination of the dimensions of the echo, 
could also prove useful in studies of small 
scale turbulence, particularly in the upper 
atmosphere where other methods of measure- 
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ment are difficult. This type of investigation 
may, however, have to await the availability 
of a different type of echoing material which 


falls vertically without marked dispersal in 
still air. 
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Abstract 


The following paper is the first of a series of two relating to the problem of internal oscilla- 
tions of a fluid in a gravity field with vertical gradients of density and velocity. The theoretical 
analysis in this paper will be supplemented by a report on an experimental investigation along 
the same lines to be published in a future issue of this journal. ‘ 

The exact, steady-state equations of motion and continuity of a perfect liquid moving two- 
dimensionally, with an arbitrary vertical distribution of density and velocity, are integrated 
once to yield a second-order differential equation. This equation is examined with regard to 
uniqueness and stability of the motion. A criterion is developed giving a sufficient condition 
for the motion to be uniquely determined by the configuration of the topography over which 
the fluid moves. It appears, further, that the condition of uniqueness is also a condition that a 
certain integrated quantity, called the kinetic potential of the motion, be a maximum. The 
suggestion is offered that this may correspond to a form of fluid instability. 

A detailed study is made in the special cases of a uniform basic velocity, and a certain type of 

‘shearing flow. In either case, it is shown that an internal Froude number of about 1/3 divides 
the motion into two states, one of which is called supercritical, the other subcritical. From several 
viewpoints, these regimes are analogous to the corresponding states of flow of water in a channel. 
In the subcritical state the flow is in the form of standing wave patterns. When flowing supercrit- 


ically, conditions may be favorable for the formation of internal “hydraulic jumps’”. 


I. Introduction 


From a meteorological viewpoint, consider- 
able interest is attached to the flow of a fluid 
with a vertical gradient of density. The at- 
mosphere itself is an example of such a fluid, 
and the existence of stable, unstable, and neu- 
tral internal gravity waves in a stratified me- 
dium of this kind must be of great meteoro- 
logical importance. A striking example is the 
“Bishop Wave”, a phenomenon now under 
investigation (COLSON, 1952). 

TEPPER (1950) and FREEMAN (1948) have 
developed models of waves on an interface 
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of two fluids of different density which appear 
to have application to atmospheric phenomena. 
In their investigations primary importance is 
attached to the instability of such waves with 
respect to the ultimate development of in- 
ternal “hydraulic jumps”. If the atmospheric 
analogy is correct, these correspond to “squall 
lines” or ““pressure-jump lines”. A model of a 
steady-state jump of this kind in a three-fluid 
system has been developed in the Hydrodynam- 
ics Laboratory of this university. Figure I is a 
photograph of this phenomenon. Further ex- 
perimentation along these lines is now under 
way and the results will be published in a 
later issue of this journal. 
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If a hydraulic jump be regarded as a mani- 
festation of fluid instability, it is apparent that 
investigations that aim only to discover the 
existence and the nature of infinitesimal waves 
in a fluid will not in every case succeed in 
predicting stability and instability. For ex- 
ample, in the case of the flow of water in a 
channel, perturbation theory of the above 
type leads to the conclusion that no unstable 
waves exist. In reality, however, violent jumps 
may occur. The cause of the discrepancy is 
well known and is due to the tendency for 
finite amplitude, long waves to steepen indefi- 
nitely as they move along, and to develop 
eventually into waves of the shock type. In 
the case of a stratified fluid, a similar oversight 
results from a classical instability investigation. 
In the absence of shear, all small perturbations 
are neutral and the regime is considered on 
this basis to be stable. For this reason, while 
drawing freely on perturbation theory, an 
attempt will be made in this paper to give as 
much consideration as is possible to phate wave 
motions. 

Lyra (1943) and QUENEY (1947) have studied 
stratified flows from the viewpoint of the 
stationary waves that mountain barriers might 
set up in a uniform current of stable Gad 
passing over them. Both of these authors 


re) 
considered the fluid as infinite in extent ver- 
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Internal hydraulic jumps. One jump is located on the second interface just downstream of the obstacle. 
larger one is on the first interface somewhat further downstream. The three liquids differ in specific 
and the internal Froude number, using the overall density difference, is 0.19. 


tically. This is equivalent to the assumption 
that the wave lengths of the perturbations are 
small compared with the depth of the atmos- 
phere. Since the waves set up by mountain 
barriers will have lengths of the order of the 
horizontal dimensions of the barrier, while 
the atmosphere may be considered to be 
about ro km thick, this is not a very good 
assumption. On the other hand, Queney was 
able to incorporate the effect of the rotation 
of the earth in his study of a stratified fluid. 
This presents a considerable advantage for the 
application of the results to meteorological 
processes. 

A paper by GÖRTLER (1943) considered this 

roblem from a somewhat different viewpoint. 
In his work the linearized equations of motion 
were examined and the flows classified into 
elliptic, parabolic, and hyperbolic regimes, 
depending on the form of the resulting system 
of differential equations. Some interesting -illus- 
trations are given in his paper of an experi- 
ment using a “stable salt solution in which the 

eal characteristic lines of the hyperbolic case 
are generated by an oscillating source of 
disturbances. 

More recently, RossBy (1951) and CRAYA 
(1951) have discussed the effect of stratification 
on the vertical concentration of momentum 
in a fluid current as a result of momentum 
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losses at the underlying boundary. Rossby de- 
rived a criterion for a critical, internal Froude 
number above which the main effect is an 
increase in the depth of the current, but below 
which there is a strong tendency for the mo- 
mentum to become concentrated in the upper 
levels. 

The purpose of this paper is to present 
some results, primarily of a theoretical nature, 
pertaining to this general problem. In section 2 
the exact, steady-state equations of motion 
of a current with an arbitrary initial vertical 
density and velocity gradients are integrated 
once. The resulting second-order vorticity 
equation is, of course, non-linear in general, 
but could probably be integrated numerically 
for any case of sufficient interest. The motion 
represented by the equation is examined quite 
generally in several ways. When thrown into a 
somewhat different form, a condition may be 
derived that the motion be uniquely deter- 
mined by the form of the boundary. In two 
special cases uniqueness corresponds to insta- 
bility, and a further investigation indicates 
that this may well be true in general. Sections 
3 and 4 are devoted to a detailed study of a 
model with a density gradient and uniform 
basic velocity. Section 5 deals with a certain 
type of shearing flow. 


2. General considerations 


Figure 2 is a drawing of the theoretical 
model considered in this paper. The motion 
of the fluid is considered to be the same in all 
planes parallel to that of the printed page. 
The y-axis is directed upward and the x-axis 
is along the general direction of the current. 
The fluid, while not homogeneous, is incom- 
pressible. At great distances upstream the 
velocity is parallel to the plane bottom and 
is an arbitrary function U(y). The density far 
upstream is given by o (y). Frictional forces are 
considered to be negligible except insofar as 
they may be required to make the problem 
determinate (see section 4). 

It is to be noted that in the model the Corio- 
lis forces vanish identically. This is a conse- 
quence of the assumption that the motion is 
strictly two-dimensional. If this assumption is 
dropped, one of the equations of motion 
contains a Coriolis term which bears a ratio 
to the inertial terms, of the order of w’/Rou'. 
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Fig. 2. Stratified fluid model. The velocity U and the 
density, 0, are arbitrary functions of the vertical coor- 
dinate alone, at a sufficiently great distance upstream. 
The upper and lower boundaries are rigid surfaces, 
although not necessarily plane surfaces. 


R, is the Rossby number, U/fL (L is a measure 
of the scale of the motion), u’ is a perturbed 
velocity component in the direction of the 
basic motion, w’ is a perturbed velocity nor- 
mal to the plane of fig. 2. If w’/Row’ is small 
compared to 1, the neglect of the rotation 
terms is justifiable. As an example, assuming 
that w’ is one-tenth of w’, U is 103 cm. sec.-! 
and f is 1074 sec.=!, the scale of the motion, L, 
should be of the order of, or less than, 100 km. 
In general, however, the validity of applying 
the results of this paper to the atmosphere 
must be determined by an individual investiga- 
tion along the above lines. 

Assuming a steady state, the equations of 
motion are 
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the two-dimensional motion. Eliminating the 
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Since the motion is incompressible, there is a 
stream function w (x, y) such that 
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The density, moreover, is conserved so that 


o —o(y) The equation (3) then can be 
written 
de ro ( 9q2/2 24/2 
dt o dy ox dy | 
do d 
o ay at 


Again, since o and the derivative of o with 
respect to y are not functions of time, we may 


write this as 
zc dpeyg* | = 
o dy Ê | o)] Fe rap) 


s(t: 
dt =) hl 


As a consequence of the steady-state assump- 
tion, this may be integrated to give 


[+0] = HW). ©) 


1 do 


Me ay - 


o dy 


where H(y) is a function of y to be deter- 
mined by the conditions far upstream. Thus, if 
yo (y) is the height of the streamline y = const. 
and Ë(y) is the vorticity upstream, the final 
result is 


do (Vy)? 
vey + - Wel = f(y) + 
gdp 2 
: (9) 
+22 [+0 |: 
el, ot 


It should be noted that £, o, U, and y, are all 
known functions of y, and hence of y. If 
proper consideration be given to the boundary 
conditions and other conditions which may be 
required to insure a unique solution, it seems 
quite feasible that (9) could be integrated nu- 
merically for any given basic distribution of 
density and velocity. It is to be noted that the 
assumption that the perturbed motion vanish 
at a sufficiently great distance upstream may be 
violated by some solutions satisfying (9). Such 
solutions will nevertheless satisfy the primi- 
tive equations of motion and continuity and 
are, therefore, dynamically possible. 
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The above equation assumes a familiar 
form if the fluid is homogeneous. Since 0 
is then a constant, (9) becomes 

Vip =¢(y). (ro) 
This expresses the conservation of vorticity 
in the absence of solenoidal forces. It has 
been shown by Long (1952) that an analogous 
equation, with the right-hand side a linear 
function of y, adequately describes the basic 
features of the perturbed motion of a spherical 
shell of rotating liquid. The waves in the rota- 
ting shell are “Rossby waves” and owe their 
existence to the special character of the basic 
absolute vorticity field. Similar wave motions 
may exist that satisfy (10) but they will be of a 
more general character, depending on the 
form of £ (y). If the basic vorticity is zero or a 
constant, however, the perturbed motion will 
be irrotational and no wave motion can occur. 

Equation (9) may be thrown into a different 
and useful form by substituting y, for the 
dependent variable y. If U does not vanish 
anywhere, the relation, 


d 1, (11) 
dp U dy” 
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With the further substitution, 6 = y, — y, 
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In this form we see one obvious model for 
which the equation is linear. This is one in 
which U?o = const. and the density is linear 
in yo. Equation (13) then reduces to 


V28 +o?8=0, t= — EM = 
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This is the only case I have been able to dis- 
cover for which the differential equation 
governing the motion of a stratified fluid is 
exactly linear. 

Equation (12) may be used to derive a 
sufficient condition that a given solution for 
the motion in a channel be a unique steady 
solution. To show this, we replace our 
present model with a circular channel of very 
great radius in order to remove the necessity 
for defining the stream function at infinite 
distances along the x-axis. In the altered model 
the function y,(x, y) has cyclic continuity. 
With this understanding we may, however, 
neglect any centrifugal forces which arise, 
since the radius of curvature is made arbi- 
trarily large. 


Multiplying (12) by Ug? we obtain 


Assume now that there are two solutions ygı, 
Yoo of (12) both of which are continuous and 
have continuous first and second derivatives 
and which take on the same values at the 
(perhaps irregular) bottom and top of the 
channel. Corresponding to these functions are 
yo) and f(yos), defined by the second equa- 
tion in (15). They are obviously equal on the 
boundaries so that fy =f(yo2) —f(yo1) vanishes 
at the top and bottom of the channel. Substi- 
tuting in (15) and subtracting, we obtain 
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where the partial derivative of G(f, y) is 
with y held constant. According to the mean 
value theorem (COURANT, 1947), this derivative 
is evaluated at some value of f between f(yo;) 
and (yo). Multiplying (16) by fy and inte- 
grating over the whole region, 

Sv vor 5 6°] dedy = 0 
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In view of the cyclic continuity and the dis- 
appearance of f on the boundaries, this re- 
duces to 


ff (io) oe dxdy=o. (18) 


A sufficient condition that fp = 0 or f(yoa) = 


(Yo) is that 
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If f, == 0, however, using (15), 
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Since we have assumed that U does not van- 
ish anywhere, the integrand is of one sign 
everywhere. Hence, yo2 = Yo; and the solution 
is unique. The inequality, (19) is, therefore, 
the sufficient condition that the motion in 
the channel be uniquely determined by the 
form of the boundaries. Expanding (19), the 
condition may be written 
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This can be simplified and made somewhat 
more intelligible if we assume that the dis- 
turbed motion is small. The last term in (2r) 
may then be neglected and the condition 


becomes 
1 dU do 1.42% ur fd 2 
Uo dy, dy U dys? 40° \dyo | 
t do g do 
F 2 = 22 
20 dy? me dYo a ee 


The usefulness of this criterion is limited by 
the fact that it is not necessary and is, in fact, 
a weak condition in many cases. If, for example, 
we examine the model in which U is uniform 
and the density is © = 0, exp (— By) the in- 
equality becomes gB/U? < 2/4. In practical 
applications 6? can be neglected and the uni- 
queness is assured only if B is zero or neg- 
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ative, i.e., if the density is constant or in- 
creases upward. It will be shown in the next 
section that a much stronger condition can 
be obtained in this case if the total depth of 
the channel is limited. Physically, the present 
condition is of interest, however. If the 
upper and lower boundaries of the channel 
are plane, one solution is the basic solution, 
U=const. If the density increases upward, 
this is the only possible steady solution cor- 
responding to the above inequality. Thus, any 
disturbance superimposed on the basic motion 
must progress, die out, or grow in amplitude. 
Since the stratification is definitely unstable, 
we know in advance that the last possibility 
is the one which will be realized. 

A non-trivial condition is obtained if we 
assume the same density distribution as above 
but take as the velocity, U= U, exp ay. 
A given solution is unique then if 


Since ß/a is small in cases of interest, and Ua 
is the velocity shear, this can be written 


Ra (=) | /26. (24) 


This expression is the Richardson number. 
According to RICHARDSON (1920) if this 
number exceeds one, the flow is unstable in 
the sense that turbulence will increase. We 
therefore have another case where uniqueness 
is associated with instability. In the turbulent 
flow examined by Richardson, the shear is 
that of the mean motion, of course. 

It is possible to indicate, in general, that the 
criterion for uniqueness may also be a sufficient 
condition that a given steady motion in the 
channel is “secularly” unstable. This concept 
of stability was introduced by Thomson and 
Tait (see Lams, 1932). The distinction between 
this and “ordinary” stability is that a given 
system may prove to be stable in the ordinary 
sense when considered frictionless but un- 
stable in the secular or permanent sense when 
dissipative forces, however slight, are permit- 
ted. A simple, but beautiful, example of this 
is given by Lams (1908). He investigated the 
stability of a small ball moving freely inside a 
hollow sphere rotating about a vertical axis 
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in a uniform gravity field. If there is no fric- 
tion whatever between the ball and the bowl, 
it is obvious that the only stable position is at 
the lower pole of the sphere, since the rota- 
tion has no significance in the problem. If 
there is any friction between the two, how- 
ever, this position becomes unstable when a 
certain angular velocity is exceeded, and the 
ball will move to another latitude on the 
sphere. This second position is, then, secularly 
stable. The “practical stability” of the system 
in this example is assured by the condition 
that them neue potential, Tr Les 
minimum, where V is the potential energy 
and T, is the kinetic energy when the ball 
is at relative rest. Since 


V—Ty = —MgR cos 0 — ~ Mo? R? sin? 0, 
(25) 


we have an extremum when sin 9=o and 
when cos @ = g/m?R. The condition that the 


first of these be a minimum is that © < (g/R)?- 
If the angular velocity exceeds this value, 
Lams showed that a ball starting initially at 
O=0o will, under the action of friction with 
the bowl, spiral outward from the pole until 
it reaches its secularly stable position of equili- 
brium. This action of the ball meets the ordi- 
nary specifications of instability in that the 
system departs permanently from its equili- 
brium when disturbed slightly. 

Helmholtz (Lams, 1932) applied this con- 
cept to a fluid system consisting of a homo- 
geneous liquid in a channel with a free surface. 
He found that the condition for steady motion 
was that V—T be stationary, where V is the 
integrated potential energy and T the kinetic 
energy. Furthermore, if the flow is sub-critical, 
Le. vt Ue (gH)? this quantity is a minimum. 
If the flow is supercritical, V—T is a maxi- 
mum. A prior analysis indicates that if a small 
amount of friction is introduced, and it is 
assumed that it acts only to reduce the eleva- 
tions of the free surface, then the maximum of 
V—T means secular instability. This is some- 
what surprising at first glance. Since the 
friction does not affect the flow directly, it 
is obvious that the choice of coordinate sys- 
tems from which U is measured is immaterial. 
If we choose one that moves with the fluid, 
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then the above result seems to show that a 
fluid basically at rest is unstable! The ariswer 
to this paradox seems to lie in the nature of 
the perturbation introduced in testing the 
maximum or minimum nature of V—T. If 
we introduced a free perturbation moving 
faster than ( gH)! on a resting fluid, it would, 
of necessity, have a finite amplitude, since 
the maximum speed of an infinitesimal wave 


is (gH) >, But it is well known that such a 
finite wave will develop into a surge (moving 
hydraulic jump) so that the resting fluid will 
be unstable in the sense considered in this 
paper. The perturbations considered in varying 
V—T, however, are those which are steady 
with respect to the coordinate system from 
which U is measured (although they need 
not be dynamically possible steady motions). 
If the value of V—T is greater than in any 
disturbed state of this type, then the conclu- 
sion of Helmholtz is that the flow is unstable. 
In channel flow this instability is manifested 
in the form of hydraulic jumps. 

Such a test of stability is a logical one from 
the viewpoint of the flows considered in this 
paper. We are interested primarily in the 
behavior of a stratified fluid passing over 
obstacles at the bottom. Since the upstream 
conditions are steady, the resulting theoretical 
perturbations will be steady. Whether or not 
they are stable will be tested by an approach 
similar to that used by Helmholtz. 

The application of this principle to the 
model under consideration will be made in 
the following indirect manner: we first show 
that the dynamic differential equation (9) 
that must be satisfied by steady-state flow in 
the model can be derived by a variational 
method. To accomplish this, we seek the 
condition that the following integral be sta- 
tionary : 


1 = ff [ext ee. 


2 


+ f egdy,—e ve dxdy. (26) 


In the variation the altered streamfunction 
w(x, y) +e(x, y) is assumed to be steady, and & 
is to vanish at the top and bottom of the 
channel. Thus, both y and yp + ¢ satisfy the 
kinematic conditions at the boundaries. We 


again replace the straight channel with a 
circular one of arbitrarily large radius so that 
the assumed motions have cyclic continuity. 
There are a few points of delicacy in this 
variation problem that require a development 
from first principles. Thus, we write out the 
expression for the integral, I(y + e). This is 
accomplished by expanding y and functions of 
y in the integrand in a Taylor’s series in powers 
of &, neglecting terms of the order e?. For 
example, 


dU? aU? & 

Lia ee 
CONTE OSG ey 
| d dy, & 
Yo(y + €) = yolp) + £ + ae 2 
u Oa! dU & 
> Yo(¥) Di j U2 dy 2° 


(27) 


To the order ¢? we then obtain 


Ip + e) = I(y) 4 Frlervre. 


vy)? U: du 
{| 2 t gly | eU ldxdy 


2 2 dy | 
A Ve de | 

es: {3 eo ae 
Siegetad? 1 do 
2 = [U? aa eae er 
Ake dU\? 0 „Fu dU do 
2 \dy de dp dy 

nd \ dxd 
2) Urdu ln y 


The condition that the integral be stationary 
is that the coefficient of ¢ vanish. This yields 
equation (9). Hence the condition that the inte- 
gral I(y) be stationary is that the function y 
represent the stream function of a dynamically 
possible steady motion. This is analogous to 
Helmholtz’ result for channel flow. In his 
case the varied integral was the integrated 
difference between potential and kinetic ener- 
gies, V—T. The integral I(p) in (26). con- 
tains —T in the last term in the integrand. 
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The remaining terms may be regarded as a 
generalized form of potential energy. This 
may be inferred from the Bernoulli equation 
for > stratified model, which is casily verified 
to be, 


a ea) | 
+ S ogdy, = const. 


For comparison, in steady, potential flow the 
equation is 


Pes 0 + V’ const; 


(30) 


V’ being the potential energy per unit vol- 
ume. It is not apparent what physical inter- 
pretation is attached to the “potential energy” 
terms in the Bernoulli equation (29). Some 
insight is obtained from the above example of 
the ball moving in a rotating sphere. In that 
case the expression V—T, was varied, Ty 
being the kinetic energy of the ball by virtue 
of the basic rotation. From the viewpoint of 
an observer in the sphere, however, —T, is 
nothing but the potential energy of the cen- 
trifugal “force” set up by the rotation. V—T), 
in this case, can very appropriately be regarded 
as a generalized potential energy. 

The remaining question to be considered 
in the variation problem is the condition that 
the integral I be a maximum or minimum. 
Evidently it is a maximum if, in (28), the 
coefficient of ¢ vanishes and the coefficient of 
e? is positive. A stronger condition is ob- 
tained, however, as follows: Since its integral 
over the region vanishes, we may add to 
the integrand of the second integral in (28) a 
term, 


e? do do & 

V = = Vy= — — Ves Vwt 
4 dy ‘? y 2 H 
e do do 
N? 31 
em 60 


Combining terms, (28) becomes: 


I(y + €) = I(y) — 


I A Ning ei 
<2 ff e(ve+ a vr) xdy 
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dor (do EU? dU\2 
ne +s-2)]+e(5,) à 
AU 
dy? 


dU do g do | | 
a. . dxdy. 


U dy | 
(32) 


The integral I will again be a maximum if the 
coefficient of &? is positive. When the terms 
in the first square brackets in this coefficient 
are replaced by using equation (9) (or the 
condition that the coefficient of & vanish), we 
obtain precisely the uniqueness condition (21). 

I have not been able to demonstrate that the 
uniqueness condition implies instability in the 
secular sense, although it would seem from 
the above discussion to be very likely. It is to 
be hoped that further investigation will 
establish rigorously this suggested relationship. 


+ 0U. 


3. General observations on the flow of 
liquid with uniform velocity distribution 


One of the simplest cases included in (9) 
is one with a uniform velocity, U and a den- 
sity distribution given by @ = 0, exp (— fy). 
We obtain 


vey + Fam): = : [= F £(Yo D 
(33) 


If we introduce the perturbed stream function 
y = Uy + y, this becomes 


vi + À (vy —B 


Op ae ed 
dy U2 ce “ 


(34) 


The middle terms in this expression are ordi- 
narily exceedingly small. Thus, if the length 
scale of the phenomena is of the order of H, 
the ratio of the third to the last term is F?, 


where F is the Froude number U/(gH)?. In 
atmospheric motions, if H is taken as the 
height of the tropopause, F? is of the order: 
of 10-? to 10-8, Furthermore, we may carry 
this term in the following development, since 
it is linear, but its effect on the motion 1s 
again found to be minor. The ratio of the 


so 


second term to the last term is of the order 
of F26/H, where 6 is the deviation of the 
stream line from its equilibrium height yp. 
With a high degree of approximation, there- 
fore, we may substitute for (34) the wave 
equation 


V2y + ay’ = 0, a = an (35) 


Some important information on the type of 
motion represented by this equation can be 
obtained by assuming that the bottom of 
the channel is corrugated in the form of 
sinusoidal “bumps” of wave length 2x7/k and 
amplitude y. Writing 


y = > À,cos skx sin| (a8 um s2k2)2 (y — H) | 


co (36) 


we have a general solution which vanishes at 
the top of the channel. The bottom boundary 
condition becomes 


Uy cos kx = % À, cos skx 


s=0 


sin EG — s2k?)2(y cos kx — H) | (37) 


If y is small, we may neglect higher powers. 
Assuming this, we obtain 

Oy 
sin (a? — k?)? H 


sin [ (=)! — |, 


cos kx 


yp == 


(38) 


for a boundary form y cos kx. In this solution 
the oscillations have the same wave length as 
the bottom configuration, but the amplitude 
varies with height, and the waves will be 
inverted relative to the bottom waves after 
some vertical distance, provided («*—k?) H? 
> n?. Defining an internal Froude number, 
F; = U/(gßH?)®, this inequality becomes F;? 
<(x?+k2H?)-1 If F; exceeds 2-1, however, 
‘an inversion cannot occur, regardless of the 
wave length of the corrugations. This value, 
m1, of the internal Froude number would 
seem to be very significant. The analogous 
problem of the flow of a homogeneous fluid 
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with a free surface over a corrugated bottom 
reveals that the criterion for the inversion of 
the free surface deformations over the bottom 
deformations is the value of the Froude num- 
ber, U/ (gH)? = 1, which separates subcritical 
from supercritical flow (ROUSE, 1938). 

In the case of open channel flow of a homo- 
geneous fluid, a definition of critical flow 
more meaningful physically is the speed U, 
which coincides with the velocity of propaga- 
tion of long gravity waves. If the stream 
velocity is less than this, long disturbances 
created by irregularities at the bottom will 
travel upstream until destroyed by frictional 
losses. If the stream velocity is greater than the 
critical speed, disturbances will be swept down- 
stream, if of sufficiently small amplitude. 
Since the crests of finite long waves move 
faster than infinitesimal waves, the former 
will be swept downstream more slowly, ab- 
sorbing the lower waves that come by. If, 
through this method of growth, a finite wave 
becomes high enough, it may reach a state in 
which it can maintain its position against the 
current. Since an elevation in such a wave 
moves more rapidly than a depression, there 
will eventually be a “breaking” and a “hy- 
draulic jump” will then exist. 

The importance of long waves in this 
discussion is that they are the waves which 
steepen as they progress and, therefore, even- 
tually break, with the accompanying turbu- 
lence which characterizes jumps and surges. 
Short waves (deep-water waves), on the other 
hand, can move without change of shape when 
they have a finite amplitude. 

When applied to stratified flows, these 
concepts become more complicated. The 
approximate speeds of infinitesimal waves in 
the model considered in this section are given 
by c? =g8 H?/(n?x?+k?H?). If we then define 
as a critical speed, the speed of long waves, we 
find an infinity of such speeds given by gf H?/ 
n?r?, where n is an integer greater than zero. 
The derivation of the wave formula shows 
that n is related to the number of nodal sur- 
faces within the channel for the particular 
wave. The number of these surfaces is, in fact, 
n—1. When n = 1 the waves extend from 
bottom to top, and the corresponding critical 
speed is given by F;= 2-1. A flow with an in- 
ternal Froude number slightly less than x 
will be supercritical with respect to waves 
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having one or more nodal surfaces between 
the bottom and top of the channel, and the 
analysis of water flow in a channel seems to 
indicate that internal jumps could result from 
the growth and breaking of these waves. 
The results of a very preliminary experiment 
using an obstacle in a three-fluid system support 
this tentative conclusion. A hydraulic jump at 
the second interface occurred just above the 
barrier for a Froude number of 0.14, which is 


well below the Froude number, (1 /6), corre- 
sponding to the fastest long waves which can 
exist in such a medium (see the discussion of a 
multi-fluid system at the end of this section). 
The experimental jump had a rather small 
amplitude and it seems likely that this would 
be true of all jumps that develop from the 
higher modes of oscillation (n > 1). The above 
suggestions have not been supported by rigo- 
rous theory or by careful experimentation, 
however, and should be regarded with reserve. 
It is to be hoped that the future experiments 
will shed light on these important questions. 

It is of some interest to derive the upper 
critical Froude number by the use of another 
analogy to water flow in a channel. Thus, if 
an obstacle is inserted at the bottom of the 
channel, the resulting wave form at the free 
surface is indeterminate for subcritical flow, 
but entirely determined if the flow is super- 
critical. (The indeterminacy can be removed 
by appeal to viscosity). In the case of stratified 
flow, we may inquire into the conditions 
under which a solution is determined unique- 
ly by specifying the bottom topography, 
y=ys (i.e. the values of y’ on the boundary), 
and the surface y = H as a streamline. As in the 
previous section, we assume that the fluid 
circulates in a circular channel of very great 
radius so that there is cyclic continuity. If 
we use the same arguments as those in deriv- 
ing (21), then, assuming two solutions y’ 
and y, of (35) whose difference is po, we 
find that 


Lf oF? vo + @ yo? dxdy =o. (39) 


The integral is taken over the entire region of 
flow. Integrating by parts, 


— ff (yo)? — dy? dxdy=o. (40) 


Now a quantity, V - y2Vn (x, y), when inte- 
grated over the region, will vanish if 7 is 
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sufficiently well-behaved in the region. Add- 
ing the integral of this expression to (40) we 
obtain! 


IN (9 vo —voV 1)? — Yo? [9 0)? as 
+ Vin +a] dxdy =o. 


(41) 
The integrand of this expression will be every- 
where non-negative (i. wy must vanish 
identically) if 

(Mn a VW + a" =O, (42) 
where «,? is slightly greater than or equal to 
a, Taking 7 as a function of y alone, the solu- 
tion is 

n = In [A sin (my + C)], (43) 
A and C being integration constants. Taking C 
as — 7x, then, provided that «, H <a, the func- 
tion 7 will be regular in the region and we 
must have a unique solution. This is equiv- 
alent to F; >a, the condition for super- 
critical flow derived above. 

The work of Tepper and Freeman, men- 
tioned in the first section of this paper, is 
concerned with the behavior of a system 
composed of a discrete number of layers of 
fluid of different density, rather than with 
continuous density distributions. It would 
seem that such a system is often a good approxi- 
mation to the atmosphere and, in any case, 
is mathematically more convenient than any 
continuous distribution that could be fitted 
to the actual conditions. Furthermore, it 
appears that the most practical experimental 
model that could be set up to examine effects 
of stratification would be one with a given 
number of immiscible fluids of varying den- 
sity layered in a channel. It is of interest, 
therefore, to investigate such discontinuous 
systems theoretically to obtain criteria for 
critical regimes of motion. 

The following investigation was suggested 
by an approach due to R.R. Webb and 
reported by GREENHILL (1887). The critical 
speed will be taken as the speed of the fastest 
long waves which can exist in the stratified 
medium. In Greenhill’s paper it is shown that 


1 This method of establishing the uniqueness of a solu- 
tion is given in RIEMANN-WEBER (1925). 
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if a channel of n+r equal layers of fluid is 
bounded above and below by rigid planes, 
the speeds c of long waves are given by the 
equation | A| = 0, where | A| is the determinant: 


Ao) 0, 0 
B, Ab; 00 
Oo By Ab B, - = ¥ (44) 
Ba A Ds 
O De An 
and 


In the last equations h is the (constant) depth 
of each layer, 0, is the density of thetoplayer, 
and A is the (constant) difference of density 
from one layer to the next. The problem of 
solving the determinant is greatly simplified by 
assuming that A is very small (which it will 
be in any physical case of interest). If this 
assumption is made, it can be seen that the 
wave speeds will be given by an equation, 


DE ae (46) 


where b is constant. Substituting back in the 
determinant the equation for b is 


(2b—1) —b o 0 
—b (2b—1) —b o 


O —b (2b—1)—b- -|=o. 


(47) 


This determinant has n rows and columns. An 
equation of this type was solved by RAYLEIGH 
(1894). His method shows that 


T2 (48) 


Of all these waves, the fastest are given by 
s =n. The critical speed becomes 
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ghA 5 177 
Pie ee 49 
d 4 Q1 2(n +1) (49) 
In order to compare this to the results of 
the calculation for a continuous density dis- 


“ tribution we note that, if the density gradient 


in (22) is small, BH, in the continuous case, is 
close to 4e/o,, where Ag is the overall 
density difference from top to bottom and 
0, is the density at the bottom. The critical 
speed is then given by 


U2 me: (s6) 
2HAojo man . 


To the same order, since h = H/n + 1, and 
0, =  —nA, (46) becomes 


GC I 


gHAola, ann +1) 


el a 
2(n+1) n> co 72" (51) 


We thus obtain the same critical Froude num- 
ber for this case as for the continuous density 
distribution. The critical numbers in the cases 
of two, three, and four layers are 0.50, 0.41, 
and 0.38 respectively. 


4. Flow over an isolated barrier 


In this section we will examine the flow of a | 
stratified fluid with uniform velocity distri- 
bution over a small isolated barrier in the N 
bottom of the cannel. Using equation (35), ify 
is infinitesimal and the bottom deformations 
are given by y, = y exp (ikx), the resulting 


stream function is 
; Uy 
Ve 1 
sin («2 — k2)?H 


[Rp An) (52) 


The resulting linearity permits a build-up of 
functions satisfying the kinematic boundary 
conditions of an arbitrary bottom profile. 
Thus, if the bottom topography is 


eikx sin - 


y= S vll) et dk, (53) 


the solution is 
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y= ES ployer de’, (59) 
so that 
bi [ee] [ee] 
De J ak Sy.) - 


le —F)#(y — H)] ned. (56) 
sin [(x® — k2)® H] 


A simple, but interesting, boundary form 
may be taken as a low, rectangular obstacle 
with a height, a and width 2b. Equation (56) 


then becomes 
co 
Ua eik (x + b) 
2 Hi k 
— 009 


sin [(x2 — k2)®(y — H)] d 


; Er 
sin [(«? — k?): H] 
eye? ra eik («—b) 
MD gt i; k 
sin | («2 — eee H){ Boe a 


(57) 


The evaluation of these integrals is possible 
by contour integration in the complex plane, 
£ =k + iv. For example, 


» ax +D sin [(oe2 aes £2)2 (y — H)| 


= dé, 
“ Ë sin [(x? — &)2 H] : 


(58) 


for any closed contour which includes the 
whole of the real axis, v= 0, includes the first 
integral in (57). The choice of the rest of the 
circuit depends on the sign of «+ b. If this is 
positive, it proves advisable to adopt a con- 
tour consisting of a semi-circle, |&|=R, in 
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Fig. 3. Circuit in complex plane for evaluation of 
equation (59). 


the upper half of the complex plane, plus the 
region of the real axis from — R to R (see 
fig. 3). Then, letting R >, we would ex- 
pect the integral along the real axis to converge 
to the first integral in (57). The success of the 
method depends on the behavior of (58) on 
the semi-circle. If x +b> 0, it may be shown 
that this part of the cyclic integral tends to 
zero as R — co. For x +b <o we use a contour 
along the real axis and around a semi-circle in 
the lower half of the complex plane, &. 

The above procedure is straightforward, but 
laborious, and we will indicate below the 
evaluation of y,’ in (57), for «+6 >0, and 
merely write down the results of the other 
integrations. Setting € = Re” on the semi- 
circle, (58) may then be written 


I=i/ exp [iR(x + b)e?] - 


sin [2 — Rte} (y — HD] 
sin |(a? — R? eri) 3 H] P 


R 
eik (x + b) 
ae ii k N 


—R 


_ sin [8 — BG H) gr pr, 
sin [(a2 — k®)? H] | 


(59) 


The integrand of J, is singular at a finite num- 
ber of points, 
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27? 


k= 0, k, = ee ie n=0,1,2,...N 
(60) 


where m, is the largest positive integer for 
which k, is real. We enclose these points ın 
small. semi-circles of radius r and use the 
resulting circuit, shown in fig. 3, to evaluate J. 

The value of the last integral in (59) may 
then be taken as its principal part, or the limit, 
as r > 0, of the integrals along the parts of 
the real axis excluded by the semi-circles. There- 
fore, 


_ lim PI, _ lim I— lim I, 


Mes IR es [RES CS) 


+ x f cine +0 sin (Cr 20 — 2 > 
5; sin [(@2 — &,2)? H] 


(61) 

In (61) the symbol, P, means “the principal 
part of”, the summation is over all semi- 
circles on the real axis shown in fig. 3, and 
en == Ry An Te». 

The cyclic integral I over the modified 
circuit has a value equal to the sum of the res- 
idues at its poles. These poles are infinite in 
number and are given by 


, (nr? 2\? 
Be Mm admet 


The result is 


D Mer 2m. 


(62) 


Sn: né 


n=m+1 


272 L 
“exp -& — 2) (x + ) | sin a (y—H). 
(63) 


As remarked above, for x + b> 0, lim I, =oas 
R increases indefinitely. The integrals over 
the small circles in (61) yield the following: 


A sin «(y — H) 
sin « H 


em ann? 
Hi PL Si 2 H? — n?n? 
2\ 3 
cos] (at TT ) (x + b) js (y —H) 
(64) 


Using (61), (63), and (64), we obtain 


Ua sin «(y — H) 
2 sin «H 


pr (x +b> 0) =- 


— = S Se \tiatet = Uanıa _ 
2 2 ( 1) a? H? — n?n? 
an 3 b H 

cos | | — ie (x + b) | sin — (y —H) + 

ded = TPS Uanx e 

D 2 1) nr? — a? H? 

2 \3 

exp [= m ) +0] sin — (y—H) 
(65) 


Applying the above method to the second 
integral in (57): 


Ua sin «(y — H) 


al ER ee, | 
Pr © We ) 2 sin «H 
pee 2 Uann 
el Re 2 HET L LRO 
Ne 2 ( 1) a? AY — n°7? 
. cos] Lu (+0) ] sin — ( — H) — 
H? À H VV 
nf S He _ 2 Uann 
2 ZS 2) nagı = a2 FI 
n? 7? 4 NIT 
BE Oe 2... ee 7.0 
exp] (3 2) (x | sin av H) 
(66) 
Ua sin «(y — H) 
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Pe ( ) à sin «H 
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2 Uanx 


oo 
I 
2 = 2 Eu ie nn? — a2 H? 


es 
(67) 


We may combine (65), (66), and (67) to 
obtain the solution for the two regions: x > b, 
and—b< x<b.Thesolution for x <—b follows 
at once from the observation that (57) re- 
quires that y be symmetric about the line, 
x = 0. 

The assumption that the obstacle has a very 
small height permits the addition to the parti- 
cular solution above, of any free oscillations 
that satisfy the conditions that y’ vanish at 
y=oand y=H. Itis obvious that ee do not 
exist if aH <a. Therefore, if F;>27} the 
solution is unique. 

The indeterminacy that arises when F;< 27! 
can be removed by considering the analogous 
problem of subcritical flow over barriers in a 
water channel. An investigation by RAYLEIGH 
(1883) revealed that the indeterminacy is elim- 
inated by introducing viscosity. Then, as 
the frictional coefficient tends to zero, the 
resulting asymptotic solution is fully deter- 
mined. This ‘ ‘practical solution” differs from 
the perfect fluid motion only in that the wp- 
stream waves are absent. Adopting this to our 
case, we annul the upstream oscillations by 
adding to the whole channel a series of free 
oscillations: 


ny 


= ( 2 Uanz 


yr : > HA 7 
I - sin | x? — . 
FT a? H? — nn? Jal? 


: A HE NIE 
om (a2 —— | x sin — 
teh H 


The final solution is 


yi (x>b) = 


s > (—1) 2 Uana 
n=m+1 Bere x? H? 
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een sin.a(y | 
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; 127? q 
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a. 
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(71) 


When F; < 2! the obstacle produces down- 
stream waves. The amplitudes are propor- 
tional to n/(F;-? — n?r?), which is a maximum 
for n=n,. When n,=1, the total motion is 


Fig. 4. Streamlines at some distance downstream of 
an obstacle whose area in the plane of the motion, 
Q = H?/10 x. The internal Froude number, F? = 1/5 n°. 
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rather simply described. One part is a local 
disturbance, dying out rapidly on both sides 
of the obstacle; the other is a single wave of 
léneth = 27 NF 72)”: downstream, 
When #, > 1, two or more waves will exist 
with respective strengths indicated in (69). 
Figure 4 is an example of the flow when F? = 
1/5 2, (n, = 2). It represents a section of the 
channel sufficiently far downstream that only 
the oscillatory terms in (69) are sensible. A 
short calculation shows that the wave terms 
are already of the order of ten times the re- 
maining terms in (69) at a distance H/4 down- 
stream. 


5. The effect of shear on the motion of 


a stable fluid 


Several serious objections may be raised 
to any application of the results of the pre- 
vious section to the atmosphere. Normally 
the flow from the surface to the stratosphere 
has a very considerable shear of velocity with 
height. In the westerlies, for example, the How 
is roughly in the same direction at all levels 
but several times faster aloft than near the 
surface. 

Many studies have been made of the stability 
of shearing flow with density gradients, no- 
tably by RICHARDSON (1920), TAYLOR (1931), 
and GOLDSTEIN (1931). The stability of the 
motion is assured, according to Richardson, if 
the number 


This result is dependent on the assumption 
of the conservation of momentum of the 
perturbed motion. 

Tylor’s analysis was greatly complicated by 
mathematical difficulties. He could only state 
in fact that, for a fluid extending above a 
rigid plane to infinity, progressive waves can 
exist if R; < 4, and that no waves of any kind 
can exist if Rh 4. 

In this section an analysis similar to the 
one in section 3 will be applied to a current 
moving between two rigid walls with a den- 
sity distribution, o = 0, exp(—f yy), and with 
U= U, exp (ay). This velocity distribution is 
somewhat different than T aylor’s (linear 
shear), but the advantage is to avoid some of 


the mathematical difficulties of the latter dis- 
tribution. 

If we assume that the perturbed motion is 
small, we obtain from (13) 


20 gh 
V20 + (24 p) ay | U2 024 ©. 


(73) 


Assuming 6 = y(y)e**, the equation for y is 


Fr gb 2 | — 
BF Go k ) 0. 
(74) 


x +y (24 


The solution may be verified to be 


5 = exp |— (« == 4 y | 
4 gb 2 a: | 1 2 ke — B ; 
zu (4) ea, m 2 +(1 


(75) 


where Z, is a cylinder function, 
Zn = Za be ir Br INN: (76) 


If we consider that the fluid is bounded by a 
wall at y=o and y=H, then the possibility of 
free, steady perturbations depends on the 
existence of two zeros of the following func- 
tion: 


Am(Ri, aH) = Ne (ere Ri=) Jn(RF?) = 


=> Nn(Ri-*) ee (ene R;*). (77) 
Its zeros are given by JAHNKE and EMDE (1945). 
The tables indicate that the first zeros increase 
in size with m. The index m, however, has a 
minimum of 1 if we assume that B < <a, a 
reasonable assumption in most cases of interest. 
The largest values of R; permitting this type 
of motion are therefore those corresponding 
to the zeros of (74) with m= 1. It appears 
that the results of this investigation are much 


fe) 


clearer if, instead of seeking critical Richardson 


numbers, we consider the values of the function 


2 
a 


(e# — 1) ioe D 


K= LES 
et Ri? F; UÜ Un 


(78) 


when the function A, is zero. F; is the mean 
Froude number, U/(g8H2)3, U is the mean 


Tellus V (1953), ı 


FLOW OF STRATIFIED FLUIDS SH 


velocity of the How, U, is the velocity at the 
top of the channel, and U, is the velocity at 
the bottom. The values of K corresponding to 
the zeros of A, vary very little when aH 
changes from zero to ®. For the first zeros K 
increases from a at aH =o, to only 3.8317 


at aH= oo. F; changes even more slowly so 
that, if the velocity at the top is less than 


about 4 Up, the critical value of F; is given by 


F; = (79) 


Us | H 
Lam: 
He 
fe) 
D 
sr 


Such a range will probably include all cases of 
physical interest. When the shear is zero the 


critical value of F; is at, as we would expect 
from the results of the previous section. We 
may formulate the above results as follows: 
The criterion for supercritical flow in a channel is 
F; = U/(gBH?)* > 1/3, approximately, where U 
is the mean velocity of the flow. 

This conclusion is based, of course, on the 
assumption of a particular velocity distribution 
in the channel, namely, U = U, e™. It would 
have been preferable, no doubt, to have 
carried through the analysis with a linear 
shear, but this presents considerable mathe- 
matical difficulties. It does not seem likely, 
however, that the above result would be 
greatly altered by this change in the model. 
Certainly, if the shear is small, such a general- 
ization must be valid, since the shear will then 
be practically linear. 

The above results seem to indicate that 
critical states of shearing flow, and quite 
possibly the instability of a given flow, de- 
pend not alone on the Richardson number 
but also on the depth of the flow if the depth is 
finite. In the case of the shear considered in 
this section, a channel of infinite depth has a 
critical Richardson number of approximately 
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one (see section 2). The critical state in a finite 
channel, however, is best described by a mean 
Froude number. 


6. Meteorological implications 


The models investigated in this paper are far 
simpler than anything likely to be encoun- 
tered in the atmosphere. The imposition of a 
rigid upper surface is a very important feature 
of the problem, yet in the atmosphere the 
closest approach to such a surface is the tro- 
popause, which is more like a fluid interface. 
In addition, such effects as those due to com- 
pressibility and the earth’s rotation are neglect- 
ed. The results of the investigation, there- 
fore, are not likely to have much quantitative 
application to the atmosphere. Regarded 
simply as an investigation of internal gravity 
waves, it is probable that certain of the phe- 
nomena examined, in particular internal hy- 
draulic jumps, occur in the atmosphere and are 
qualitatively explainable by the approach used 
in this paper. If so, more reasonable models 
can be examined based on equation (9) which 
is valid quite generally. In addition, other ex- 
tensions of this kind are obvious. For example, 
the procedure used in deriving critical flows 
in a multiple-layer system could be used to 
find such criteria in a three-fluid system ap- 
proximating the atmosphere. The layers need 
not be of equal depth and the first two layers 
would represent the troposphere with a sub- 
sidence inversion. The upper layer could be 
used for the stratosphere. 
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Influence of Suspended and Dissolved Matter on the 
Transparency of Sea Water 


By N. G. JERLOV, Fishery Board of Sweden, Gothenburg 


(Manuscript received November 13, 1952) 


Abstract 


The transparency at different levels in the sea as recorded by means of a transparency-meter is 
compared with the distribution of suspended particles studied with the aid of a Tyndall-meter. 
À simple procedure is used to distinguish between different components of light extinction 
and to derive values of the content of the dissolved yellow substance. The distribution of yellow 
substance is illustrated for the transition area between the Baltic and the Skagerrak and for the 
outflow of the Göta Älv. It is emphasized that the content of yellow substance may often be 
considered as a characteristic property of a water mass. 


The transparency of the sea is generally stud- 
ied in situ by means of a transparency-meter 
which, when standardized (JERLOV, 1944), 
records the true light extinction, i.e., the sum 
of absorption and scattering. The extinction 
is caused by the water itself, by suspended 
particles, and by dissolved substances, chiefly 
yellow substance. Studying the extinction 
by unfiltered and Berkefeld-filtered samples 
of sea water it is possible to distinguish be- 
tween the extinction due to particles and that 
due to dissolved matter (James and BIRGE, 
1938; CLARKE and JAMES, 1939). This proce- 
dure is, however, too laborious to be used 
conveniently in routine work at sca. 

A great deal of interest has been directed to 
the particle component in the sea, in partic- 
ular from a biological point of view (see 
KREY, 1949; ARMSTRONG and ATKINS, 1950). 
In an earlier report the present author has 
described a method of studying the distribution 
of particles in the sea with the aid of a Tyndall- 
meter (JERLOV, 1951). The method offers the 
possibility of characterizing water masses and 
studying horizontal and vertical movements 
in the sea by means of the particle content, 
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which is expressed as a scattering coeflicient 
(conveniently in km~?). 

On the other hand, the transparency of the 
sea as measured by means of a transparency- 
meter can be used for identifying a water 
mass (JOSEPH, 1950; WYRTKI, 1950; JERLOV, 
1951; Krey, 1952). Alone transparency meas- 
urements, even if made for different colours, 
do not furnish sufficient data to distinguish 
between the extinction due to particles and 
that due to dissolved substances, as the ab- 
sorption by particles as well as that by yellow 
substance exerts a selective action on the trans- 
mitted light. 

It is obvious, however, that if the particle 
content is studied together with the selective 
light extinction the task of distinguishing 
between different components of extinction 
would be highly facilitated. An attempt to 
realize this is made below. 


Particle distribution 


Coastal waters, as a rule, contain more 
suspended particles than do ocean waters, 
due to a richer plankton population and a 
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larger amount of detritus and of inorganic 1951) whereas the Kattegatt and the Skagerrak 
matter derived from land or from the bottom. are seldom lower than 100 in particles. In 
The clearness of deep ocean samples is illus- order to define our problem of distinguishing 
trated by a particle content of 20 (JERLOV, between different components of light ex- 
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Fig. 1. Particle distribution in the Gullmar Fjord on the 18th July 1952. 
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Fig. 2. Particle distribution in the Gullmar Fjord on the 20th November 1952. 
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tinction it is pertinent to give some aspects 
on the particle distribution in the ‘stratified 
water at the Swedish west coast. 

A section from the Gullmar Fjord run on 
the 18th July 1952 is shown in fig. 1. On this 
occasion, after a period of strong northwest- 
ern winds, the fjord was filled with high- 
salinity water, the upper layers were well 
mixed, and the thermocline was as deep as at 
40 m. The particle content at the surface was 
low but increased toward the innermost part 
of the fjord. The main body of fjord water was 
characterized by a particle value of 200— 
300. Only near the bottom at the shallowest 
stations particles were more abundant, which 
is attributed to sediments being whirled up. 
Owing to this effect the clearest water is on 
an average present at mid-depths in the 
fjords (JERLOV and Koczy, 1947). 

Quite another situation, which sometimes 
occurs in the autumn or in winter in the 
Gullmar Fjord, is illustrated in fig. 2. This 
section was run in November 1950. The 
fjord was covered with a turbid layer of low 
salinity, due to river water and precipitation, 
and spread over the salt water under the pre- 
vailing calm conditions. The layer decreased 
in thickness toward the mouth of the fjord 
but its terrigenous material supplied the entire 
basin with an abundance of particles (the scale 
is here twice as large as in fig. 1) which, 
judging from the Tyndall observations, were 
definitely brown at all levels. 

As first observed by Hans PETTERSSON 
(1934) with the aid of optical instruments an 
accumulation of suspended matter is generally 
present in or near a density discontinuity 
layer in the sea. Such a distribution appears 
in fig. 1 at 40 m but more clearly in fig. 3 
for the international P-section, Pater Noster- 
Skagen (at the boundary between the Katte- 
gatt and the Skagerrak). In particular, the 
discontinuity, when close to bottom, is connect- 
ed ‚with an abundance of particles probably 
stirred up by vertical oscillations in the layer, 
due to the tides. 

The rivers carry considerable amounts of 
sediment, most of which flocculates and settles 
rapidly. An example of these conditions is 
given in fig. 8 showing the outflow from the 
Göta Alv. The penetration of salt water as a 
deep layer from outside into the Göteborg 
Harbour is also clear from this fig. The par- 
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Fig. 3. Particle distribution in the P-section (boundary 
between the Kattegatt and Skagerrak) on the 2sth July 
1952. 


ticle contamination is as high as 3,000 at the 
Free Port (Frihamnen) where the river enters. 
Further it is observed that sediment is whirled 
up from the bottom especially in the shallow 
waters of the inner harbour. Accordingly, 
the clearest water is present as a tongue which 
reaches bottom at the Free Port. 


Components of light extinction 


Adopting the following notations. 


c = total extinction coefficient 
Cy = extinction coefficient by the water itself 
a, = absorption coefficient by suspended par- 
ticles 
a, = absorption coefficient by yellow sub- 
stance 
s = scattering coefficient by suspended par- 
ticles 
We have 
C— Cy = dy + ay +5, (1) 


where ¢, is known (CLARKE and JAMES, 1939), 
c is measured by means of the transparency- 
meter, and s by means of the Tyndall-meter. 

As a matter of routine the extinction was 
determined with different colour filters, viz. 
RG ı (655 mu) and UG 1 + BG 12 (380 my), 
sometimes also with BG 12 (440 mu). The 
absorption of red light by yellow substance 
in concentrations found in the sea is practically 
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Fig. 4. Relationship between scattering by particles, s,. 


and selective absorption, Aa, + 4 


nil. Hence, we can write the difference A 
between coefficients for 380 mu and for 655 mu 
in the following way. 


A (c—cw) cel A ices: VA dp sé (ay) 380 mu (2) 


The factor As is small, if the particles are 
not strongly coloured, and does not exceed 
the experimental errors in the determination 
of c. Nevertheless the factor is taken into 
account below. 

The absorption of red light by particles 
comes out directly from equation (1) as ay = 0. 
Apparently it is proportional to the scattering 
by the particles but the values are too low to 
permit any definite conclusions. But there 
is every reason to believe, that the absorption, 
dy, for a given wave-length is a function of 
scattering, s, even if the relationship may be 
different for a different average size and 
colour of particles in the sea. We consider the 
difference Aa, and let the experiments de- 
cide whether Aa, is a function of s or not. 
Therefore, the left member in equation (2) 
containing the measured quantities should be 
plotted against the particle scattering s for 
depth intervals where the content of yellow 
substance is expected to be constant. 

The fact that yellow substance apparently is a 
stable soluble product points to it being a con- 
servative concentration. This could not hold 
true in a strict sense, as a production of yellow 
substance normally takes place wherever or- 
ganic matter in the sea decomposes to which 


should be added that a certain consumption 
would seem probable, but the assumption is 
fully acceptable for Swedish coastal waters 
which are rich in yellow substance largely 
supplied by fresh water outflow. The conse- 
quence would be that content of yellow 
substance and, for instance, salinity are inti- 
mately connected. KALLE (1949) has proved 
this relationship to exist by studing the yellow 
substance in filtered water-samples, collected 
in the Gulfs of Bothnia and of Finland. The 
work by JosepH (1950) and by Krey (1952) 
points in the same direction. 

Therefore, in water columns of constant 
salinity the yellow substance is likely to show 
small variations. Fig. 4 is prepared in view 
of this. At different localities in off-shore as 
well as in fjord waters, two levels are selected 
which display almost the same salinity but 
different particle content. Straight lines con- 
nect the pairs. Most of these depths are above 
the thermocline. Surface observations have 
been avoided as they are often effected by 
contaminations from the ship. 

The diagram in fig. 4 shows the selective 
absorption to be strongly influenced by par- 
ticle scattering. The factor A a, is responsible 
for this selective action. As the lines are prac- 
tically parallel, Aa, can be expressed by the 
relation 


Ad = 1.5.8 (3) 


This simple formula is found for the west- 
ern Swedish areas at the end of July 1952. 
It is not expected to be of general application 
as the relationship between absorption and 
scattering by particles is dependent on their 
average size and colour. 

This attempt to interpret the role of particle 
absorption enables us to find from equations 
(2) and (3) values of yellow substance ex- 
pressed as absorption coefficient for 380 my. 
Then, all components involved in the extinc- 
tion are obtained from equation (1), for ultra- 
violet light as well as for red light. 

The resolved components for 380 and for 
655 my and their variations with depth are 
set out in fig. 5. The long wave-lengths are 
to a higher degree scattered than absorbed, 
whereas the case is opposite for the short 
wave-lengths. Because of particle accumula- 
tion in the discontinuity layer below 30 m, 
all factors increase toward the bottom except 
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yellow substance, which shows a slight de- 
crease. This brings us over to the distribution 
of yellow substance. 


Vertical and horizontal distribution of 
yellow substance 


Another illustration of the fall off in con- 


centration of yellow substance in the discon- 
tinuity layer is given in Table ı for the inter- 


Table 1. Extinction components. Section PI—P4 
on the 25th July 1952 


Depth 655 mu 380 mu 
M = 
Et, Ss | 4, Cait S | oF a, 
0202 7210.277 | 0.04 | 0.56 | 6.18 | 0.28 |'0:10 
Io 19 T2 N 55 15 29 II 
20 20 15 5 57 18 2 10 
30 15 I4 il 49 16 2 10 
40 18 17 I 52 19 2 6 
50 26 24 2 67 26 39 3 


national P 1—P 4 section. The upper layers, 
which are mixed with Baltic water, display a 
yellow value of 0.10 which drops to 0.03 at 


Fig. 5. Extinction components for 380 my and for 
655 my in the Kattegatt (57° N; 12° E) on the 24th 
July 1951. Total extinction coefficient = c; extinction 
coefficient by the water itself = c,,; absorption coef- 


ficient by suspended particles — ap; absorption coef- 
ficient by yellow substance = a,; scattering coefficient 


by suspended particles = s. 
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Fig. 6. Decrease of the content of yellow substance in 
the upper layers, toward the north in the Transition 
Area. 


so m where North Sea water is present. 
The density increase between 30 and so m is 
associated with the usual richness in particles 
as already shown (fig. 3). 

The considerable difference between the 
Baltic and the Skagerrak as to the content of 
yellow substance is also demonstrated in fig. 6 
which refers to measurements at depths above 
the thermocline carried out at the end of 
July 1952. The Baltic water becomes rapidly 
intermixed in the Sound. There is a marked 
bending of the curve where the Sound opens 
out into the Kattegatt. The decrease of yellow 
substance is more gradual in the Kattegatt 
where an additional amount of yellow sub- 
stance is furnished by rivers. On the whole, 
the distribution of yellow substance charac- 
terizes the hydrographic situation, and the 
yellow colour identifies the Baltic water when 
it spreads toward the north. 

An intensification of the yellow colour is 
often a local phenomenon due to outflow of 
fresh water. Particularly in the fjords, where 
the low turbulence favours a very pronounced 
stratification, a layer of low salinity is fre- 
quently formed over salter water. This layer 
does not only contains an abundance of par- 
ticles but has also a high content of yellow 
substance. Values of a, = 1.0are often observed 
at Bornö Station. 
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Fig. 7. Salinity distribution at the outflow of the Göta Alv. 
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Fig. 8. Particle distribution at the outflow of the Göta Älv, 
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Fig. 9. Distribution of yellow substance at the outflow of the Géta Aly. 


When we discuss the role of rivers for the measurements were carried out for blue light 
colour of the sea it is worth while to ascertain as well as for red light, and the content of 
whether the yellow substance may be used as yellow substance is expressed as the absor- 
an indication of water from the Göta Aly. In tion coefficient for 440 mu. The values of a, 
view of the high extinction in the river, the found in the inner harbour are somewhat un- 
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certain as equation (3) possibly does not hold 
true for particles abundant in the river water 
(fig. 8). This is characterized by a yellow 
value of 0.8 which corresponds to about 1.6 
for 380 my. The penetration of sea water as a 
deep layer into the innermost harbour and 
the mixing boundary appear clearly in fig. 9. 
The inflowing deep current retains its colour 
remarkably well, whereas the yellow sub- 
stance in the outflowing river water is rapidly 
reduced seawards by the addition of salt water 
from below and by horizontal diffusion. At 
a distance of 20 km from the mouth of the 
river the content of yellow substance is not 
far from normal for the off shore water. 
The trend of the isopleths for yellow sub- 
stance and for salinity (fig. 7) are in good 
agreement. 


The distribution of yellow substance may 
be used as an indication of fresh water, in par- 
ticular where it flows out into sea water of low 
salinity which cannot be determined with 
accuracy. This discussion and also the German 
investigations emphasize the point that the 
yellow substance in certain cases may be con- 
sidered as a characteristic property of a water 
mass. It is observed that the particle scattering 
found by the Tyndall-meter needs to be known 
only in relative units. Transparency measure- 
ments in coastal waters should preferably be 
made with the aid of a transparency-meter 
with a light path of two metres and with an 
automatic device for changing colour filter. 
They do not require any additional time at 
the stations as they can be conducted at the 
same time as the routine serial observations. 


RÉPERENCES 


ARMSTRONG, F. A. J., and Arxins, W. R. G., 1950: The 
suspended matter of sea water. Journ. Marine Biol. As- 
606,29, p. 139. 

CLARKE, G. L., and JAMES, H. R., 1939: Laboratory ana- 
lysis of the selective absorption of light by sea water. 
Jour. Opt. Soc. Amer., 29, p. 43. 

JAMES, H. R., and Bircz, E. A., 1938: A laboratory 
study of the absorption of light by lake waters. Trans. 
Wisconsin Acad., 2I, 

JERLOV, N. G., 1951: Optical studies of ocean waters. 
Reports of the Swedish Deep Sea Expedition, 3, No. 1. 

— 1951: Optical measurement of particle distribution in 
the sea. Tellus, 3, p. 122. 

JERLOV, (JOHNSON) N. G., and Koczy, F., 1947: Re- 
cords of transparency in Gullmar Fjord. Conseil 
Perm. Internat. p. ? Explor. de la Mer, Jour. du Conseil, 
15, p. 39. 

JosepH, J., 1950: Durchsichtigkeitsregistrierungen als 
ozeanographische Untersuchungsmethode. Deutsche 
Hydrogr. ZS, 3, p. 70. 

— 1950: Quantitative Durchsichtigkeitsmessungen im 
Meere. Deutsche Hydrogr. ZS, 3, p. 214. 


Tellus V (1953), I 
5—205429 . 


Katie, K., 1939: Die chemischen Arbeiten auf der 
”’Meteor’’-Fahrt Januar bis Mai 1938. Ann d. Hydr. u. 
Mar. Met. Beiheft zum Januarheft, p. 23. 

— 1947: Der Stoffhaushalt des Meeres. Leipzig. 

Katte, K., 1949: Fluoreszenz und Gelbstoff in Bottni- 
schen und Finnischen Meerbusen. Deutsche Hydr. 
Zeits. 2, p- 117. 

Krey. J., 1949: Über Art und Menge des Seston im 
Meere. Verhandl. d. Deutschen Zoologen in Mainz, 
1949, Pp. 295. 

— 1952: Die Charakterisierung von Wasserkôrpern 
durch optische Messungen. Archiv f. Hydrobiol., 49, p.1. 

PETTERSSON, HANS, 1923: Göta Älvmynningens hydro- 
grafi. Ur Göteborgtraktens natur. p. 73. 

— 1934: Scattering and extinction of light in sea-water. 
Medd. Göteborgs Högskolas Oceanogr. Inst., 9. 

WYRTKI, K., 1950: Uber die Beziehungen zwischen 
Trübung und ozeanographischen Aufbau. Kieler 
Meeresforsch. 7. 


A Further Study on the Relation between the Jet Stream 


and 


Cyclone Formation 


HERBERT RIEHL? 
University of Chicago 


and 


SIDNEY TEWELES, TR? 
U.S. Weather Bureau, Washington, D.C. 


(Manuscript received October 5, 1952) 


Abstract 


Events at sea level and aloft over the United States are related to the approach of a speed 
maximum in the jet stream. In the left hand portion of the area downstream from the jet maxi- 
mum where air at the jet level is decelerating, frontogenesis, cyclogenesis, and spread of precipi- 
tation occur. Other indications of high level divergence to the left of the advancing jet maxi- 
mum are given by changes in the structure of a nearby cold dome. Insofar as the changes in the 
cold dome are precedent, they constitute a means of forecasting cyclogenesis. 


I. Introduction 


In a previous paper (RIEHL, 1948) the follow- 
ing comment concluded a study on the rela- 
tions between the jet streams in the high 
troposphere and cyclone formation: 

“It is suggested that the jet stream appearing 
in connection with a pattern of very long 
waves in the westerlies provides a mechanism 
for the initiation of cyclone development 
and an increase of wave number. But it should 
be emphasized again that only one species of 
cyclone formation has been considered here, 


1) Participated under research contracts between the 
Office of Naval Research and the University of Chicago. 
- ?) Participated while on assigment to the Advanced 
Forecasters Course sponsored by the Weather Bureau 
and the University of Chicago. 


namely, that associated with initial westerly 
flow aloft without pronounced streamline 
curvatuie. Nor is it suggested that the jet 
alone can create cyclones. It is evident that the 
jet is effective only if it is superimposed on a 
disturbance of the lower atmosphere. Cleariy, 
the amount of cyclonic development to be 
expected, depends in large measure on this 
factor. Therefore, when jet stream, long wave 
pattern, and low tropospheric disturbance coin- 
cide in a favorable sense, ensuing cyclone devel- 
opments will attain the greatest intensity.” 

It is the objective of this report to deter- 
mine more precisely some of the favorable 
circumstances mentioned. As in the previous 
paper, we shall do this with an example whose 
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Fig. 1a. Surface isobars (mb) and fronts, November 12, 
1951, 1830 GMT. Area of steady precipitation is shaded. 


7); 


Fig. ıc. Surface isobars and fronts, November 13, 1951, 
1830 GMT. 


salient features are typical of a large group 
of cases, though not of all cases. 

In the course of an experiment in fore- 
casting carried out jointly by the University 
of Chicago and the U.S. Weather Bureau we 
observed that strong cyclogenesis frequently 
followed the appearance of elongated and 
nearly isolated cold domes aloft as shown in 
figure 2. The observation of these domes as 
such is not new, and few synoptic meteorolo- 
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Fig. 1b. Surface isobars and fronts, November 13, 1951, 
0630 GMT. 


Fig. rd. Surface isobars and fronts, November 14, 1951, 
0630 GMT. 


gists would dispute that, in the absence of the 
domes, the ensuing surface events would be 
quite different. But the laws that determine 
the course of the surface developments in 
relation to the thermal and wind structure 
aloft are far from obvious. We propose to 
bring out some pertinent facts of the case to 
be discussed that may provide some suitable 
hints concerning the routes to be followed in 
the search for the correct laws. 
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Fig. 2a. Winds and isotherms 
(°C.) at soo mb, November 
12, 1951, 1500 GMT. Pro- 
gression of cold dome in 12- 
hourly steps, November 12, 
0300 GMT to November 14, 
0300 GMT, is marked by 
heavy line. Axes of minimum 
and maximum temperature 
gradient are marked by dashed 
lines. ‘*W’’ denotes warm 
center and ‘‘C’’ denotes cold 
center. On wind vectors a 
long barb denotes 10 knots, 
a short barb $ knots, and a 
heavy tringular barb 50 knots. 


Fig. 2b. Winds and isotherms 
at 500 mb, November 13, 1951, 
0300 GMT. Axes of mini- 
mum and maximum tempera- 
ture gradient are marked by 
dashed lines. “W’’ denotes 
warm center and ‘‘C’’ denotes 
cold center. On wind vectors 
a long barb denotes 10 knots, 
a short barb 5 knots, and a 
heavy triangular barb so knots. 
Heavy dot marks position of 
surface low pressure center. 
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Fig.2c. Winds and isotherms 
at soo mb, November 13, 
1951, 1500 GMT. Axes of 
minimum and maximum 
temperature gradient are 
marked by dashed lines. 
““W? denotes warm center 
and *C’’ denotes cold center. 
On wind vectors along barb 
denotes Ioknots, a short barb 
5 knots, and a heavy triangu- 
lar barb so knots. Heavy dot 
marks position of surface low 
| pressure center. 


| 2. Surface events 


1 The setting of the period studied, Novem- 
ber 12—14, 1951, in the hemispheric picture 
is at a time when, subsequent to several days 
à of a westerly circulation with weak amplitude 
(high index or Index Stage NII, Rrent et al. 
1952), a relative maximum in the westerlies 
4 was shifting northward into the higher lati- 
i tudes. Over the United States two typical 
high index troughs has passed eastward prior to 
{ November 11 when a large low pressure area 
X formed over the western half of the conti- 
{nent and when, with increasing southerly flow, 
a precipitation area spread rapidly from Texas 
{to the Great Lakes. 
On November 12, 1830 GMT (fig. 1a) we 
1 observe this large low pressure area east of the 
Rocky Mountains. The circulation is quite 
‘disorganized. There are several weak centers, 
land the frontal analysis is complex and rather 
luncertain. Organization of this diffuse pattern 
» begins within 12 hours (fig. 1b) and proceeds 
lrapidly around 1830 GMT on November 13 
(fig. 1c). The low pressure center that emerged 
"from Kansas travels slowly toward the Great 
Lakes on a path with strong counterclockwise 
’curvature. Deepening steadily, it moves to- 
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ward the NNW late on November 13 and 
finally becomes nearly stationary south of 
Lake Superior as a great vortex with central 
pressure near 975 mb (fig. 1d). 


3. Thermal structure of the middle tro- 
posphere. 


At soo mb, the isotherms over North Amer- 
ica initially exhibit the relatively unorganized 
pattern (fig. 2a), typical of “high index” 
conditions. Then there travels across the con- 
tinent a line of maximum spread of the iso- 
therms (minimum temperature gradient), fol- 
lowed by a well marked and nearly isolated 
elongated cold dome with a great isotherm 
concentration on its south side (for convenient 
reference on studies of such concentrations 
see PALMEN, 1948, PALMEN and NAGLER, 1948, 
and PALMÉN and Newron, 1948). 

These two features of the isotherm field 
propagate eastward at a mean rate of about 
25 knots. The surface low pressure area lies 
intermediate between them and initially is 
closer to the line of minimum temperature 
gradient. We observe three interesting facts. 

(1) The surface low organizes a great 
distance from the cold dome (figs. 2 b, c)—as 


0300 
Nov. 14 


1500 
Nov 13 


1500 
Nov. 12 


0300 
Nov 13 


Fig. 3a. Vertical cross-section of potential temperature 

(°A) following cold dome on track given in figure 2a. 

Heavy line denotes tropopause and shaded area indicates 

adiabatic layer. Note that isentropes in the stratosphere 
are drawn for intervals of 10 °A. 


much as 1,000 miles distant. Whatever the 
role of the dome — which is pre-existent — 
in helping to organize the vortex at the sur- 
face, it cannot be simple superposition by 
means of height falls aloft. Actually, the 12- 
hour height changes at 300 mb were zero 
over the area of largest surface pressure falls. 
A more subtle mechanism is required, capable 
of downstream transmission over long dis- 
tances. 

(2) The Low deepens not in the region 
where the Polar front at 500 mb, as seen from 
the soo-mb isotherm gradients, is strongest, 
but far downstream in an area where the s500- 
mb temperature gradient is weak and nearly 
uniform along an axis drawn through the 
low center at the time when organization be- 
gins (fig. 2b). It is not possible to regard the 
development as occurring within the zone of 
strongest baroclinity of the troposphere. In- 
deed, it has been our observation that only 
weak and stable frontal waves form when an 
intense frontal zone at soo mb extends across 
North America with nearly uniform strength. 

These observations corroborate the work of 
Ryp (1923, 1927), POGADE (1938) and SCHER- 
HAG (1948) who point out that deepening 
preferably occurs in the left hand portion of 
“delta” regions of the upper wind and ther- 
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Fig. 3b. Tephigram showing soundings near center of 

cold dome. Horizontal lines are isentropes, vertical lines 

isotherms, sloping lines isobars (mb) and dashed line 
shows the moist adiabatic lapse rate. 


mal field. SUTCLIFFE (1947) derived an ana- 
lytical expression suggesting cyclogenesis in 
regions where there is maximum advection 
of cyclonic thermal vorticity. 

(3) It is possible to track the cold dome over 
the two-day interval. Figure 2 a shows its 
path. The velocity of the dome is approxi- 
mately the same as that of the winds in its 
center. Since these winds varied little with 
height above the mountains, a vertical time 
section following the center of the dome 
(fig. 3a) shows substantially the same air 
parcels as time progresses. Construction of 
this section presented some difficulty since the 
center was never situated precisely over a 
radiosonde station. Nevertheless, we preferred 
to draw the section with use of the nearest 
sounding at each observation time rather than 
interpolate from analyses. We consider this 
procedure as most straightforward in this 
case, since at each sounding period there was 
at least one alternate station with precisely 
the same indications. 

The section shows a rise of the tropopause, | 
in this area defined as the top of the cold | 
dome, by 30—40 mb until November 13,, 
1500 GMT, when it levels off. In the upper! 
troposphere, the isentropes ascend through the? 
period and the lapse rate, initially moist adiabatic 
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Table 1. Ascent of isentropic surfaces in 
the cold dome. 


(0) A p (mb) 
290 —30 
292 —60 
294 ——90 
296 05 
298 —95 
300 —65 
305 —10 


(fig. 3b), steepens to become dry adiabatic. Table 1 
shows, in millibars, the amount of ascent of 
several isentropic surfaces during the 36 hours 
from November 12, 1500 GMT, to November 
14, 0300 GMT. 

The observed cooling, particularly in the 
lowest part of the stratosphere, is too large to 
be explained by non-conservative processes. 
It follows that the cold dome is rising through the 
period as may also be confirmed by noting 
the increase of area covered by the —30° C. 
isotherm on successive soo-mb charts. In 
view of the steep lapse rates, the changes in 
thickness between isentropic surfaces, which 
can be computed from table 1, may not be 
sufficiently accurate to warrant a diagram 
showing divergence and convergence as a 
function of height. But the table does suggest 
that in the troposphere there is a deep layer of 
gradually converging air, that this air ascends 
at rates with magnitude of about one centi- 
meter per second and that it is evacuated 
laterally in a narrow layer under the tropo- 
pause. 

The foregoing observations do not accord 
with the viewpoint (MARGULES, 1903) that 
cold domes must sink as a whole during sur- 
face cyclogenesis. As just seen, far from sub- 
siding, the cold dome center actually spreads up- 
ward during the period of deepening. We are not 
suggesting, of course, that sinking of cold 
domes does not take place in the great majority 
of cases. But we do obtain the impression that 
such sinking is not uniquely necessary for 
cyclonic development (SPAR, 1950) and that 


! the role of the dome in the cyclogenetic 


mechanism may be other, at least initially, 
than the simple sinking usually visualized. 

The point has come up here that the ascent 
of the dome takes place while it crosses the 
mountainous regions of western North Amer- 
ica. A possible “mountain effect” enters into 
almost all detailed synoptic studies that can 
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Fig. 4a. Contours of the 300-mb surface (100’s feet, 
first digit omitted) and 12-hour height changes (100’s 
feet), November 12, 1951, 1500 GMT. 
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Fig. 4 b. Contours of the 300-mb surface and 12-hour 
height changes, November 13, 1951, 0300 GMT. 


Fig. 4c. Contours of the 300-mb surface and 12-hour 
height changes, November 13, 1951, 1500 GMT. 


be made over the continent. We do not see 
how this can invalidate the inferences just 
drawn. The point is that the cyclone deep- 
ens while the cold dome ascends, irrespective 
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Fig. 5a. Isotachs at soo mb (knots), November 12, 1951, 
1500 GMT. Heavy lines mark jet stream axes. Tick 
marks along axis of principal current serve to identify 
the horizontal axis of fig. 10a. Precipitation areas are 
shaded. Dashed lines give 3-hour surface isallobars (mb). 


N 
ED) 
Fig. 5 b. Isotachs at 300 mb, November 13, 1951, 0300 
GMT. Heavy lines mark jet stream axes. Tick marks 


correspond to those along horizontal axis of fig. 10 b. 
Heavy dot marks position of surface low pressure center. 
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Fig. sc. Isotachs at 300 mb, November 13, 1951, 1500 
GMT. Heavy lines mark jet strem axis. Tick marks 
correspond to those along the horizontal axis of fig. roc. 
Heavy dot marks position of surface low pressure center. 
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of the reason for such ascent which is not a 
topic of investigation in this report. 


4. Structure of the Jet Stream 


We shall now investigate the structure of 
the high-tropospheric wind field as a possible 
connecting link. Figure 4 shows the 300-mb 
contours for the period and figure 5 the 300- 
mb isotach analyses, prepared with aid of the 
observed winds and with computations from 
the contour field. As is generally the case, 
this analysis is least certain in the area of 
strongest winds. We cannot claim that we 
know the strength of the wind maximum 
and the gradients around it with precision. 
However, we believe that the analysis repre- 
sents as fair an approximation to true condi- 
tions as can be secured with the available ob- 
servations. Many features of the wind field 
as analyzed accord closely with previous de- 
scriptions. The vorticity distribution near the 
maximum, for instance, is very similar to that 
computed by PALMÉN and NEWTON (1948). To 
the left of the maximum the absolute vorticity 
is between 2f and 3f, where f is the Coriolis 
parameter. To the right, the anticyclonic shear 
amounts to —1.5 f; but this is offset by a cur- 
vature term of +0.5 f, so that the total abso- 
lute vorticity 1s very nearly zero. 

On November 12, 1500 GMT, a strong jet 
maximum is situated near the West Coast 
of the United States, centered about 600 miles 
south of the cold dome. Farther downstream 
the organized jet decomposes into the “fingery” 
structure common in “delta” zones (RIEHL et 
al., 1952). This accords with the open and irreg- 
ular isotherm pattern at 500 mb (figs. 2 a—b) 
over central North America. We also ob- 
serve that in this area the soo-mb winds 
cross the isotherms at large angles ranging up 
to 90°, and that winds and isotherms are 
nearly parallel only in the zone of great iso- 
therm concentration in the west. 

On the subsequent maps, the jet maximum 
propagates eastward parallel to the cold dome 
and, in the mean, at the same rate. The leading 
edge of the jet maximum, initially situated a 
little west of 102° W. (fig. sa), reaches the 
states just south of the Great Lakes on Novem- 
ber 13, 1500 GMT (fig. 5 c). It is at the time of 
arrival of this leading edge that the surface cy- 
clone organizes on its left hand margin. We can 
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see this clearly from the surface maps and from 
the three-hourly surface pressure tendency 
centers entered with dashed lines in figures sb, c. 

In consequence, we can regard the downstream 
propagation of the leading edge of the western jet 
stream as a link connecting the cold dome in the 
west and the cyclone formation in the Central 
States. There may be other such links. If so, 
we have failed to notice them. But the cor- 
ollary evidence which follows suggests that 
the proposed connecting link is realistic. 


5. Advection of Vorticity 


As brought out in the literature of recent 
years (for reference see PALMÉN, 1948, and 
RıEHL et al., 1952), it is likely that the surface 
pressure falls in areas where advection of more 
cyclonic absolute vorticity takes place in the 
upper troposphere. This statement is based 
on mixed dynamic-empirical reasoning that 
horizontal mass divergence is occurring in 
regions where higher absolute vorticity is 
imported from upstream and that this mass 
divergence exceeds any compensating con- 
vergence in the lower levels. 

Although vorticity advection charts are not 
presented, we can deduce readily from figures 


2 4—5 that the relative vorticity decreases down- 


stream at 300 mb to the left of the jet stream 
axis over the area where the surface pressure 
is falling November 12—13, and that there- 
with higher vorticity is brought into this 
area by the wind. The curvature of the con- 
tours changes from cyclonic to anti-cyclonic 
as we go downstream at a particular time 
(cf. figs. 4c, 5c), and the cyclonic shear weakens 
as the gradient of the isotachs becomes less. 
The downstream variation of the Coriolis 
parameter is small and may be neglected. We 
see then that the advection of absolute vortic- 


ity has the requisite sign for pressure fall. 
| If the 300-mb surface may be taken as repre- 


sentative for the upper troposphere, as is 


| generally the case, it follows that the law stated 


initially holds for the present case. 


6. Distribution of Vertical Motion 


Although we have just mentioned that a 
good relationship appears to exist between the 
signs of surface pressure change and high 
‚level vorticity advection, the correlation be- 
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tween the magnitude of these quantities is poor, 
except for short period fluctuations. The 
surface pressure fall is a small residual between 
low level mass convergence and high level 
mass divergence. Mass continuity is provided 
by vertical motions which under the condi- 
tions described must be directed upward over 
areas of surface pressure fall. Such broad scale 
ascent should be reflected in the precipitation 
pattern, except perhaps in the lee of mountain 
ranges. We observe the heaviest precipitation 
near the threehourly surface pressure fall centers 
and this precipitation, in the mean, lies to the left 
of the 300-mb jet stream axis. Thus the relative 
geographic positions of jet stream core, region 
of surface pressure fall, and region of precipita- 
tion correspond to that observed by STARRETT, 
1949 and that demanded previously on the 
basis of dynamic reasoning (UNIVERSITY OF 
CHICAGO, 1947 and RiEHL, NORQUEST, and 
SUGG, 1952). 


7. Cross-Stream Circulation 


It is one of the assumptions in the derivation 
of the relation between mass divergence and 
advection of vorticity, stated carlier, that the 
flow and vorticity patterns move much more 
slowly than the wind and that we can neglect 
local changes compared to advective changes. 
The validity of this assumption is borne out 
by figures sa—c. Although the winds blow at 
(computed) speeds near 200 knots in the jet 
stream core, the pattern propagates at little 
more than 25 knots. The air very rapidly 
moves through the pattern, and it must suffer 
extreme deceleration when leaving the area of 
highest wind. (Earlier case studies, Wosus, 
1950 and TEWELES, 1950, have described this 
phenomenon.) This, following the equations 
of motion, is accomplished mainly by motion 
toward higher pressure, i.e., by a clockwise 
cross-stream circulation looking downstream 
along the jet core. To the left of the axis- this 
high level cross-stream circulation is likely to 
be associated with mass divergence, as we 
have already established for the present case 
from the vorticity advection pattern. We see 
the cumulative effects of this circulation in 
figure 6a, a vertical cross section taken nor- 
mal to the jet stream axis from Lake Charles, 
La. (LCH) northward to International Falls, 
Minn. (INL). The time is 0300 GMT, No- 
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Fig. 6a. Vertical cross-section of isentropes (solid lines, °A.) and isotachs (dashed lines, knots), 

from Little Rock, Ark. to International Falls, Minn., November 14, 1951, 0300 GMT. Heavy 

lines denote fronts and tropopauses. ‘‘W”’ stands for center of westerly current, and ‘‘E”’ for 
center of easterly current. 


vember 14—12 hours subsequent to the time 
of figure sc. The principal wind maximum still 
is upstream from the section, and the pertinent 
portion of the isotach pattern has not changed 
appreciably except for continued gradual 
downstream displacement of the jet center. 
Choice of the section shown, rather than an 
earlier one, was prompted by a very suitable 
station distribution and availability of data. 

On the section, which cuts through the 
forward edge of the cold dome near Omaha 
(OMA), we locate the core of the jet at 260 
mb, slightly north of Little Rock (LIT), and 
about 400 miles south of the deepest portion 
of the cold dome. The section is drawn so 
that we face upstream. The winds which blow 
out of the section toward us are decelerating 
near the jet core as mentioned. According to 
computations made from the equations of 


motion, the angle between contours and 
streamlines in this region must have attained 
20° to allow the observed deceleration of the 
air. From there northward, the rate of decele- 
ration must decrease and eventually become 
small near the cold dome since, as seen earlier, 
the dome moves with the speed of the winds. 
Even here, however, some divergence has 
been taking place as brought out earlier. 
Computing an approximate value of this 
divergence from figure 3a and table 1 for the 
layer between the isentropes 298° A. and 
305° A. with the formula 


1 dAp 


Ap dt. 


= — div, v 


we find that div, vs: 10 sec-1, a rather 
small value. 
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Fig. 6b. “Probable meridional and vertical displace- 
ments associated with intensification of zonal wind 
maximum. Arrows indicate direction of displacement 
of isentropes. Isentropes given by thin lines, tropopause 
by heavy line.” (From UNIVERSITY OF CHICAGO, 1947.) 
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If we combine all the evidence adduced— 
the deceleration of air near the jet core, the 
decrease of this deceleration toward the north, 
and the presence of divergence in the upper 
troposphere some distance north of the jet 
core as given by the cold dome computa- 
tion and the vorticity advection pattern—it 
becomes plausible to suggest that the configu- 
ration of the isentropes in figure 6a indicates 
in part the accumulated effect of vertical dis- 
placements upstream. The sense of the cross- 
stream circulation pattern (fig. 6b) would 
be that suggested in an earlier publication 
(UNIVERSITY OF CHICAGO, 1947). We should 
like to emphasize, however, that the arrows of 
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figure 6 b are not meant to suggest closed 
circulation orbits. If air evacuated laterally 
from the top of the cold dome were to pass 
to the other side of the jet core with the cross- 
stream circulation, it would first have to 
assume the very high vorticities of the zone 
just north of the jet center, then the very low 
vorticities to its south. Clearly, this is most 
unlikely. Besides, we note in the present case 
that the thick isentropic layers south of the 
jet lie between 324° and 340° A., those north 
of the cold dome between 312° and 320° A., 
whereas the potential temperature of the air 
evacuated from the cold dome ranges from 
298° to 305° A. It is more likely that the part 
of the jet stream core downstream from a maxi- 
mum is gradually displaced toward higher con- 
tours on any isobaric surface and this indeed is 
observed in many cases (SAWYER, 1950). 


8. Vertical Variation of Jet Stream Axis 


In the jet stream publication mentioned 
initially, much emphasis was placed on the 
marked reversal of temperature gradient 
across the jet stream axis above the level of 
strongest wind. It was shown that a band of 
warm air extends along this axis on the pole- 
ward side at 200 mb, and that a narrow band 
of very cold air parallels its equatorward 
margin (fig. 7). The axis itself was situated 


Fig. 7. “Isotherms (°C.) at 200 mb, January 28, 1947, 0300 GMT. Jet stream center at 
300 mb marked by heavy line.” (From RIEHI, 1948.) 
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Fig. 8. Isotherms (°C.) at 200 mb, November 13, 1951, 

1500 GMT. Jet stream center at 300 mb marked by 

heavy line. Tick marks along jet axis serve to identify 
[the horizontal axis of fig. 9. 


within the zone of strongest 200-mb tempera- 
ture gradient. Such a position is requisite if 
the geostrophic component of the wind is to 
decrease with height. Farther poleward and 
farther equatorward the temperature gradient 
again reversed, thereby proving that the 200- 
mb temperature field found in the jet stream 
zone could not have advective origin but that 
vertical motions as shown in figure 6 b had to 
account for its existence. PALMEN and NAGLER 
(1948) have reached the same conclusions. 

Figure 8 shows the 200-mb isotherms at the 
time of strongest deepening, November 13, 
1500 GMT. In several respects, this chart 
verifies the description of the 200-mb tempe- 
rature field given earlier. Relatively warm air is 
in evidence everywhere at the tropical margin 
of the chart, and from there the temperature 
decreases toward the jet axis. We find very 
warm air on the poleward side of the jet 
center (cf. fig. 5 c) and very cold air on its 
equatorward side. Indeed there is a suggestion, 
particularly north of the jet maximum, that 
the centers of greatest 200-mb temperature 
anomaly are closely associated with the area 
of highest wind. In view of the discussion of 
figures 6 a, b this suggestion appears quite 
reasonable. 

We also note some interesting differences 
between figures 7 and 8. Along the jet axis 
the 200-mb temperature varies much more in 
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the November 1951 than in the January 
1947 case. In fact, the temperature gradient 
reverses along the axis downstream from the 
zone of strongest winds. This is also brought 
out clearly in graphical form in figure 9. 
It follows that the level of strongest wind 
must rise downstream above the region of 
surface deepening. Southeast of Lake Michi- 
gan it reaches the 200-mb surface. Over the 
eastern Great Lakes, it must actually lie above 
200 mb. 

As the intense 200-mb temperature gradient 
north of the jet axis over the north central 
Plains is directed nearly parallel to the con- 
tours (not reproduced but similar to those 
of fig. 4c), and as the air even at some distance 
from the axis moves at a rate several times 
greater than the speed of the system, we can 
safely conclude that ascending motion is 
taking place through the 200-mb surface. 
Here, the schematic vertical motion picture 
of figure 6b cannot hold entirely. The ascent 
in the rain area of figure 5 c appears to extend 
to very high levels, a conclusion similar to 
that reached by FLEAGLE (1947) in his studies 
of upper troughs and ridges. South of the jet 
axis, the temperature gradients are much weak- 
er. We can only state that the region where 
subsidence compensating for the upward mass 
transport through the 200-mb surface takes 
place, is not completely delineated by the 
charts presented here. 

We can verify some of the conclusions 
drawn from figures 8—9 by constructing iso- 
tach cross-sections along the axis of the jet 
stream, a representation which we have not 
yet encountered in the literature. This is done 
in figures 10 a—c. The sections shown in these 
figures follow the axis of the jet stream, and 
tick marks refer to points along the jet at 


140 V300 
120 Knots 


J 60 


—> S$ 


Fig. 9. Wind speed at 300 mb (knots), and temperature 
gradient (°C.) at 200 mb taken over distance of 250 km 
normal to jet axis (heavy line) of fig. 8. Marks at 
bottom correspond to tick marks along jet axis of fig. 8. 
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300 mb as given in figures sa—c. It is to be 
noted that the sections do not portray condi- 
tions along the vertical but that they pass 
through the axes of strongest wind at all levels. 
For the construction we first analyzed isotach 
charts at 700, 500, 300, and 200 mb, and in 
part also at 250 mb. We then drew lines con- 
necting the points of highest speed on each 
surface and projected these lines to coincide 
with the 300-mb axis. Finally we plotted the 
wind values so oriented on cross-section paper 
and drew isotachs. The distance projected did 
not exceed 200 miles and from 400—5$00 mb 
upward the axes nearly coincided, i.e., the 
jet stream was almost vertical as is commonly 
the case. 

The first impression that one gets from 
figures 10 a—c is that they resemble sections 
taken normal to the upper westerlies. Varia- 
tions of wind structure along the current have 
the same order of magnitude in the present 
case as variations normal to the current, 
although the intense gradient of baroclinity 
seen in figure 6a of course is not present. 
Figures 10 a—c verify the upward displace- 
ment of the jet axis downstream from the 
wind center (marked J) as inferred above. It 

ill be an interesting problem for the future 
to draw corresponding cross-sections of isen- 
tropes and attempt to determine the actual 
rate of deceleration of the air particles. This 
deceleration would be less than indicated by 
constant pressure charts if the air ascends 
substantially in the region where the isotachs 
trend upward. 

Figures 10 a—c reveals another curious fea- 
ture. We have already noted on figures $a—c 
that a secondary jet maximum formed at 
300 mb downstream from the main maximum 
on November 13 and that it was this maximum 
which was most directly connected with the 
cyclone formation. The longitudinal sections 
indicate how this center began to form early 
on November 13 (fig. 10b) and then became a 
separate entity later on that day (fig. 10 c). 
Qualitatively, one gets the impression that an 
“impulse” becomes detached from the main 
maximum and propagates forward at a more 
rapid rate than the parent center. If future 
synoptic studies should. establish a general 
connection between such “impulses” and 
deepening, a new approach to the problem 
of the dynamics of cyclogenesis would indeed 
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be provided (cf. RIEHL and JENISTA, 1952). 
At this time, further speculation on this topic 
is not warranted. 


9. Conclusion 


In an attempt to enlarge on previous de- 
scriptions of the relations between cyclogenesis 
and the structure of the upper atmosphere, the 
following has been noted. 

(1) The developing surface pressure fall 
in the case studied could not be explained by 
simple superposition of an upper pressure fall. 

(2) The cyclone did not develop on an 
intense frontal zone but far downstream from 
this zone. 

(3) While the deepening progressed, the 
center of the cold dome upstream ascended. 
Therefore simple sinking of the cold air 
could not account for the deepening. 

(4) An intense jet stream maximum located 
to the right of the cold dome, looking down- 
stream, elongated rapidly from the central 
Rocky Mountain area toward the Great Lakes 
by means of sending forward an “impulse” 
which could be observed forming a new wind 
maximum on November 13. 

(5) Downstream from the main jet center 
the axis of strongest wind ascended to reach 
levels above 200 mb over the zone of deep- 
ening. The ascent of the axis was coupled 
with upward motion of the high tropospheric 
air to its left. 

The scheme of cross-stream circulation in 
the jet stream zone proposed earlier (UNIVER- 
SITY OF CHICAGO, 1947) is supported by the 
data in the vicinity of the main maximum. 
In the region where the jet axis ascends, some 
modification is required. Nevertheless, the 
major part of the precipitation area is observed 
to lie to the left of the axis, in accordance with 
the earlier findings. 

(6) At the arrival of the forward edge of 
the “impulse” above a pre-existing weak sur- 
face Low, the latter began to deepen strongly. 
Dynamically, this deepening could be related 
to the observed pattern of vorticity advection 
aloft. 
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Fig. 10a. Cross-section of 
isotachs (knots) along axis 
following jet stream core, 
November 12, 1951, 1500 
500! 120 GMT. Marks at bottom 
correspond to tick marks 
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Fig. 10b. Cross-section of isotachs, November 13, 1951, 0300 GMT. Marks at bottom correspond 
to tick marks of fig. 5 b. For details of construction of figure see text. ‘‘J’’ denotes jet stream center. 
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Fig. 10 c. Cross-section of isotachs, November 13, 1951, 1500 GMT. Marks at bottom correspond to 
tick marks of fig. sc. For details of construction of figure see text. “J” denotes jet stream center. 
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Application of Barotropic Tendency Equation 


to Medium-Range Forecasting 
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Abstract 


The simple barotropic vorticity equation has been used to compute height tendencies on 
a few s-day mean soo mb charts. The results are approximately as good as similar com- 
putations made on synoptic flow patterns despite the fact that the time interval between 
observed charts is 5 times as great. The order of magnitude of the computed tendencies 
may be considerably improved provided some account is taken of the effects of trans- 
formation of potential into kinetic energy, surface friction and non-adiabatic heating. 
A first approximation to these factors may be obtained by setting them equal to their 


climatological or normal values. 


A simple method for integrating the barotropic equation is derived upon assuming that 
the anomalies of the large-scale flow patterns may be approximated by a relatively small 


number of Rankine vortices. 


I. Introduction 


In the preparation of weather forecasts for 
time-intervals longer than one or two days, 
it is Customary to isolate in some way ‘the 
larger-scale features of the atmospheric circu- 
lation. This procedure is common to most of 
the present schools of long-range forecasting. 
One way to accomplish this is to prepare 
time-averaged or mean circulation patterns 
at selected levels above the earth’s surface. 
This averaging process not only reveals sys- 
tems having large space dimensions but also 
those having large time dimensions as well. 
Such mean patterns have very slow periods 
of evolution and may be extrapolated for 
longer periods in advance. 

It is beyond the scope of this article to 
discuss in detail the pros and cons of this 
particular approach to the problem of extended- 
range forecasting. The reader must be re- 


ferred to two sources of the extensive litera- 
ture on this subject (Namias and Crapp, 
1951; Namias, 1951). The purpose of the 
present paper is to report on further investiga- 
tions of the use of physical-dynamical methods 
in forecasting the evolution of mean patterns. 
In particular, tests have been made of the 
simple barotropic vorticity equation which 
already promises to be reasonably successful 
when applied to short-range forecasts using 
synoptic charts (BOLIN and CHARNEY, 1951; 
STAFF MEMBERS, UNIVERSITY OF STOCKHOLM, 
1952). A very great advantage of this method 
is that it can be carried out in an extremely 
short time by the use of high-speed electronic 
computing machines. 

In this paper, only the evolution of 5-day 
mean circulation patterns at the s00 mb 
level will be discussed. 
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2. The vorticity equation for time-averaged 
flow 


tionless) flow and for constant pressure sur- 
faces may be expressed in the following form 
(ELIASSEN, 1952): 


ag on 
— = —v:-Vn—- nV v—-w— (I 
on one 


| 
| 
| 
| 
The vorticity equation for laminar (fric- 


at 


Here Ë is the component of relative vorticity 


dv du 
normal to the pressure surface | — — — }; v, 
OX» Oy» 
the horizontal wind vector; 7, the absolute vor- 
ticity (+ f); f the Coriolis force; and w 
| the vertical “wind” normal to the pressure 


| surface (3). The symbol V refers to quasi- 


horizontal gradients as measured on the pres- 
sure surfaces. The only approximation made 
in this equation for laminar flow is to neglect 


Sow Ow 7 
the term in 7 and Er which represents the 
Xp 


Pp 
turning into the vertical of the absolute vor- 
ticity vector. 


In applying this equation to turbulent or 
time-averaged flow, a procedure is used exactly 
similar to that employed in developing the 
well-known Reynolds stresses. First, each term 
in (1) is averaged for the desired time interval. 
| The resulting arithmetic mean values may be 
represented symbolically by placing a “bar” 
! over each term. Now let each of the mete- 
orological quantities—for example the abso- 
| lute vorticity — be represented by its arithmetic 
) mean value (7) for the period of averaging 
i plus a departure from this mean (7): 


nt (2) 
| 


| Equation (2) is now substituted in the aver- 
| aged expression for (1), together with identical 
| equations for the other quantities, and use is 
made of the fact that in all cases the sum of 
» che departures from the mean (27) is zero. 
| This follows from the definition of the arith- 
+ metic mean. Then, after some algebraic mani- 
pores the following equation is obtained: 
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The bar over each term indicates that it is 
measured directly on the mean flow pattern. 
Thus it is found that the vorticity equation 
for laminar flow may be applied unchanged 
to mean flow provided that a term is added 
(the last term in (3) ) which represents local 
convergence of vorticity due to the average 
vorticity transport by the eddies which have 
been smoothed out of the mean flow. An 
additional approximation has been made, 


however. The term w’ au has been neglected, 


so that the last term in (3) represents the con- 
vergence of vorticity due to horizontal eddies 
only. It is assumed that in the mid-troposphere, 
where this equation is applied, the eddy 
transport of meteorological quantities takes 
place in quasi-horizontal surfaces. The third 
and fourth terms in (3) are determined in- 
directly by such physical processes as non- 
adiabatic heat sources and sinks, surface friction, 
flow over mountain barriers, and transforma- 
tions between potential and kinetic energy. 
Since (3) applies to any mean flow it is 
assumed to apply also to “normal” circulation 
patterns. These are defined as mean flow 
patterns for a very long number of years 
for one particular month or season. To dis- 
tinguish between “normal” and “mean” quan- 
tities the subscript N will be used to replace 
the bar when referring to normals. The 


quantity ed is usually quite small compared 
Ü 


to the other terms (less than 5 % of vw: V nn)» 
and so for convenience it will be omitted 
here. Thus, with some re-arrangement the 
vorticity equation for normal flow becomes: 


0 
—VN° Vin = INV VN + Wy ae u 
+ Vv) x (4) 


It is assumed, as in applying (3) to synoptic 
charts, that the first two terms may be ex- 
pressed by the height contour field of a con- 
stant pressure surface by making the following 
geostrophic approximations: 
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In order to keep the theoretical develop- 
ment within one section it is necessary to 
anticipate the results a little. It has been found 
in applying (3) to longer-interval mean charts 
that it is not possible to neglect the last three 
terms (as has been done on synoptic charts). 
However, to a first approximation these terms 
may be replaced by their normal equivalents, 
given by (4). Making this substitution, and 
replacing the relative vorticity in the first 
term of (3) by its geostrophic value (6) one 
obtains a very simple equation for the local 
height tendency: 


av = (— v . V 7) —(— vn: VAN) (7) 


The last term may be evaluated once and 
for all from existing normal flow patterns, 
and so in practice this equation is not much 
more difficult to apply than the simple “baro- 
tropic tendency equation” obtained by drop- 
ping this last term. 

The principal assumption involved in the 
derivation of (7) implies that the sum of the 
last three terms of (3) does not depart far from 
its normal value, while the vorticity advec- 
tion (first term on the right) may vary within 
wide limits. Much argument, both theoretical 
and practical, may be brought to bear against 
this hypothesis. For example Starr and his 
co-workers (LORENZ, 1952) have shown fairly 
conclusively that the eddy transport of mo- 
mentum, equivalent to the last term of (3), 
is not only important for maintaining the nor- 
mal zonal flow in the Northern Hemisphere 
but in determining its day-to-day changes as 
well. Nevertheless evidence will be presented 
in the present paper that such an hypothesis as 
that given above must be made if the simple 
barotropic formula is to be applied at all to 
longer-interval mean flow patterns, and also 
if it is to be applied to synoptic charts on 
a hemisphere scale. 


3. Tests of the barotropic tendency 
equation 


Equation (7) with the last term omitted is 
called the “barotropic tendency equation”. 
As has already been mentioned, it has been 
applied with fair success to synoptic weather 
charts, which represent mean flow patterns 


for time intervals of several hours. The first 
attempt to apply it to longer-interval mean 
patterns was made by Auserr (1951), who 
obtained some encouraging results in com- 
puting tendencies on a single 5-day mean 
700 mb chart. The next test was made by 
Berson (1952) using s-day mean 500 mb 
patterns for the summer season. To a certain 
extent the research reported on here repre- 
sents an addition to that of Berson, as I had 
the benefit of working very closely with him 
during the latter portion of his study. 

In the present study a preliminary test of 
the barotropic formula was made on three 
successive 5-day mean soo mb charts for 
November 1951 and for the eastern Atlantic 
and Europe (Area 2 of fig. 8). These were 
chosen for the same area and dates used in a 
previous investigation employing synoptic 
charts. (STAFF MEMBERS, UNIVERSITY OF STOCK- 
HOLM, 1952). Also, the same procedure was 
employed. The results indicate that the aver- 
age correlation coefticient between computed 
and observed tendencies (0.6) is about the 
same despite the fact that the “tendency in- 
terval” (interval of 2 or 3 days between 
successive mean charts) is 4 to 6 times that 
for the synoptic tests (12 hours). This indicates 
at once that further tests on mean patterns are 
justified. However, it was also found that the 
computed mean tendencies were about 40% 
too large, as compared to only 10 % for the 
synoptic tendencies. The explanation for this 
discrepancy probably lies in the three last 
terms of (3) which have been omitted from 
the calculations. One possible explanation is 
that the sum of these three terms remains 
constant (perhaps close to its normal value) 
while the vorticity advection term is known 
to decrease rapidly in magnitude with the 
length of the averaging interval. Considera- 
tions of this kind have led to the develop- 
ment of (7). 

Using these same cases, tendencies were also 
estimated from (7) using the normal vorticity 
advection term as well. While this showed a 
distinct improvement over the barotropic 
formula, especially with regard to the magni- 
tude of the tendencies, the test was considered 
unsatisfactory mainly because of the smallness 
of the chosen area. It was decided that a further 
test should be made using a winter situation 
for the entire Northern Hemisphere. 
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4. Procedure for evaluating the barotropic 
equation with normal vorticity advection 


In computing mean tendencies for the 
hemispheric case from the barotropic equation 
the procedure used is the same as that de- 
scribed by Born and CHARNEY (1951) with 
only minor variations. Briefly, this consists of 
analysing the mean contours on a polar- 
stereographic projection. Horizontal gradients 
in contour height used in estimating wind 
speeds and relative vorticities from (s) and 
(6), as well as in obtaining the mean vorticity 
advection, are estimated by the method of 
finite differences using a rectangular grid 
having a spacing whose absolute magnitude 
is 620 km at 50° N. In contrast to the study 
of Bolin and Charney the local vorticity and 
vorticity advection is obtained for every 
point in the grid instead of for every other 
point. The tendencies are then computed for 
each point by the method of “relaxation”. 
It should be noted in the example which 
follows that the tendencies so obtained usually 
have the opposite sign to that of the vorticity 
advection. This is due to the fact that the 
latter determines the curvature of the tendency 


profile. 


In making the computations some assump- 
tion must be made regarding the tendencies 
for the gridpoints at the boundaries of the 
chosen region. In this study it has been assumed 
that these tendencies are zero. 


The verification of the tendencies is also 
slightly different from that used in the synoptic 
study. They are compared to an observed 
“centered tendency’ which is defined as the 
difference in contour height between mean 
charts one day later than and one day pre- 
vious to that on which the computations are 
made. 


In computing mean tendencies from (7) it 
makes no difference theoretically whether 
the mean and normal vorticity advection 
terms are first added and then relaxed, or 
whether they are relaxed separately and the 
results then added. In practice it would be 
preferable to adopt the former method since 
then only one relaxation is necessary. However, 
in the present study the second procedure was 
chosen, since it was desirable to investigate 
separately the importance of the two terms. 
The results of relaxing the normal vorticity 
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advection yields a fictitious normal tendency 
which in general has no direct relation to the 


ot 
assumed that (7) represents an improvement 
of the simple barotropic formula, then it is 
clear that the error in the latter should be 
related directly to the relaxed normal advec- 
tion. For this reason this will be referred to as 
the “normal correction term”. 


true normal tendency (Fe). If it may be 


5. Test of tendency equation for a winter 
situation and for the Northern Hemisphere 


The case chosen was the s-day mean 500 
mb circulation pattern for Feb. 7—11, 1949, 
shown as the thin solid lines in fig. 1. These 
lines are drawn for every 300 feet and are 
labeled in 100’s of feet. This mean pattern 
represents the average of 5 0300 Z synoptic 
charts obtained from the printed historical 
map series (U.S. WEATHER BUREAU, 1949). 
The synoptic situation was one of rapid 
breakdown of a large blocking high over 
northern Europe with a return to strong 
zonal westerlies over Britain and Scandinavia. 
Equally rapid changes, but in the opposite 
direction, were taking place over the Pacific 
and North America, where a very strong 
belt of westerlies was giving way to an exten- 
sive ridge in the Gulf of Alaska and a deep 
trough over the western United States. A deep 
low over western Asia was filling slowly. 
Thus this situation was quite abnormal, 
especially over the eastern Pacific and North 
America, and should provide a severe test of 
(7). The observed 2-day height changes cen- 
tered on the mean chart are shown as the 
dashed lines in fig. 4. These are drawn for 
every 5 units and expressed in 10’s of feet 
per 12 hours. 

Before discussing the synoptic situation 
further it is desirable to consider first the 
normal chart for February together with the 
derived normal correction term. The normal 
contours for February, obtained from the 
latest available normal chart series (U. S. 
WEATHER BUREAU, 1952) are shown in fig. 2 as 
the thin solid lines. The reader is referred to 
the above mentioned publication for a com- 
plete discussion of these normal contours, but 
it must be mentioned here that they differ 
in some important respects from previously 
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Fig. 1. Five-day mean s00 mb chart for Feb. 7—11, 1949. Contours (thin solid lines) labeled 


in 100’s of feet and drawn for every 300 feet. Vorticity advection (dashed lines) in units which 
on relaxation give tendencies in 10’s of feet per 12 hours, and drawn for every 10 units. 


published monthly normals (JomTt METEORO- 
LOGICAL COMMITTEE, 1944; SCHERHAG, 1948) 
particularly in Canada and northern Asia, 
where upper air data have increased. signifi- 
cantly in recent years. Because of the limita- 
tions of the data with which these charts 
were constructed it is safe to say that these 
contours are most accurate over North America, 
the North Atlantic and Europe; of less value 
over the Pacific; and least reliable in the Asiatic 
region. The same remarks, with only slight 
modification, may be said to apply to the 
individual s-day mean chart shown in fig. 1. 

The normal vorticity advection term is 
shown as the dashed lines in fig. 2. The values 
shown on the chart, with isolines drawn for 


every Io units, are in such units that when 
relaxed they yield tendencies in 10’s of feet 
per 12 hours. It is natural to attempt to in- 
terpret this field in terms of the physical pro- 
cesses leading to the establishment of the 
normal patterns. However, the task is ren- 
dered extremely difficult for two principal 
reasons: there is no way to tell from (4) alone 
which of the three terms on the right hand 
side is responsible for the normal vorticity 
advection; and more important, the various 
physical processes (such as non-adiabatic 
heating and cooling and flow over moun- 
tains) may each contribute to all three of the 
terms. To illustrate the diversified interpre- 
tations that have been placed on the normal 
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Fig. 2. Normal chart for February. For explanation of units, see fig. I. 


vorticity advection term it may be mentioned 
that WIPPERMANN (1952) has considered it to 
be largely the result of surface friction, non- 
adiabatic heat sources and sinks, or forced 
flow over mountain barriers. BERSON (1952) 
has presented evidence indicating that the 
frictional transfer term is considerably larger 
than the others for long-term averaged flow. 
Finally, Namias and Crapp (1946) have in- 
terpreted the normal advection at 700 mbs 
in terms of fields of divergence which must 
accompany stationary waves in a baroclinic 
current. 

In spite of these difficulties the question is of 
great importance since its answer may deter- 
mine to a large extent the emphasis which 
should be placed on attempts to approximate a 
true 3-dimensional atmospheric model for 
numerical computations. To date such at- 
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tempts have been confined to introducing 
some of the aspects of transformation between 
potential and kinetic energy (PHILLIPS, 1951). 
For this reason a few remarks will be made 
regarding the broad features of fig. 2. 

There seems little doubt that the meridional 
bands of positive and negative vorticity ad- 
vection in western North America and the 
eastern Pacific are due to forced flow over 
the mountain ranges. For example, the ex- 
tensive area of positive advection just to the 
west of the continental divide in the United 
States is clearly due to lifting, which produces 
divergence and upward motion through a 
thick layer in the lower troposphere. Both of 
the first two terms to the right of (4) make 
positive contributions in this case. This opinion 
is strengthened when it is noted that the max- 
imum geostrophic wind speed normal to the 
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Tendencies from barotropic formula without normal vorticity advection (dashed lines) in 10’s 
of feet per 12 hours, drawn for every 5 units. 


mountain crests is located at about 40°N, 
just to the south of the maximum vorticity 
advection (+ 7 units in fig. 2). 

The prevailingly negative values over the 
oceans and positive values over the relatively 
flat continental areas are most likely due to 
heating and cooling from below. This is 
supported by the fact that the largest values 
are located near the principal heat and cold 
sources. Thus, maximum negative values are 
located east of Greenland, in the central por- 
tions of both oceans, and over the Mediter- 
ranean, Black, and Caspian Seas; while maxi- 
mum. positive values are located near the 
known cold sources in northern Canada and 
Siberia. However, the exact way in which 
these large-scale heat sources and sinks con- 


tribute to each of the three last terms in (4) 
is not yet understood. 

Many of the smaller-scale features of fig. 2 
are no doubt due to local uncertainties in the 
normal charts, although the reader’s attention 
is called to the geographical similarity of cer- 
tain features over the two oceans. 

Figure 6 shows the normal correction term, 
or relaxed values of the normal vorticity 
advection. This illustrates the very extensive 
smoothing which results from the relaxation 
procedure..Most of the minor details of fig. 2 
have disappeared leaving only the large-scale 
features which were discussed previously. The 
magnitude of the normal correction term is 
quite large. Comparison with fig. 4 shows 
that it is of the same order as the observed 
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Fig. 4. Five-day mean chart for Feb. 7—11, 


1949. Contours (thin solid lines) and observed 


DR en er 


two-day centered height change (dashed lines). For explanation of units see fig. 3. 


tendencies. It is approximately one fourth 
the size of 12-hour height changes observed 
on synoptic charts. The reason for the great 
smoothing and large tendencies is due of 
course to the extensive fields of predominantly 
positive or negative vorticity advection over 
the continents and oceans. Since in accord- 
ance with (6) the vorticity advection gives 
the curvature of the tendency rather than the 
tendency itself, it is clear that uniform fields 
of small values of the advective term will 
produce very large tendencies in the centers 
of these fields. Such extensive fields also lead 
to great difficulties in completing the relaxa- 
tion procedure, as anyone familiar with this 
process will know. In the present case time 
limitations prevented complete relaxation of 
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fig. 2, but this cannot affect the sign and gen- 
eral order of magnitude of the tendencies. 
The advective fields obtained from the s-day 
mean hemisphere chart of Feb. 7—11, 1949 
are shown in fig. 1. In its large-scale features 
and even in some of the minor details the 
pattern bears much resemblance to the nor- 
mal fields of fig. 2. The values are mainly 
negative over ocean areas and positive over 
land areas. However, the positive values due 
to lifting over the Rockies are much larger 
than in fig. 2, due of course to the much 
greater wind speeds. It is of interest to note 
that the largest value in this area (+ 41) is 
located a little to the north of the maximum 
wind normal to the mountain crests, which is 
again located near 40° N. This feature is in 
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Fig. s. Error in tendencies computed from barotropic formula without normal vorticity advec- 
tion, defined as computed minus observed tendencies. See fig. 3 for explanation of units. 


general agreement with the distribution of 
geostrophic winds normal to the mountain 
ranges. These increase rapidly to their peak 
at 40° N but then appear to level off into a 
broad maximum with but slowly decreasing 
speeds until they decrease rapidly in magni- 
tude north of 47° N. While the maximum 
s-day mean winds in this area are about twice 
their normal value the mean vorticity advection 
is almost six times as large as normal. 

It is not surprising from what has been said 
that the computed barotropic tendencies (fig. 
3) obtained by relaxing the field of fig. ı 
should bear a marked resemblance in many 
areas to the normal correction term of fig. 6. 
However, over the Rocky Mountains the 
mean tendencies are 7 times as large as their 


normal counterparts, and the accompanying 
large field of negative tendencies seems to 
spread out over a large area on either side of 
the ranges. 

Comparing these computed tendencies with 
those observed (fig. 4) it will be seen that 
over Europe and the Atlantic the general 
pattern is in fair agreement. However, as found 
in the preliminary test, the values are in 
a general much too large, in some cases ex- 
ceeding the observed by several hundred per 
cent. Over the relatively flat areas in the 
eastern United States and western Atlantic 
the agreement is not very good, although the 
rises in the central United States are estimated 
fairly well. The worst error in this region is 
the failure to predict the broad area of falls 
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Fig. 6. February normal correction term, defined as relaxation of normal vorticity advection. 
For explanation of units, see fig. 3. 


extending from northeastern Canada to the 
Antilles, with extremes east of Hatteras and 
southeast of Hudson Bay. However, near these 
two areas there are weak relative minima in 
the field of computed positive tendencies. 
In the western United States and extreme 
eastern Pacific the general pattern of com- 
puted tendencies is fair although the negative 
values are all several times too large. It is 
interesting to note that an observed fall area 
is located near the minimum in the computed 
field. The marked rises which occurred in the 
Gulf of Alaska were not computed, in spite 
of the fact that fig. ı shows a rather broad 
area of negative vorticity advection in this 
region. This result is due to the profound 
influence of the mountain effect, which upon 
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relaxation “spread” to the east and west and 
produced large errors in adjacent regions. 

Over the Pacific the agreement is in general 
unsatisfactory although the eastward motion 
of the relatively sharp trough east of northern 
Asia is qualitatively anticipated. The most 
conspicuous error in this area is the failure of 
the “cut-off” low northwest of the Hawaiian 
Islands to move north-eastward as expected. 
This may be due in part to poor analysis 
since during the period in question there 
were only three upper-air stations over this 
broad area of the Pacific, and these were to 
the south of the low center. 

Over Asia the agreement is poor, although 
there is some indication of the slow filling of 
the deep vortex in western Russia. 
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As expected from (7) the pattern of the 
errors in the computed tendencies (fig. 5) is 
in fair agreement with the normal correction 
term (fig. 6) over the areas of plentiful data 
from eastern Europe to central North Amer- 
ica. The magnitude of the errors is also in 
fair agreement although somewhat too large. 
Over the Rocky Mountains where data is 
also plentiful, the agreement in pattern is 
again fair, although the magnitudes are poor. 
Some agreement in general pattern can be 
noted in the Pacific area while over Asia, 
where the data is least reliable, the results are 
poor. 

As a final test of the method, the complete 
equation (7) was evaluated (fig. 7) by sub- 
tracting the field of fig. 6 from that of fig. 3. 


. Tendencies from barotropic formula with “normal correction term’’ added. For explana- 
tion of units, see fig. 3. 


Figure 7 should be compared to the observed 
tendencies (fig. 4). The marked improvement 
over Europe and the Atlantic can be noted, 
although the order of magnitude of the 
computed tendencies, while greatly improved, 
is still somewhat too large. Over North Amer- 
ica and the western Atlantic there is again 
some improvement, especially in the magni- 
tudes, although the negative tendencies in 
eastern North America are still lacking. How- 
ever, there is some sharpening of the relative 
minima close to those observed. The errors 
in the Rocky mountain region are still quite 
large, and no marked improvement can be 
noted over the Pacific and Asia. 

These results have been summarized statis- 
tically in Table 1. The Northern Hemisphere 
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Fig. 8. Areas of Northern Hemisphere used in com- 
puting statistics of Table ı. 


was divided into 7 areas, shown in fig. 8. 
Statistics were computed for only the middle 
latitude regions (Areas 2 to 6). The bounda- 
ries between the areas were chosen more or 
less arbitrarily although the northern and 
southern boundaries as well as those defining 
the Asiatic region were chosen on the basis 
of the distribution of available data. The sta- 
tistics shown in the table are defined as follows: 
X, in the first and second columns under each 
area refers to the observed centered tendency, 
while in the last column it represents the nor- 
mal correction term; X, in the first to third 
columns is respectively the computed baro- 
tropic tendency, the ‘“corrected” tendency 
from (7) and the error in the barotropic tend- 
ency; Ox and Gx, are the corresponding 


standard deviations; and X, and X,, the arith- 
metic means. The correlation coefficient is an 
indication of how well the observed tendencies 
may be estimated from those computed, 
while the standard deviations are a measure 
of the gross error or “‘bias” of the computed 
values. 

In general the numbers in Table 1 confirm 
the qualitative conclusions reached in previous 
paragraphs. In areas 2 and 3 the results agree 
with what one would expect from (7) while 
in areas 5 and 6 they do not. This result is 
certainly in agreement with what one would 
expect from the distribution of data over the 
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globe. It is interesting to note that although 
fair correlations are obtained relating errors to 
normal correction term in areas 2 to 4 there 
is no corresponding improvement in correla- 
tion for the corrected tendencies. However, 
in areas 2 and 3 there is marked improvement 
in the order of magnitude of computed tend- 
encies. 

The results in eastern North America are 
rather disappointing in view of the fact that 
this is a relatively flat region with plentiful 
data. However, in all probability this is an 
exceptionally poor case because tests on 
synoptic charts over the same area (BOLIN 
and CHARNEY, 1951) average somewhat better 
than those over Europe (STAFF MEMBERS, UNI- 
VERSITY OF STOCKHOLM, 1952), while all tests 
on mean charts made so far indicate about the 
same degree of skill as corresponding daily 
computations for the latter area. 


6. Short-cut methods of computation 


Solution of (7) by the method of relaxation 
is of course designed for use with high speed 
electronic computers. Because of difficulties 
such as the thousands of arithmetical opera- 
tions (each of which must be carefully checked), 
and the slowness of the relaxation procedure 
when dealing with large areas of uniform vor- 
ticity advection, the length of time necessary 
to carry out a hemispheric analysis by hand, is 
prohibitive except for very limited research 
studies. It is out of the question for use in 
practical forecasting. For this reason, and be- 
cause electronic machines will probably re- 
main generally unavailable for several years, 
it is desirable to develop short-cut methods for 
solving (7). Fortunately several such methods 
have been proposed in recent months. One of 
these, developed in connection with the pres- 
ent study, will now be briefly outlined. 

Let us designate the arithmetic mean value 
of any meteorological quantity, say absolute 
vorticity (7) by its normal value (74) plus a 
departure from normal (7°): 


D NN tn (8) 


Substitution of (8) in the first term on the 
right of (7) yields four new terms, one of 
which cancels with the normal advective term 
in (7). The following equation results: 
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LE =v". 79° vy: nv nn 
f Ot (0) 

The first term on the right of (9) represents 
the advection of departure from normal (d. n. 
vortex elements by the d.n. flow field. The 
practical solution of this term has been dis- 


cussed in a previous paper (CLAPP, 1952) where - 


it was assumed that one may represent the 
d.n. flow pattern as a collection of simple 
symmetrical Rankine vortices. A glance at 
the 5-day mean d.n. contour patterns shown 
in the above paper will suggest that this as- 
sumption seems quite reasonable. It was shown 
that based on this assumption the displacement 
of one individual vortex center can be ob- 


tained in Is to 20 minutes. This represents 
the local solution of the first term on the right 
of (9). The second term represents the advec- 
tion of the d.n. vortex centers by the normal 
flow field. This term is solved locally by 
simply computing the normal wind speed and 
direction at the center of a given d.n. vortex. 
Since the normal winds are zonal in character 
this will give displacements mainly towards 
the east. However, as pointed out (Crapp, 
1952) these displacements will be about two 
times as large as those observed as well as 
those computed from the first term on the 
right of (9). The reason for this becomes 
clear on considering the last term in (9) which 
represents the advection of the normal vorti- 


Table 1. Relation between 5-day mean computed and observed tendencies, and between estimated 

and observed errors in tendency. Units: 10’s of feet per 12 hours. Number of cases indicated by 

n. X, = observed tendency in columns 1 and 2 under each area, and ‘‘normal correction term’’ in 

column 3. X,=respectively computed barotropic tendency, ‘‘corrected’’ tendency from (7) and error 
in barotropic tendency. 


Area 2: Eastern Atlantic and Europe. 


Area 3: Eastern North America and western Atlantic. 


n= SI Be di 
Computed | Corrected Error Computed | Corrected Error 
VS. VS. VS. VS. VS. VS. 
observed observed | correction observed observed | correction 
tendency tendency term tendency tendency term 
RAT 10, oe Ors 0.67 0.73 0.86 XX 2 wenns — 0.20 — 0.12 0.84 
OPA DE ae 13 9 Io On ee II 6 12 
OR oy mover 6 6 6 Oe Bit San 3 3 6 
Ox, | Oxy oe AP 1.5 1.7 Ox, | Ox, ae, 1.8 RO 
Xp eres eees 4 fe) 4 owner Io 8 10 
Xie ee. Be o o 4 RIO EN I | I | 3 
Area 4: Mountain region Area s: Pacific. 
n = 30 i = So) 
TT oon 0.61 0.58 0.68 Pro anes 2a | 0.72 0.14 0.12 
Cats Da 12 9 10 Ol e 8 9 8 
Oase 4 4 3 OT ee oR 5 | 5 6 
Ox2 | Ox; 2.8 2.3 2.9 0%/[Ox..... 1.6 2.0 1.4 
Xp... Sag Fee fot EX. cee. PRE 2 5 3 
CAM Se 5 5 — I ee aT — it 4 
Area 6: Asia. 
n 67 
TXiXg sve ste ee — 0.13 — 0.33 — 0.46 
RO rn 5 7 7 
EN paire 4 4 3 
Ox, | Oxy ER 143 1.7 22 
Xac........ — 3 8 en 
MOON SORTE — I — I on 
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city field by the d.n. flow pattern. If the as- 
sumption is again made that the d.n. vortex 
centers are simple Rankine vortices, and if it 
may be assumed that 7, & f then a very 
simple expression can be obtained for the 
displacement of an individual vortex due to 
the last term in (9). This expression will 
contain B, the north-south gradient of f, as 
well as the dimensions of the system; and will 
produce a displacement towards the west. It 


—_ BL? 
is equivalent to the dynamic term in 
470? 


Rossby’s well-known wave equation (RossBy, 
1939), while the second term to the right of 
(9) represents the advective term (U). It is 
known that for large-scale systems these two 
terms nearly balance one another. 

In the manner outlined above a forecast may 
be obtained for the displacement of a single 
d.n. vortex center in approximately 25 min- 
utes. Since these centers often cover a very 
large area of the globe, such a forecast if 
successful would be of material value to the 
forecaster. 


7. Conclusions 


The basic assumption in the derivation of 
(7) and (9) is that the main role of terrestrial 
heat sources and sinks and the transformations 
between potential and kinetic energy is to 
maintain the normal circulation against the 
dissipative effects of friction. Otherwise the 
motion and development of the large-scale 
atmospheric systems is determined largely by 
the internal redistribution of vorticity within 
the atmosphere itself. It is clear that if this 
hypothesis is even partly true then the way 
is open to forecasting much longer weather 
changes, such as secular trends. Indeed, upon 


making similar assumptions STEWART (1943) 
showed theoretically that large-scale vortex 
systems having the general character of the 
atmospheric centers of action could have very 
long periods of evolution, perhaps several 
hundred years. 

It is not difficult to extend this hypothesis to 
include variations from normal of the physical 
processes mentioned at the beginning of the 
last paragraph provided it can be assumed 
that these processes are a function of the cir- 
culation itself. To give a simple example of 
this, it is possible to account entirely for the 
“mountain effect” since this is a function of 
such factors as the wind velocity and direction 
relative to mountain crests. 

It has been mentioned that the large areas 
of small but uniform values of the normal 
vorticity advection, which seem to act as a 
“background” to the much larger but less 
regular values of the individual advection, 
have a very important effect on the computed 
tendencies when these are carried out on a 
hemisphere-wide scale. Prior tests of the 
barotropic formula, especially with respect 
to synoptic circulation patterns, have been 
made for relatively small areas. If this formula 
is to be applied to synoptic charts on a hemi- 
spheric scale, an equation similar to (7) must 


be used. 
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On Temperature Frequency Distribution in the Free Atmosphere 


and a Proposed Model for Frontal Analysis 


By ROY BERGGREN, University of Uppsalat 
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Abstract 


The possibility of tracing baroclinic fields in the lower stratosphere by means of temperature 
frequency distributions in the main isobaric surfaces is discussed. The existence of a stratospheric 


front is indicated by the diagrams. 


In the study of the temperature frequency 
distribution in the main isobaric surfaces of the 
troposphere there is the opportunity to trace 
baroclinic fields. Starting with the sloping 
baroclinic field connected with frontal layers, 
the importance of which BERGERON and Swo- 
BODA pointed out already 1924, we may pursue 
the following reasoning. 

Let us begin with the simplest case, that of 
the whole temperature difference between 
two air masses being concentrated in the frontal 
layer (sec fig. 1). If this front is oscillating over a 
station, the temperature in a certain isobaric 
surface will be constant within the warm air as 
well as within the cold air. In the frontal 
layer, however, the station will show rapidly 
changing temperature. The temperature dis- 
tribution in an isobaric surface will therefore 
have two maxima, corresponding to the 
temperatures of the warm and cold air respec- 
tively. The temperatures between these two 
maxima will be less frequent. In the real 
atmosphere with its baroclinity also outside 
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the frontal layer, the above-mentioned picture 
will be less distinctive, but will certainly keep 
its main features. 

In the case of a sloping large-scale front in 
the atmosphere the warmer air is always 
situated above the front. This means that if a 
station has warm air in, say, the 700 mb sur- 
face, it must be in the warm air in the soo mb 
surface, too. This gives the same temperature 
frequencies for the two surfaces. In addition 
to these cases, however, we have cold air at 
the lower surface and warm air at the upper 
surface in all the cases, when the front is 
situated between these two isobaric surfaces. 
Thus, with increasing height the frequency 
of warm air will increase at the expense of 
the frequency of the cold air. 

According to McIntyre (1950, p. 102) the 
reasoning in the above paragraphs gives the 
necessary and sufficient conditions for the 
existence of a sloping baroclinic layer in the 
atmosphere. This cannot be quite true, be- 
cause the case of a horizontal inversion swing- 
ing up and down above a station would 
give the same characteristic temperature fre- 
quency distribution. From an empirical study 


96 ROY BERGGREN 


mb 
100 


= 150 


150 - 


200 - 


- 200 


250 - 


+ 250 


300 + 


+ 300 


400 4 


-- 400 


Tropical air 
wa L 500 


L 600 


lise te 


NZ 
Polar air 


IF 700 


+ 800 


+ 900 


Ee | 


Schematic vertical cross section based on two mean soundings (representative for western Europe and 


the eastern part of the middle Atlantic in November). Isotherms are given as thin solid (C°, intermediate 
values are given by dotted lines), isentropes as thin dashed lines (°K), and frontal boundaries and tropopauses 
as heavy lines. The inclination of the front is 1/100. (After BERGGREN, 1952 b, fig. 22.) 


of several stations at different latitudes and 
their relative frequencies of, say, Tropical 
air one can, however, state that the baroclinic 
layer must be sloping. From synoptic experi- 
ence and from theory it is obvious. 

The discussion of the temperature distribu- 
tion in isobaric surfaces can also be applied to 
the nearly vertical baroclinic field in the 
middle and upper troposphere connected with 
the jet stream. This baroclinic field appears 
as the strong baroclinity in the warm air im- 
mediately above the frontal layer (see e. g. 
BERGGREN, 1952 a, fig. 6). 

Up to this point a review of the above- 
mentioned paper by McIntyre (1950)has been 
given. When the present author’s attention 
was drawn to this paper, he was investigating 
some problems connected with the frontal anal- 


ysis in the upper troposphere and the lower 
stratosphere. In two papers (1952a and 1952b) 
I have tried to show that the tropospheric 
Polar front in many cases may be prolonged 
into the stratosphere. Fig. 1 in the present 
paper gives a schematic vertical cross section 
through such a front. If this idea is correct it 
ought to be possible to apply McIntyre’s 
reasoning to the frequency distributions above 
the 300 mb level. 

As a first example the British radiosonde 
station Lerwick on the Shetland Islands has 
been chosen. In order to compare the two 
papers, the period used in this study is the 
same as that used by McIntyre. From the 
Upper Air Section of the Daily Weather Report 
of the Meteorological Office, London, five winter 
months were chosen, viz. January, February, 
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Fig. 2. Winter temperature frequency polygons for 
Lerwick (60° N, o1° W). The dashed line gives the 
| mean values for the different levels. 


and December 1947, January and February 
1948. Furthermore only the oo GCT or 06GCT 
soundings have been used. In Mclntyre’s 
paper it is not clearly pointed out, whether 
the investigation is restricted to the night or 
carly morning soundings or not. In our case, 
however, we prefer to do so since we are 
studying stratospheric conditions, and in that 
case the radiation error may cause trouble. 
The frequencies in the 700, 500, 400 and 
300 mb surfaces are given in McIntyre’s paper 
as fig. 7. In our fig. 2 the temperature frequency 
distribution of the 250, 200, 170, and 150 mb 
levels are shown. The 250 mb distribution 
agrees with the 300 mb distribution given by 
MclIntyre, who relates the fact that this curve 
is nearly Gaussian to the observed fact that at 
about the 250 mb level there are almost baro- 
tropic conditions along a meridian. In my 
papers mentioned above I have pointed out 
that the temperatures of the Tropical and 
Polar air are about the same at this level. 


The fact that the frequency distribution for the 
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Fig. 3. Winter temperature frequency polygons for 
Aldergrove (55°N, 06° W). The dashed line gives the 
mean values for the different levels. 


250 and 300 mb levels shows a marked single maxi- 
mum, does not imply that we have only one air 
mass at this level. The maxima of the two 
distributions become so pronounced because 
the Polar and Tropical air acquire about the same 
temperatures at these levels. 

If we continue to the 200 mb level (see fig. 2) 
we get the double maxima once more. We 
are now well in the stratosphere as far as the 
Polar air is concerned but near the tropopause 
level in the Tropical air. The 200 mb surface 
is well above the level where the front changes 
its character from warm front to cold front 
or vice versa (see fig. 1). 

The Polar air maximum at —54° C and 
the Tropical air maximum at —61°C are 
easily detectable. If the reasoning in the first 
paragraphs is applied to this case, and if the 
front model given in fig. 1 is right, then the 
frequency distributions of the 170 mb and 
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the 150 mb surfaces must show an increase 
of the Polar air peak at the expense of the 
Tropical air peak. Fig. 2 shows that this is 
really the case. (A discussion of the reality 
of the double maxima will be given below.) 

Before going further I must strongly 
emphasize that it is not my intention to postu- 
late that this high-tropospheric and low- 
stratospheric front, seperating tropospheric 
Tropical air and stratospheric Polar air, always 
can be found. A short discussion of the applica- 
bility of the proposed model is given in the 
above-mentioned paper (BERGGREN, 1952b). 

The frequency diagrams for the British 
radiosonde station Aldergrove in northern 
Ireland were also constructed. In fact, they 
give the same picture as the one from Lerwick, 
and they are reproduced as fig. 3 without any 
special remark. Especially interesting are the 
extreme low temperatures down to: —70° C 
or even lower, attained in the Lerwick fre- 
quency diagram for 150 mb, but the dis- 
cussion of this item must be considered outside 
the scope of this paper. 

In order to be able to compare the results 
with those of McIntyre’s I have tried to choose 
the same period and the same stations as he 
did. But unfortunately the North American 
radiosonde material for the winter months of 
1947 and 1948 were not available here. For 
this reason five other winter months have been 
chosen. The data from December 1948 and 
January 1951 have been taken from Daily 
Upper Air Bulletin, the winter months of Jan- 
uary, February, and December 1949 from 
Daily Series Synoptic Weather Maps. 

In fig. 4 the temperature frequency diagrams 
for the 700, 500, 300, 200, and 150 mb sur- 
faces for the radiosonde station Tatoosh Island 
in northwestern USA are given. 

Before going further we must compare 
the mean values for the winter months in 


Table I. 
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=20) -10° ORG 
Fig. 4. Winter temperature frequency polygons for 
Tatoosh Island (48°N, 124°W). The dashed line gives 
the mean values for the different levels. 


question. These are given in table I. The fre- 
quency diagrams given by McIntyre for Ta- 
toosh Island have two maxima in the 700 mb 
diagram, corresponding to —2,5° C and 


The mean temperatures (°C) of different isobaric surfaces and for selected months for 


Tatoosh Island (48° N, 125° W). (From Monthly Weather Review and Climatological Data.) 
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jan. — 47|feb. —— 47|dec. — 47|jan. — 48|feb. —48 dec.—48|jan. — 49|feb. — 49|dec. — go]jan. — 51 
700 mb} — 9.2 | — 5.1 | — 8.6 | — 7.2 | — 10.9 — 10.3 | —12.0,| —12.8 | —10.2 | —11.9 
500 mb} — 23.9 | — 21.0 | — 24.0 | — 22.6 | — 26.6 || — 24.9 | — 26.9 | — 28.8 — 25.4 | 26.6 
300 mb] — 47.2 | — 46.7 | — 47.5 | — 46.8 | — 48.7 | — 48.1 — 49.9 | — 49.9 | — 47.7 | — 50.1 
200 mb] — 56.0 | — 61.3 | — 55.3 | — 53.9 Nu Al — 5) 55.0 54-4 53.4 | — 54.6 
100 mb — — — —59.8 | —51.9 || —54.5 | — 52.9. | — 53.3 — — 56.0 
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—14,5° C, respectively. As can be seen from 
our fig. 4 we have three peaks for the same level, 
the warmest of which is situated at —6,5° C. 
This discrepancy can be explained if we look 
at table I. The winter months chosen by 
McIntyre are definitely warmer than the ones 
chosen by us. This difference of the mean 
temperatures is brought about by the lower 
frequency of warm Tropical air as well as 
by higher frequency of rather cold Polar air. 

The other two maxima correspond to old 


Polar air and fresh Polar air. The latter of these . 


two air masses has come down more directly 
from the cold air reservoir in the north, being 
recently transformed Arctic air. 

This lack of warm Tropical air seems to 
have almost disappeared at the soo mb level, 
where we get a maximum at —22,5° C, com- 
pared with McIntyre’s at —21,5° C. The higher 
frequency of cold Polar air in our case is 
clearly seen in fig. 4, which has an extra 
maximum at —36,5°C, whereas McIntyre 
has only one cold maximum at — 31,5°C. 

Anyhow there is no question about the 
existence in both cases of a strong baroclinic 
field. The warm air maximum in fig. 4 has 
increased compared to the distribution at 
700 mb at the expense of the cold air maxima. 

The single maximum distribution of the 
300 mb surface temperatures is established for 
both periods. In our case the temperature 
corresponding to maximum frequency is 
about 3°C colder than in the other case, and 
the distributions have a rather different shape. 

If we now continue higher up in the atmos- 
phere, we once more get the three maxima at 
the 200 mb level. If there were a sloping 
baroclinic field above this station between the 
200 and 150 mb surfaces, then the maximum 
corresponding to the warm (Polar) air should 
increase at the expense of the colder peak. 
As can be seen from fig. 4, this is indicated in 
the upper diagram. The question whether the 
maximum corresponding to Tropical air is 
real or not will be discussed below. There a 
discussion of the validity of the three maxima 
in fig. 4 also will be given. 

We have thus come to the conclusion that 
for this North American station as well as 
for the British stations there exists a sloping 
baroclinic field in the stratosphere. 

As to a station situated far to the south on 
the North American continent one would 
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Fig. 5. Winter temperature frequency polygons for 
Miami (26°N, 80° W). The dashed line gives the mean 
values for the different levels. 


expect that the sloping baroclinic field seldom 
exists. As an example McIntyre has chosen 
Santa Maria in California. The frequency 
distribution of this station (see his fig. 2) 
shows a very marked peak corresponding to 
Tropical air. The few cases with Polar air 
form a tail to the left in the frequency diagram, 
more distinct at the 700 mb level than at 
500 mb level as is to be expected. 

The present author has chosen the station 
Miami in Florida (see fig. 5), instead of Santa 
Maria in southern California. This change 
has been done because of the fact that here 
the Polar air very seldom reaches above the 
700 mb surface and the Tropical air thus 
dominates the frequency distribution. The 
one-maximum curve must be found at all 
levels if the reasoning of the first paragraphs 
is right. Fig. 5 gives the frequency distributions 
with very marked peaks and the disturbances 
come out only as a few lower values in the 
700, soo and 300 mb diagrams. Higher up in 
the 200 and 100 mb surfaces these disturbing 
values are naturally higher than the tempera- 
ture corresponding to the maximum frequency. 


IOoo 


When studying the temperature frequency 
distribution in the main isobaric surfaces for a 
single station, it is generally very difficult to 
decide whether the different maxima are real 
or not. The double maxima of the diagram 
for 170 mb in fig. 2 and the triple maxima of 
the diagram for 200 mb in fig. 4, for instance, 
may be the result of this special sampling. 
The populations from which these samples 
are drawn might just as well be normally 
distributed. When assuming that this is the 
case and investigating the deviation of the 
sample distributions by standard statistical 
methods (y?-methods at the 5 % level) one 
obtains no definite results. In some cases the 
hypothesis of normal distribution must be 
rejected, in others not. Furthermore it may be 
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so that the population is not normally distri- 
buted. 

Instead of making a loose assumption about 
the distribution and investigating the existence 
of a significant deviation the present author 
prefers to use an indirect proof. In Melntyre’s 
paper as well as in this paper quite a lot of 
examples are given and in all cases there are 
at least indications of the phenomena on 
which the whole reasoning in these papers 
are based. The probability that this is a false 


. result of the ramdom sampling must be 


considered small. 

My conclusion is thus that the (sloping) 
baroclinic field in the troposphere is shown 
(and indirectly proved) by McIntyre’s figures, 
and the corresponding field in the lower 
stratosphere by my figures. 
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Expectancy of Clear Days in the Central Zone of the Total 
Eclipse of the Sun June 30, 1954, in Sweden 


By F. LINDHOLM, Swedish Meteorological and Hydrological Institute 


(Manuscript received Dec. 9, 1952) 


Abstract 


The occurrence of the total solar eclipse in 1954, visible in southern Sweden, has raised the 
question regarding the cloudiness conditions which may be expected in this part of the country 
at the time of the eclipse. A detailed investigation of the average sunshine conditions at Swedish 
stations situated along the track of the totality has therefore been made. From this is concluded 
that the coastal regions are decidedly more favourable than the inland. A comparison between 
the western and eastern ceasts leads to the conclusion that the average cloudiness at 14h, the hour 
of observation nearest to the totality, is slightly higher on the West coast than at locations on 
northern Öland and southern Gotland. Regarding the number of clear sky days at different places 
the coasts are also more favourable than the inland, and of those the northern part of Öland 
seems to be considerably more favourable even than the outer part of the West coast. The 
statistics refer to observations at 14h June 20—July 10 in the years 1931— 1950. 


It is evident from the maps of average daily 
cloudiness and sunshine in Scandinavia 1887— 
1900 by H. E. HAMBERG (1908) as well as 
from recent maps based on observations 1901— 
1930 that the coasts and the sea off both the 
West- and the East coasts of Götaland along 
the track of the eclipse on June 30, 1954 
(GRÖNSTRAND, 1950) are more favourable, or 
statistically seen, less unfavourable for astro- 
nomic observations than the interior of the 
country. 

The reversed daily variation of the cloudiness 
at inland stations and at coastal stations de- 
pending upon orographic effects on cloud 
genesis over land tends to increase the average 
cloudiness (isonephs) gradient at the coasts 
at noon, when the totality of the eclipse 
occurs. This is seen in fig. ı which is based 
upon observations at 14% M.E.T. in June and 
July during the period 1901—1930 from 
which monthly mean values have been 
computed. 
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The corresponding distribution of the dura- 
tion of sunshine (isohels), given in fig. 2, 
shows similar features. The monthly mean 
values of sunshine in per cent of the possible 
amount between 13 and 14? S.T. (solar 
time) have been calculated from the values 
of cloudiness in fig. 1 by application of the 
regression function s = 36 + 0.745 C,. This 
equation was deduced from simultaneous 
records of duration of sunshine between 13" 
and 14" in per cent (s) and observations of 
cloudiness N in scale o—roo at 14" during the 
period 1927—1949 at the meteorological sta- 
tion Vinga Lighthouse on the Westcoast 
(C, is per cent of clear sky: C, = 100 — N). 

The statistics in table r has been prepared 
in order to establish more clearly whether 
the West- and East coasts of Gétaland differ 
significantly in their mean cloudiness and 
sunshine duration and therefore whether any 
of these two parts of the country would be 
preferable with regard to cloudiness conditions 
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Fig. 1. Isonephs of average cloudiness at 14h M.E.T. 


June—July 1901—1920. 


for astronomic observations. The cloud obser- 
vations at 145 during the period June 20— 
Juli 10 1931—1950 have been treated for Vä- 
deröbod, Hällö and Mäseskär on the Skager- 
rack coast and Västervik, Olands Norra Udde 
and Hoburg on the East coast of Sweden. All 
these stations are located near the central 
line of the eclipse. The whole material has 
been divided in two ten-years groups. 

As is seen from table 1 the mean cloudiness 
at the west coast stations and Hoburg on 
southern Gotland has increased in the later 10 


Fig. 2. Isohels of average relative duration of sunshine 
(94) 13h—14h S.T. June—July 1901—1930. 


years period 1941—1950 corresponding to 
both a reduced number of clear days and an 
increase of the number of overcast days. On 
the other hand Västervik and Olands Norra 
Udde show a marked decrease of the mean 
cloudiness during this period chiefly due to a 
decrease of the number of overcast days. The 
number of clear days has remained at a constant 
value at these two stations. 

Regarding the whole period 1931—1950 we 
find the lowest average cloudiness of the 6 
stations at Olands Norra Udde. The difference 


Table 1. Average cloudiness (N, scale o—10) and frequency of clear or almost clear days and of overcast 
or near overcast days at 14h M.E.T. 


X, = mean number of days N = 0—2, June 20—July Io 
%, = mean number of days N = 8—10, June 20—July Io 
4 1931—1940 1941—1950 1931—1950 
o 
Stations Be Lat 
N | ¥1 Je@ıl %2 left) N | 1 fe(x)] Xe Je(%a)} N | %1 Je@ıl Xe [e)%) 
Väderöbod| 11° 2’|58° 31°| 4.1| 7.7| 0.9| 4.9 | 0.5 | 5.5! 4.8| 0.5 | 7.6| 0.7| 4.8| 5.9| 0.5 | 6.3 | 0.5 
Hallé ....| 11° 13’| 58° 20° | 4.3 | 6.7 | 1.0] 4.5 | 0.6| 5.2| 4.9| 0.7 | 5.8] 0.7 | 4.7| 5.8| 0.6| 5.2| 0.5 
Mäseskär .| 11° 21’| 58° 6°) 4.8| 5.9| 0.9| 5.0| 0.6] 5.6] 5.4] 0.8| 7.8] 0.8] 5.2| 5.7 | 0.6| 5.7 | 0.6 
Västervik.| 16° 36°| 57° 47’| 5.7| 6.3| 0.9| 9.3 | 0.8| 4.8 | 6.4 | 0.9 | 5.7 | 0.6| 5.2 6.4| 0.6| 7.5 | 0.6 
Ol. n. udd.| 17° 6°| 57° 22’| 4.6| 9.7| 1.2] 7.5| 0.8| 4.2| 9.8| 1.3 | 5.8] 0.9 | 4.4| 9.7 | 0.9 | 6.7 | 0.6 
Hoburg ..[18° 9°| 56° 55’| 4.6] 7.8| 1.0] 6.6| 0.9] 5.4 | 5.3| 0.7 | 7.6| 1.0| 5.0] 6.6| 0.7 | 7.1| 0.6 
Visingsö ..| 14° 2071587 57) — || — | — | = | 6.6] 1.2 | 0.5 | 9.0] 2.4 — 
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between Olands Norra Udde and the other 
stations is much more pronounced for the 
number of favourable days with cloudiness 
Be >, 

The table shows a marked predominance 

of clear days also at Hoburg and Västervik; 
thus the eastern stations are more favourable 
than the western ones, Väderöbod, Hällö and 
Mäseskär. 
‚The frequency of clear days is 46% at 
Olands Norra Udde, about 30 % at Hoburg 
and Västervik; at Väderöbod, Hallé and 
Mäseskär clear days occur in 27 à 28 per cent 
of the total number of days in the period 
June 20—July 10. 

The frequency of overcast days, cloudiness 
8—10, is of the same order of magnitude as 
that of clear days. We find that overcast days 
occur in 32 to 36 per cent at the eastern sta- 
tions and 25 to 30 per cent at the stations on 
the West coast. 

From the calculated standard errors e (x) = 
= —— 2 , o being the standard deviation and 

Vn—ı 
n the number of years, we conclude that only 
the mean for clear days at Olands Norra Udde 
differs significantly from those of the other 
stations. If we divide the stations according to 
their means in three groups; Olands Norra 
Udde x = 9.7 + 0.9, the westernmost stations 
X = 5.8 + 0.6 and the easternmost stations 
except Olands Norra Udde x = 6.5 + 0.7 
clear days of a total number of 21 days, we 
obtain the following values of the differences 
of these mean values x, — x, and their standard 
errors € (X1—X2) = Ve: (x,) + €? (x2) (table 2). 

The differences between the mean at Olands 
Norra Udde and the means of the westernmost 
and also of the other easternmost stations are 
three times as great as the standard error of 
the differences. The excess of clear days at 
Olands Norra Udde is thus statistically signifi- 
cant. One cannot, however, exclude the 
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Table 2 
Stations X — X | E(X1 —%;) 
Ölands Norra Udde—western- 
COE CLEMO RE PAR once 3.9 1.08 
Olands Norra Udde—eastern- ’ 
IHOSUISTAON EE RE ET 62 arr 
Easternmost—westernmost 

SAONE PSM Ae rer 0.7 0.88 
Easternmost stations—Vi- 

SIN SOS dre hare ee 5.5 0.86 


possibility that this predominance partly or 
completely may be due to the personal equa- 
tion in the daily estimation of the cloud 
amount at this station. 


Nothing can be said about the reality of 
the difference of the means between the 
eastern stations Västervik and Hoburg on 
one side and the western station on the other 
side, the standard errors being as large as the 
difference itself. 


Professor T. Bergeron has suggested that 
the convection over the island of Visingsö in 
the lake of Vättern ought to show similar 
characteristics as over the islands of Oland 
and Gotland. Ten years of observations at 
Visingsô 1941—1950 do not support this 
idea. The frequency of clear days at Visingsö 
only amounts to 5 % which should be com- 
pared with 30 % at the East coast. The average 
cloudiness of 6.6 is high and the frequency of 
cloudy days amounts to 43 %. 


Records of duration of sunshine using the 
Campbell Stokes autographs give a more 
objective estimation of the cloudiness. Un- 
fortunately there is only one station with 
records available in the interesting eastern 
part of the track, viz. Olands Södra Udde to 
the south of the central line and these observa- 
tions cover only the last ten years. The values 
taken from the records at Olands Södra Udde 
and at Vinga at the West coast are summarized 


Table 3. Average duration of sunshine and frequency of clear days, duration 0.8—1.0 hours and of 
overcast days, duration 0.2—0.0 hours between 13h—14h S.T. June 20—July 10 


Average Mean number of Mean number of 
Stations oe Lat. duration days duration days duration 
We 0.8—1.0 hours 0.2—0.0 hours 
Vierer 030% 38% 73.8 | 14.2 + 0.9 3.9 + 0.5 
Ölands Södra Udde ..| 16° 24’ lo 1% 76.5 147-1007 3.2 + 0.5 
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in table 3 for the same 21 days period of nine 
years between 1940 and 1952. 

_ The means 74 % at Vinga and 77% at 
Olands Södra Udde in the table agree well 
with the values on the map (fig. 2) showing 
the relative duration of sunshine between 
135142 S.T. calculated for the period 
I90I—1930. 

The mean number of days with sunshine 
for 0.8—1.0 hours is somewhat greater at 
Olands Södra Udde than at Vinga as well as 
the average duration; the frequencies being 
7o per cent of the total number of days at 
Olands Södra Udde against 67.8 per cent at 
Vinga. The mean number of overcast days, 
0.0—0.2 hours of sunshine, is also somewhat 
greater at Vinga than at the station in the 
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eastern part. But, the differences between the 
means of these two stations are too small com- 
pared with their standard error for securing 
a statistically certainty both with regard to 
clear and overcast days. 
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SHORTER, CONTRIBUTION 


Note on the Classic Stability Problem 


By VICTOR P. STARR 
Massachusetts Institute of Technology 


(Manuscript received November 14, 1952) 


It has come to the writer’s attention (private 
communication) that a question has arisen con- 
cerning his comments (STARR 1950) in regard to 
the stability of zonal flow. In terms of the problem 
as outlined previously, and in terms of the same 
symbolism, the question relates to the statement (1) 
that a monotone distribution of the vorticity £ with 
respect to yo implies stability of zonal motion. The 
point at issue is whether an incipient small disturb- 
ance of otherwise arbitrary nature imposed upon 
such zonal flow can so alter the circumstances of the 
problem as to render the previous conclusion 
unjustified. The purpose of this second note is 
to show that such small arbitrary perturbations 
cannot lead to instability and hence need not be 
considered explicitly in the case of a sensibly mon- 
tone vorticity distribution. 

A disturbed vorticity pattern having one local 
extreme of vorticity superposed on a monotone 
distribution is shown in the figure. The entire zonal 
belt A may now be divided as shown into three 
regions Ar, Ar, A3, the second being simply con- 
nected while the other two are doubly connected. 
Obviously the vorticity lines in A» do not encircle 
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Figure T. See text. 
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the zonal belt. It is to be observed that yo which i 
the mean value of y for the lines of constant ¢ 
changes value abruptly from b to c as one passes 
from region Ar to region A; along the common 
boundary. We may now form the integral /\ over 
the regions Ar and A; where the £ lines encircle 
the whole belt. This quantity is then related to the 
invariant I for the entire region A as follows: 


b d 
dé dé 
— WP) e2d}d — D eddy, = 
ei io Yo 
= Pepar + “ ff reaa—5 I + invariant. 
2 As 


This equation may be verified by the methods 
previously used. The two integrals on the left 
could be written as a single one with proper inter- 
pretations as regards the discontinuity in yo between 
A; and 43. The first term on the right is the line 
integral of Cy? around the region A: with respect 
to A, while the second involves the area integral of 
yf over this same region. Here according to our 
assumptions d&/dyo is of the same algebraic sign in 
the regions Ax and Az. 

If we are to deal with an incipient small disturb- 
ance in the initial state, it is essential that each of 
of the four integrals mentioned above should at 
most be infinitesimals as compared with 21/r. 
Otherwise the vorticity distribution would be 
sensibly disturbed already at the initial moment, 
and the problem would not be one relating to the 
stability of sensibly zonal flow. Generally speaking 
this means that the area Az should be small initially 
(and therefore also later) as compared with A, and 
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that & in regions Ar and A; should be infinitesimal 
as compared with d—a. 

Since d&/dt= o and the range of y is limited, the 
first two integrals on the right can never become 
finite in the subsequent motions of the fluid so 
that as a result, due to the invariance of I, the inte- 
grals on the left can never become finite. It thus 
follows that sensible departures from a zonal dis- 
tribution of vorticity cannot develop later. Several 
further points should be noted. 

(a) Any finite number of closed regions similar 
to A: may be included by generalizing the treat- 
ment. Likewise some of these regions may be 
contiguous with the latitude walls without producing 
a change in the result. In the end the treatment can 
include all possible topological configurations of 
the incipient disturbances, bringing them under 
the scope of the arguments set forth above. 

(b) In all probability other and perhaps better 
arrangements of the argument given may be made. 
However in any case the crux of the problem lies 
in the fact that the development of finite anomalies 
in the vorticity distribution would require a finite 
change in the average value of y for all the vortices 
present. This is precluded by the invariance of I. 
On the other hand, if there were substantial fractions 
of the area A over which the vorticity distribution 
is not monotone, the constraint imposed by the 
invariance of I would no longer prevent the develop- 
ment of such anomalies, since contributions from 
different regions to such a change of the mean 
value of y for the vortices may then cancel. 

(c) In most discussions of stability the question 
ultimately arises as to what is to be regarded as a 
suitable definition of this term. For the present 
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problem such questions are to a large part due to 
the fact that we do not have sufficiently general 
explicit solutions of the vorticity equation at our 
disposal. Thus if a sufficiently complete and accurate 
discription of the physical behavior of hydrodynamic 
systems such as the one here treated were available, 
the question would be primarily one of terminology. 

In may be noted, however, that the treatment 
given contains one particular instance where the 
meaning of stability is considerably simpler, namely 
the case of rotation as a solid (dö/dyo = constant). 
This may serve as an illustration of the essential 
properties of a stable system, although from view- 
points other than the present one it possesses special 
characteristics differing from those of stable systems 
in general. 

(d) It is interesting to speculate as to whether it is 
permissible to simplify vorticity distributions such 
as the one shown in the figure by the use of what 
may be termed vorticiry CUTS. By this is meant the 
artificial introduction of lines along which the 
vorticity lines are extremely closely packed but 
along which the vorticity remains finite in intensity. 
By the appropriate use of such devices all possible 
vorticity distributions in the present problem may 
be reduced to the class of cases for which all the 
vorticity lines extend completely around the zonal 
belt, without altering the corresponding motions 
significantly. Conclusions concerning stability could 
then be reached more directly. 
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Starr, V. P., 1950: Note on recent study of stability 
by R. Fjortoft. Tellus, 2, pp. 321—322. 
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Letters to 


toth December, 1952. 


Dear Sir, 


In a paper published in this journal CHARNEY 
and ELIASSEN (1949) have discussed the influence 
of a mountain barrier on a westerly current 
flowing over it. In their one-dimensional treatment 
of the problem the effect of horizontal convergence 
in the flow in the lee of the mountains is shown 
to be a negative pressure (contour-height) tendency 
over a limited distance downstream from the peak. 
A series of actual tendency computations, in the 
soo mb pressure surface, made by Borm and 
CHARNEY (1951) on the basis of the two-dimen- 
sional, barotropic non-divergent vorticity equation 
gave a. o. an apparent verification of the existence 
of such a convergence area in the lee of the 
Rocky Mountains. In a preliminary note (see 
Bou and H. NEWTON, 1951) the present writer 
purported to have verified that topopraphic effect 
in a case of similar large-scale flow over North 
America, through computation on a series of 
five-day mean soo mb contour charts. 

Meanwhile a check of signs involved in the 
interpretation of the difference between the com- 
puted and observed tendency has shown that an 
identical oversight has occurred in both the latter 
authors’ and the present writer’s study. Correc- 
ting for the mistake the respective area, in either 
case in the lee of the Rocky Mountains and 
downstream from a broad westerly current crossing 
it, would have to be regarded as the seat of 
horizontal divergence in the ligth of the limited 
interpretation of that difference in tendency. It 
must then be assumed that during those two-week 
periods treated the topographic effect, if any, 
was overshadowed by other processes such as a 
soleneoidal or frictional source or else a large- 
scale horizontal eddy-transfer of vorticity in the 
region. 

In view of the general importance which is 
attached to the significance of semi-permanent 
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fields of divergence and convergence, it may be 
indicated to give here what is contended to be 
the right interpretation of the numerical findings 
by Borm and CHARNEY. | 

The quasi geostrophic vorticity equation for 
the purpose of the present discussion has the 
(reduced) form 


où 
2 V-Vn=—y7 div, V 
where the usual notation is used, © being the 
geostrophic relative vorticity and n the absolute 
vorticity (7=¢€-+f). Since in the notation of 
the authors J¢/dt = gf— \72 x, and consistently 
with t—V-VYn=gf"1W?x, the above 
equation becomes 
V7? (x, — x) = = div, V 
8 


Applying the authors’ regression equation 
want b of X= = —b+ x, (t—a), 
it follows 


nu div V — 72 


\/2b (x3 (:—a)] 


In considering the distribution of the constant b 
in relation to divergence, the regression coeffi- 
cient a has to be put equal to unity, whence 


WE b 1 div, V 


The quantity fn =f? + f¢ may be approximated 
by f? so that fng! = A, a positive quantity. 
Therefore 

NAIL D = A div, V 


In Fig. 7 of the paper by Bouin and CHARNEY 
the factor b is negative over the area in question 
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and has a pronounced minimum. The Laplacian 
in the vicinity of that minimum is certainly 
positive as it is also over a large area surrounding 
the minimum. Thus the divergence is positive 
in that area. 

Although the authors do not explicitly identify 
that area of negative b with a region av con- 
vergence (div, V<o), they draw attention to 
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the resemblance with the negative tendency field 
corresponding to a topographically produced con- 
vergence area (Fig. 8). Naturally they did not 
imply resemblance in position and shape only, 
as is clear from the context. 


Yours sincerely, 


F. A. BERSON 
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REPLY 


Dear Sir, 

We should like to thank Dr. Berson for 
pointing out an inconsistancy in our article on 
the influence of the mountains on the flow pattern 
at the 500 mb surface. Checking the original 
computations once more we have found that 


the inconsistancy is explained by a drafting error 
in the sign of the isopleth labels of figure 7. 
Our conclusions therefore remain as before. 


Sincerely yours 


Bert Bouin, JULE CHARNEY 
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